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Introduction 


This  is  our  first  annual  progress  report  representing  work  done  on  the  grant  in  the  past  year. 
Some  of  the  data  has  already  been  reported  in  one  of  our  three  Quarterly  Reports  submitted 
since  the  grant  was  awarded.  This  report  is  more  comprehensive  and  includes  updates  from 
previous  reports  as  well  as  information  from  other  investigators  not  previously  reported.  We  will 
relate  our  work  progress  to  the  milestones  outlined  in  the  grant  proposal. 

We  hypothesized  that  the  differential  susceptibility  to  DoD-priority  pathogens  is  the  result  of  host 
genetic  variability,  and  that  these  discreet  loci  and/or  gene  pathways  can  be  identified  using 
BXD  recombinant  inbred  mice.  The  DoD  priority  pathogens  are  naturally  occurring  endemic 
diseases,  emerging  infectious  diseases,  and  potential  biowarfare  agents:  Leishmania  major, 
Burkholderia  pseudomallei,  Severe  Acute  Respiratory  Syndrome  (SARS-CoV),  highly 
pathogenic  H5N1  Avian  Influenza  virus,  Francisella  tularensis  and  multidrug  resistant 
Acinetobacter  baumannii. 

We  outlined  a  three-step  sequence  for  studying  each  pathogen  with  the  following  specific  aims: 

Objective  1  (18  months):  To  identify  specific  phenotypic  differences  in  the  responses  of  the 
BXD  parental  strains  following  infection  with  each  of  the  DoD-priority  pathogens. 

Objective  2  (6  to  30  months):  To  identify  host  genetic  loci  and  pathways  that  correlate  with  the 
differential  susceptibility/resistance  phenotype(s)  of  the  parental  mouse  strains  to  the  DoD- 
prioirty  pathogens. 

Objective  3  (12  to  36  months):  To  define  and  validate  candidate  genes  and  gene  networks 
responsible  for  differential  susceptibility/resistance  phenotype(s)  of  the  parental  mouse  strains 
for  each  DoD-priority  pathogen. 

Each  project  is  performed  in  parallel  to  achieve  the  milestones  set  forth  in  Objectives  1-3.  As  we 
stated  in  our  proposal  the  infection  models  are  at  various  levels  of  development,  and  therefore 
the  timelines  for  the  objectives  differ  according  to  the  pathogen.  We  are  reporting  the  progress 
of  each  pathogen  model  in  the  one-year  period  since  the  grant  was  initiated.  In  year  1 ,  most  of 
the  work  performed  was  related  to  Objective  1  under  the  following  subheadings: 

Objective  la.  Infection  of  parental  strains  with  DoD-priority  pathogens  and  sample 
collection. 

We  will  report  on  ongoing  studies  for  the  different  pathogens  related  to  this  objective. 

Objective  1b.  Infection  of  primary  cells  with  DoD-priority  pathogens  and  sample 
collection. 

We  will  report  on  ongoing  studies  for  the  different  pathogens  related  to  this  objective 

The  format  of  this  report  will  present  the  data  of  the  ongoing  projects  for  the  different  pathogens, 
as  well  as  reports  from  the  animal  core  and  bioinformatics  core. 
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Key  research  accomplishments,  reportable  outcomes,  conclusion,  references 
1 .  Projects  with  DoD-priority  Bacterial  Pathogens 

Task  1:  Infection  of  parental  strains  with  DoD-priority  pathogens  and  sample  collection  (Year  1: 
12  months) 


1 .1  Studies  with  Francisella  tularensis  (FT) 


1.1.1  Overview 

We  have  been  using  the  live  vaccine  strain  (LVS)  to 
identify  phenotypic  differences  in  the  responsiveness  of  the 
two  parental  strains  to  Francisella  infection.  We  have 
identified  several  phenotypic  differences  with  the  LVS 
strain  that  could  be  used  as  readouts  for  BxD  analyses. 

Through  a  series  of  studies  in  which  immune  cell  influx  into 
the  lungs  of  mice  was  determined  following  infection  with 
FT  we  have  found  that  FT  delays  neutrophil  recruitment 
into  the  lung  by  several  days  (data  not  shown).  We  have 
also  found  that  there  is  a  significant  difference  in  the 
number  of  neutrophils  that  are  recruited  into  the  lungs  of  B6  vs.  D2  mice  (Figure  1)  that 
correlates  well  with  the  differences  observed  in  the  bacterial  burdens  measured  in  spleens  and 
livers  of  infected  mice  (Figure  2).  Interestingly,  there  appears  to  be  an  inverse  correlation 
between  the  ability  of  D2  mice  to  recruit  neutrophils  (which  in  part  is  responsible  for  reduced 
bacterial  burden  in  the  spleen  and  liver)  and  susceptibility  to  infection.  In  our  hands,  D2  mice 
appear  to  be  marginally  more  sensitive  to  FTLVS  infection,  especially  when  dosed  with  5  x  104 
CFU  via  the  intranasal  route  (Figure  3). 


Figure  1  B6  or  D2  mice  (5/group)  were  challenged 
with  5  x  10e4  CFU  FTLVS  via  the  intranasal  route 
Five  days  later,  brochioalveolar  lavage  was 
performed  and  the  neutrophil  populations  were 
detected  via  flow  cytometric  analysis. 


Bacterial  Load  in  organs  at  5  days  p.i. 
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Figure  2.  C57BI6J  or  DBA5J  mice  (5tjroup)  were  challenged  with  5  x  10e4 
CFU  of  FTLVS  via  the  intanasal  route  five  days  post  infection,  the  indicated 
organs  weie  collected  via  sterile  dissection,  homogenized,  and  then  the  FT 
burden  was  determined  by  dilutuion  plating  on  Mueller- Hinton  agar-  plates 
"=  p.0.05;  n.s.=  not  significant. 


A 


mica 


B 


D2  mice 


C57BL*\M  ’  Cr.»H-V  <»<••? 


♦-  LwiWi'/'i  -*?- 


a 


10U 


50 


0  - - - 

0  5  10  IS 

Hnys  po«l  Infnr.bon 


Figure  3.  Female  age-matched  C57BL6  and 
DBA2J  mice  were  infected  intranasally  with  FTLVS 
at  the  indicated  doses.  Survival  was  recorded  daily 
over  a  period  of  20  days. 
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We  have  also  performed  several  preliminary  studies  in  an  attempt  to  identify  potential  readouts 
of  differences  between  B6  vs.  D2  responsiveness  to  FT  Schu  S4.  We  have  found  that  D2  mice 
express  significantly  higher  amounts  of  TNF-a  and  SOCS3  than  B6  mice,  while  B6  mice 
express  larger  quantities  of  both  IFN-  y  and  MIP-2  than  do  D2  mice  (Figure  4). 


Figure  4.  Age  matched  female  mice  were  infected  with  1X106 
cfu  of  SchuS4  via  the  intraperitoneal  route  and  sacrificed  at  48 
hours  post  infection.  Total  RNA  was  extracted  from  liver  tissue 
and  analyzed  for  gene  expression.  Relative  levels  of  A)  TNF 
alpha,  B)  IFN  gamma.  C)  SOCS3  and  D)  MIP2  mRNA  was 
analyzed  by  real-time  PCR,  normalized  to  (5-actin  levels  and 
expressed  as  fold  increase  compared  to  RNA  from  PBS 
treated  group.  Data  is  derived  from  a  single  representative 
experiment  and  presented  as  mean  +  SEM. 


Figure  5.  FT  SchuS4  was  transformed  with  pXBI  78-LUX  and  cultured 
overnight  on  modilied  Mueller-Hinton  agar.  The  image  was  captured  using  an 
IV  IS  Spectrum  whole  animal  imaging  system. 


We  have  also  made  progress  in  our  efforts  to  monitor  FT  infection  via  whole  animal  IVIS 
imaging.  First,  we  have  verified  that  the  pXB  178-LUX  reporter  vector  can  be  used  to  drive 
photon  production  in  FT  Schu  S4  (Figure  5).  In  addition,  using  FTLVS  transformed  with  the 
pXB  178-LUX  reporter  vector,  we  are  able  to  visualize  dissemination  of  FT  infection  in  real  time 
following  intranasal  challenge  (Figure  6).  Each  of  the  mice  in  this  study  were  challenged  with 
approximately  1  x  106  FTLVS  in  either  5pl,  lOpI,  25pl,  or  50pl/nare.  The  mice  that  received  their 
challenge  dose  in  only  5pl/nare  each  experienced  a  non-lethal  upper  airway/superficial  cervical 
node/salivary  gland  infection  that  peaked  in  intensity  between  48  and  72  hours  post-infection 
and  never  disseminated  to  the  lungs,  spleen  or  liver.  All  of  the  other  animals  in  the  study 
experienced  a  fully  disseminated  infection  that  resulted  in  mortality  2-3  days  after  bacteria  were 
detectable  in  the  liver  via  IVIS  imaging.  We  look  forward  to  evaluating  the  ability  of  FT  Schu  S4 
to  disseminate  from  the  upper  airways  following  intranasal  dosing  in  5pl/nare  (which  fails  to 
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3  Hrs  Post-challenge 


Figure  6  BALB/c  mice  (3/group)  were  challenged  with  1  x  10e6  CFU  of  FTLVS  bearing  the  pXBI  78-LUX  reporter  plasmid  via  the  intranasal 
route  in  either  5pl,  lOpI,  25pl,  or  50pl/nare.  Dissemination  of  FTLVS  was  then  monitored  using  an  IVIS  Spectrum  whole  animal  imaging  system 
Images  were  collected  at  the  indicated  timepoints  post-infection.  The  images  collected  from  all  mice  were  normalizied  to  reflect  photons  per 
second  per  cmA2/sr. 

deliver  the  bacteria  efficiently  to  the  lungs).  These  studies  will  be  initiated  as  soon  as  we  are 
able  to  begin  FT  work  in  the  RBL.  All  indications  are  that  the  RBL  will  be  available  for  these 
studies  by  the  end  of  July  at  the  latest. 


We  have  also  generated  a  number  of  attenuated  mutant  strains  of  FT  LVS  and  have 
characterized  a  couple  of  them  in  some  detail.  We  are  currently  working  to  generate  clean 
deletions  of  the  same  genes  (if  possible)  in  the  Schu  S4  background.  If  the  mutations  are 
attenuating  in  the  Schu  S4  strain,  we  will  be  able  to  use  these  strains  for  additional  BXD  studies 
in  FT.  Their  reduced  virulence  will  likely  make  it  possible  for  us  to  identify  additional  readouts  of 
resistance/susceptibility  differences  between  the  parental  mouse  strains  that  would  otherwise 
be  masked  by  the  extreme  virulence  of  the  wild-type  Schu  S4  strain.  Any  new  readouts  will 
increase  our  ability  to  identify  host  genetic  elements  that  correlate  with  susceptibility  to 
tularemia. 


1.1.2  Reportable  outcomes 


Publications 


1.  Clinton,  S.  R.,  J.  E.  Bina,  T.  P.  Hatch,  M.  A.  Whitt,  and  M.  A.  Miller.  2010.  Binding  and 
activation  of  host  plasminogen  on  the  surface  of  Francisella  tularensis.  BMC  Microbiol 
10:76. 
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1.2  Studies  with  Burkholderia  pseudomallei. 


We  have  not  performed  any  experiments  relating  to  the  Burkholderia  component  of  the  project. 
We  have  received  approval  from  the  CDC  to  add  Burkholderia  to  our  select  agent  program,  and 
we  are  in  the  process  of  acquiring  the  bacterium.  Our  efforts  have  been  derailed  to  date  due  to 
a  failure  of  UTHSC  to  reach  a  material  transfer  agreement  with  BEI  Resources.  We  expect  this 
to  be  rectified  within  the  next  month.  As  soon  as  we  acquire  the  strain(s)  we  will  begin  studies  to 
define  resistance/susceptibility  readouts  in  the  parental  strains,  and  then  we  will  move  into  BXD 
stains  to  identify  host  genetic  elements  that  correlate  with  strain-specific  differences  in  disease 
phenotypes. 
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1 .3  Acinetobacter  baumanii  project. 

Our  initial  goals  for  the  in  vivo  Acinetobacter  work  will  include  defining  the  baseline  for  challenge 
dosing  in  the  two  parental  mouse  strains.  Ideally,  we  would  identify  a  challenge  dose  that  is 
lethal  to  one  of  the  parental  strains,  but  not  lethal  to  the  other  strain.  Our  initial  studies  have 
suggested  that  survival  as  a  readout  for  resistance/susceptibility  differences  of 
immunocompetent  mice  will  not  be  possible  with  this  pathogen.  We  are  currently  attempting  to 
define  more-subtle  resistance/susceptibility  readouts  that  can  be  used  with  this  pathogen  (influx 
of  immune  cell  types  into  infected  lungs,  cytokine/chemokine  production  profiles,  innate 
signaling  differences,  etc).  Because  Acinetobacter- related  disease  is  typically  observed  as  a 
nosocomial  infection  of  immunologically  suppressed  patients,  we  have  also  begun  to  evaluate 
methods  of  immunosuppresion  for  the  mouse  model  (cyclophosphamide  treatment  or  gamma 
irradiation)  prior  to  infection.  Immunosuppression  prior  to  challenge  may  enable  us  to  more 
clearly  define  host  genetic  elements  that  relate  to  differences  in  resistance/susceptibility 
observed  between  BXD  mouse  strains. 

We  have  generated  a  luminescence  reporter  vector  for  Acinetobacter  and  are  now  able  to 
visualize  the  bacterium  in  vivo  via  IVIS  imaging.  In  fact,  we  are  currently  using  Acinetobacter- 
plux  to  develop  our  procedures  for  using  our  newly  acquired  aerosol  delivery  device.  These 
studies  will  help  us  to  define  ideal  inocula  concentrations  vs.  exposure  durations  for  reliable 
delivery  of  desired  CFU  to  the  lungs  of  mice. 

We  have  also  initiated  the  in  vitro  studies  with  Acinetobacter.  In  order  to  identify  specific 
phenotypic  differences  in  the  response  of  the  BXD  parental  strains  to  infection  with  the  DOD- 
priority  pathogen  Acinetobacter  baumanii  we  analyzed  the  innate  immune  response  of  bone 
marrow-derived  dendritic  cells  of  C57BL/6J  and  DBA/2J  strains  stimulated  in  vitro  with  live  or 
heat  killed  A.  baumanii.  As  shown  in  Figure  1 ,  A.  baumanii  induced  a  pattern  of  cytokines 
characteristic  of  Gram-negative  bacteria.  In  particular,  induction  of  the  chemokines  RANTES 
and  IP-10  suggests  that  A.  baumanii  activates  the  TLR4  pathway,  which  controls  the  expression 
of  RANTES  and  IP-10  through  the  TRIF  pathway.  Most  cytokine  and  chemokine  transcripts 
were  induced  to  the  same  extent  in  the  two  parental  strains.  The  mRNA  for  IL-1 2p40  however, 
was  induced  more  efficiently  in  the  C57BL/6J  cells  than  in  DBA/2J.  This  effect  may  reflect  the 
known  polarization  of  the  C57BL/6J  strain  toward  the  TH1  lineage. 

In  order  to  analyze  whether  A.  baumanii  is  capable  of  activating  the  inflammasome  leading  to 
secretion  of  IL-1  (3,  mouse  bone  marrow-derived  dendritic  cells  were  stimulated  with  A.  baumanii 
for  18  hours  and  IL-1 13  was  measured  in  supernatants  by  ELISA  (figure  2).  These  experiments 
revealed  that  A.  baumanii  is  a  poor  stimulator  of  the  inflammasome.  This  result  is  consistent 
with  the  notion  that  the  inflammasome  is  primarily  activated  by  intracellular  bacteria  or  bacteria 
that  possess  type  three  or  type  four  secretion  systems.  Future  experiments  will  determine 
whether  inflammasome  activation  occurs  in  vivo  during  A.  baumanii  infection. 
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Figure  1.  Comparison  of  innate  immune  response  to  Acinetobacter  baumanii  in  C57BL/6j 
and  DBA/2J  bone  marrow-derived  DC.  BM-DC  were  stimulated  for  four  hours  with  live  or 
heat-killed  (HK)  Acinetobacter  (A. b.)( moi  10)  or  for  3  hours  or  with  Listeria  monocytogenes 
(L.m.) or  Salmonella  thyphimurium  ( S.t .).  Total  RNA  was  extracted  from  cells  and  expression  of 
a  panel  of  cytokines  (A)  and  chemokines  (B)  was  analyzed  by  RNase  protection  assay. 
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Figure  2.  A.  baumanii  is  a  poor  stimulator  of  the 
inflammasome.  BMDC  were  stimulated  with  A.  baumanii 
(A.b.)  at  the  indicated  MOI  for  18  hours  or  with  Listeria 
monocytogenes  (L.m.)  and  IL-lb  was  measured  in 
supernatants  by  ELISA. 
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2.  Influenza  virus  (H5N1)  project 

2.1  Conclusions  reported  in  Quarterly  report  of  December  2009 

•  Inbred  mouse  strains  vary  in  their  response  to  inoculation  with  highly  pathogenic  H5N1 
influenza  virus  infection.  Several  strains  are  very  susceptible  with  infection  ultimately 
resulting  in  death,  while  other  mouse  strains  overcome  and  survive  the  infection.  The 
observed  differences  in  pathogenicity  are  caused  by  the  genetic  differences  between  the 
mouse  strains. 

•  Increased  susceptibility  to  highly  pathogenic  H5N1  virus  is  associated  with  higher  viral  loads 
in  the  susceptible  animals  as  well  as  increased  production  of  pro-inflammatory  cytokines. 
Both  items  are  hallmarks  of  H5N1  infections  in  humans. 

•  Recombinant  Inbred  mouse  strains,  derived  from  a  susceptible  and  a  resistant  mouse  strain 
identified  multiple  genetic  loci  associated  with  resistance  to  H5N1  pathology.  These  loci  will 
be  referred  to  as  QTL  for  influenza  virus  resistance  (Qivr). 

•  RNA  expression  analysis  in  lungs  of  H5N1  infected  susceptible  and  resistant  mice  also 
found  increased  pro-inflammatory  cytokine  mRNA  expression  in  the  susceptible  strains.  In 
addition,  RNA  profiling  identified  several  candidate  genes,  located  within  the  Qivr,  linked  to 
increased  resistance  following  infection. 

2.2  Research  Accomplishments 

Since  the  last  quarterly  progress  report  we  have  made  significant  progress  with  the  congenic 
mouse  lines.  These  lines  are  genetically  almost  identical  to  C57BL/6  (resistant)  mice  except  for 
a  small  region  that  is  derived  from  DBA/2.  These  regions  correspond  to  the  previously  identified 
loci  associated  with  susceptibility  to  H5N1  influenza  virus  infection.  Three  lines  were  obtained 
from  various  sources  and 
expanded  at  our  facility.  The 
B6.D2.1 1 D  congenic  line  was 
homozygous  upon  arrival  and 
after  expansion  was  tested  for 
its  susceptibility  to  A/Puerto 
Rico/8/34,  a  mouse  adapted 
H1N1  influenza  A  virus. 

Although  we  need  to  do  larger 
numbers,  the  data  suggests 
that  there  are  no  large 
difference  in  weight-loss  and 
mortality  of  these  mice 
compared  to  their  wild-type 
controls  (Figure  1).  Future 
experiments  will  also 
determine  the  susceptibility  to 
highly  pathogenic  H5N1  virus. 

Two  other  congenic  mouse  lines,  B6.D2.7C  and  B6.D2.17D,  arrived  at  our  facility  as 
heterozygous  mice.  After  extensive  breeding  and  genotyping  we  have  been  able  to  generate 
several  breeding  pairs  of  homozygous  B6.D2.17D  mice.  The  B6.D2.7C  mice  have  been  difficult 
to  breed  and  therefore  we  have  not  been  able  to  generate  homozygous  mice.  Hopefully  in  the 
next  progress  report  on  H5N1  influenza  we  are  able  to  report  on  these  mice  and  show 
significant  differences  in  disease  progression  for  one  or  both  of  these  congenic  lines. 
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Figure  1:  Survival  curve  after  infection  with  different  doses  (EID50) 
of  A/Puerto  Rico/8/34  virus. 
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Because  of  the  extreme  susceptibility  of  DBA/2  mice  to  highly  pathogenic  H5N1  virus,  we 
hypothesized  that  other  influenza  A  viruses  could  also  induce  severe  disease  in  these  mice.  A 
set  of  18  different  human  or  avian  influenza  A  virus  isolates  belonging  to  8  different  subtypes  of 
influenza  virus  were  pathotyped  in  DBA/2  and  C57BL/6  mice.  We  found  that  more  than  70%  of 
them  were  pathogenic  for  DBA/2  mice  while  only  two  viruses  were  pathogenic  for  C57BL/6  mice 
(Table  1). 


Table  1:  Mortality  after  intranasal  inoculation  of  DBA/2  and  C57BL/6  mice  with  various  human  and  avian 
influenza  A  virus  isolates 

Influenza  A  virus  isolate 

Subtype 

DBA/2 

106EIDso  104  EIDSo 

C57BL/6 

io6eid5(I 

A/Memphis/3/2008  1 

H1N1 

100 

88 

0 

A/California/4/2009  2 

H1N1 

100 

100 

60 

A/mallard/ Alberta/7 9/2003 

H2N3 

33 

0 

ND4 

A/mallard/ Alberta/33/2004 

H2N4 

0 

0 

0 

A/pintail  duck/ Alberta/66/2005 

H4N1 

20 

0 

0 

A/mallard/  Alberta/ 147/2007 

H4N6 

0 

ND 

ND 

A/shorebird/Dela  ware/ 101/2004 

H5N7 

20 

0 

0 

A/ruddy  turnstone/Delaware/103/2007 

H5N9 

100 

25 

ND 

A/teal/Hong  Kong/W3 12/1 997 

H6N1 

100 

100 

40 

A/mallard/ Alberta/1 54/2003 

H6N5 

100 

40 

0 

A/shorebird/Delaware/22/2006 

H7N3 

100 

100 

0 

A/mallard/ Alberta/177/2004 

H7N9 

0 

0 

0 

A/quail/Hong  Kong/G  1/1997 

H9N2 

100 

0 

0 

A/duck/Hong  Kong/Y280/1997 

H9N2 

100 

20 

0 

A/mallard/ Alberta/1 62/2007 

H9N5 

0 

0 

ND 

A/mallard/ Alberta/22 1/2006 

H9N6 

0 

0 

ND 

A/blue-winged  teal/Alberta/27 1/2007 

H10N5 

100 

75 

0 

A/mallard/  Alberta/5 6/2004 

H10N7 

100 

75 

0 

1  Human  influenza  A  virus  isolate,  2  Human  2009  pandemic  H1N1  vims  isolate,  4  ND,  not  done. 

Because  2009  heralded  the  first  influenza  pandemic  of  the  21st  century,  we  also  pathotyped  a 
pandemic  H1N1  influenza  virus  strain  called  A/California/04/09.  The  DBA/2  mice  were 
extremely  susceptible  to  this  virus  with  a  LD50  of  30  infectious  particles,  while  only  a  few 
C57BL/6  mice  succumbed  to  infection  with  106  particles.  Combined  these  data  show  that  DBA/2 
mice  are  very  susceptible  to  most  natural  non-adapted  influenza  A  viruses,  and  therefore  these 
mice  can  be  used  as  a  new  small  animal  model  in  vaccine  or  therapeutic  efficacy  studies. 
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We  used  the  model  to  understand  one  of  the  fundamental  questions  of  the  2009  influenza 
pandemic;  why  are  the  elderly  protected  from  infection  with  the  pandemic  H1N1  virus?  In  vitro 
hemagglutination  inhibition  and  virus 
neutralization  assays  suggested  that  sera 
from  elderly  contain  cross-reactive 
antibodies;  however  nobody  had  shown 
that  these  antibodies  were  fully  functional 
in  vivo.  For  our  study  we  obtained  human 
sera  from  60-80  year  old  individuals  and 
tested  the  sera  for  the  presence  of  cross¬ 
reactive  antibodies  against  the  pandemic 
H1N1  virus.  Sera  that  contained  virus 
neutralizing  antibodies  were  pooled  and 
injected  intraperitoneally  prior  to  challenge 
with  a  lethal  dose  of  A/California/04/09.  As 
a  control  we  used  sera  from  aged  matched 
individuals  that  lacked  any  detectable 
cross-reactive  antibodies. 

As  shown  in  Figure  2,  the  cross-reactive 
antibodies  in  sera  from  unexposed 
individuals  protected  the  mice  from  excess 
weight-loss  and  mortality  compared  to 
negative  sera  or  PBS  controls.  Based  on 
these  findings  we  suggest  that  the  elderly 
are  partially  protected  against  the 
pandemic  virus  and  that  these  antibodies 
were  most  likely  derived  from  a  prior 
infection  with  a  related  H1N1  virus. 


*** 
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Figure  2  Survival  curve  (A)  and  weight  loss  (B)  after 
intranasal  challenge  with  a  lethal  dose  of 
A/California/4/09  pandemic  H1N1  virus. 
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2.3  Conclusions 

•  DBA/2  mice  are  highly  susceptible  to  many  non-adapted  isolates  of  influenza  A  virus 

•  DBA/2  mice  can  be  used  as  a  sensitive  small  animal  model  in  vaccine  and  therapeutic 
efficacy  studies  specific  for  influenza  virus 

•  Cross-reactive  antibodies  identified  in  sera  from  elderly  individuals  is  fully  functional  in  vivo 
and  protects  against  severe  morbidity  and  mortality 

2.4  Reportable  outcomes 

Publications 

1.  Boon  AC,  Debeauchamp  J,  Krauss  S,  Rubrum  A,  Webb  AD,  Webster  RG,  McElhaney  J, 
Webby  RJ.  Cross-reactive  neutralizing  antibodies  directed  against  pandemic  H1N1  2009 
virus  are  protective  in  a  highly  sensitive  DBA/2  influenza  mouse  model.  J  Virol.  2010;  in 
print. 

Presentations 

•  Lectures  at  University  of  Maryland,  College  Park,  Maryland  and  Washington  University, 
St  Louis,  Missouri. 

•  Presentation  at  the  BXD  World  meeting  in  Braunschweig,  Germany  2009 
Abstracts 

•  Abstract  at  the  BXD  World  meeting  in  Braunschweig,  Germany,  2009 

•  Abstract  for  the  Options  for  the  Control  of  Influenza  VII,  Hong  Kong,  2010 

Animal  models 

Developed  a  new  sensitive  small  animal  model  to  study  human  and  avian  influenza  A  viruses 
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3. 


Leishmania  major  project 


3.1  Overview 

We  set  out  in  year  1  of  this 
study  to  determine  the 
feasibility  of  using  the  BXD 
resource  to  study  the  host 
response  to  Leishmania  major 
( Lm ),  a  protozoal  parasite  that 
causes  a  disease  in  mice 
similar  to  cutaneous 
leishmaniasis  in  humans.  The 
murine  model  of  Leishmania 
infection  has  been  very  well 
characterized.  For  our  studies 
purified  metacyclic 
promastigote  stage  Lm 
(WHOM/IR/-/173)  was  injected 
at  several  different  doses  into 
both  hind  footpads  (ifp)  of 
cohorts  of  mice  (aged  60±10 
days).  Just  prior  to  and  at 
weekly  intervals  post-infection, 
we  followed  the  clinical 
infection  course  by  using  a 
metric  caliper  to  monitor 
changes  in  the  diameter  of 
both  hind  footpads.  In  addition, 
at  regular  intervals  we 
recovered  ~50  pi  of  orbital 
sinus  blood  for  analysis  of 
serum  Ig  isotype  and 
circulating  cytokine  levels.  At 
selected  intermediate  times 
and  at  the  endpoint  of 
experiments  (~8-1 1  weeks), 
animals  were  sacrificed  and 
analyzed  as  follows.  Local  and 
disseminated  parasite  burdens 
were  determined  by  scoring 
cultures  of  serially  diluted 
spleen  and  footpad 
homogenates  for  outgrowth  of 
parasites.  Global  gene 
expression  profiles  will  be 
determined  at  a  later  date  if 
additional  funding  is  approved. 
This  will  utilize  RNA  recovered 
from  snap-frozen  draining 
popliteal  lymph  nodes  stored  in 


170 


<s> 


weeks  postinfection 


B 


weeks  postinfection 


Figure  1.  Pathology  of  L.  major  in  C57BL/6  (open  circles)  and  DBA/2  (solid  squares). 
Mice  were  infected  subcutaneously  in  hind  footpads  with  2  x  106  parasites.  (A)  Shown  is 
weekly  footpad  thickness  as  measured  by  metric  caliper.  Error  bars  indicate  standard 
error.  (B)  Shown  are  parasites  recovered  from  the  hind  feet  of  infected  mice  as 
determined  by  serial  dilution  of  feet  homogenates  cultured  on  blood  agar  plates  for  7-10 
days.  (C)  Shown  are  total  mononuclear  cell  numbers  per  draining  popliteal  lymph  node, 
as  determined  by  trypan  blue  exclusion  of  single  cell  suspensions. 


15 


our  tissue  bank.  Cell  suspensions  made  from  popliteal  lymph  node  were  used  to  determine 
antigen-specific  recall  cytokine  responses  to  soluble  Lm  antigen  (SLA).  As  shown  in  Figure  1A 
the  data  we  obtained  from  clinical  course  we  observed  was  consistent  with  published  reports 
with  DBA/2  being  less  resistant  than  C57BL/6.  Despite  the  differential  clinical  course,  parasite 
burdens  were  similar  between  the  BXD  parental  strains  (Figure  IB).  These  results  suggested 
that  differences  in  the  inflammatory  response  to  Lm  might  be  responsible  for  the  contrasting 
clinical  outcomes.  In  support  of  this  interpretation,  draining  popliteal  lymph  node  cellularity  was 
significantly  higher  in  DBA/2  versus  C57BL/6  (Figure  1C). 

To  investigate  the  nature  of  the  differential  inflammation  exhibited  by  DBA/2  and 
C57BL/6  in  response  to  Lm  infection,  we  performed  FACS  analysis  of  draining  popliteal  lymph 
node  cells.  We  found  that  the  number  of  MHC  class  II  positive  cells  (dendritic  cells  and 
macrophages)  was  significantly  elevated  in  DBA/2  versus  C57BL/6  (Figure  2A).  By  contrast,  the 
number  of  CD4  positive  cells  (T  helper  and  NKT  cells)  was  similar  between  BXD  parental 
strains,  consistent  with  the  similar  kinetics  of  parasite  clearance  and  the  known  dependence  of 
this  on  the  adaptive  T  helper  immune  response  (Figure  2B).  Interestingly,  we  also  noted  that  the 
number  of  CD4/CD25/Foxp3  positive  T  regulatory  cells  was  also  significantly  elevated  in  DBA/2 
versus  C57BL/6  (Figure  2C).  The  significance  of  this  finding  has  yet  to  be  determined. 
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Figure  2.  Immune  response  in  L.  major  infected  C57BL/6  and  DBA/2  mice.  Shown  are  FACS  plots  of  draining  popliteal  lymph 
node  cells.  Each  dot  represents  the  number  of  (A)  MHC  class  II+,  (B)  CD4+  and  (C)  CD4+CD25+Foxp3+  (Treg)  cells  per  LN 
from  individual  C57BL/6  (open  circles)  and  DBA/2  (closed  squares)  mice.  (D)  Shown  is  cytokine  and  chemokine  content  (as 
measured  by  Milliplex  MAP  kit)  in  the  supernatant  of  draining  popliteal  lymph  node  cells  cultured  for  24  hours  with  soluble  L. 
major  antigen  (SLA).  Open  bar:  C57BL/6  mice,  Closed  bar:  DBA/2.  Asterisks  indicate  significant  differences  (P  <  .05). 
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Finally,  to  determine  whether  popliteal  lymph  node  cells  from  infected  DBA/2  and 
C57BL/6  mice  exhibited  functional  differences  that  could  potentially  explain  the  differential 
inflammatory  response,  we  cultured  popliteal  lymph  node  cells  from  mice  infected  for  8  weeks 
with  soluble  leishmania  antigen  (SLA)  for  24  hours  and  then  used  a  Milliplex  bead  assay  to 
measure  cytokine  and  chemokine  content  in  culture  supernatants.  As  shown  in  Figure  2D,  with 
just  a  few  exceptions,  the  soluble  mediator  response  was  significantly  more  vigorous  from 
DBA/2  than  from  C57BL/6.  Taken  together,  our  results  indicate  that  the  BXD  parental  strains 
exhibit  a  very  strong  genetically  programmed  difference  in  their  inflammatory  response  to  Lm. 
Further,  our  results  demonstrate  that  phenotypic  analysis  of  the  BXD  resource  will  be  a  fruitful 
endeavor.  Finally,  our  results  establish  the  optimal  conditions  for  reproducing  a  stereotypic 
infection  course,  optimal  time  points  for  phenotypic  measurement  and  validation  of  an 
informative  set  of  phenotypic  measures  with  which  to  characterize  the  genetic  contribution  to  the 
host  response  to  Lm  infection. 

3.2  Key  Research  Accomplishments 

1.  Developed  an  optimal  dosing  regime  for  reproducible  Lm  infection-dependent  clinical 
course. 

2.  Selected  and  optimized  a  set  of  in  vivo  and  ex  vivo  phenotypic  readouts  with  which  to 
study  the  immunological,  parasitological  and  clinical  aspects  of  Lm  infection  in  mice. 

3.  Identified  optimal  time  points  for  phenotypic  assessments. 

4.  Together,  these  accomplishments  set  the  stage  in  year  2  for  successful  completion  of 
the  phenotypic  analysis  of  the  BXD  resource. 

3.3  Reportable  Outcomes 

Thus  far,  our  studies  have  revealed  that  parental  strains  C57BL/6  and  DBA/2  differ  in  their 
immunological  and  clinical  responses  to  Lm.  Interestingly,  these  differences  are  not  correlated 
with  parasitological  score,  suggesting  underlying  genetic  differences  in  inflammatory  responses 
to  Lm.  Thus,  genetic  dissection  of  these  differences  using  the  BXD  resource  is  likely  to  be  highly 
informative  with  regard  to  fundamental  aspects  of  the  inflammatory  response. 

3.4  Conclusions 

Our  initial  studies  have  revealed  that  the  host  response  to  Lm  is  well  suited  for  genetic  analysis 
using  the  BXD  mouse  resource.  Parental  strains  C57BL/6  and  DBA/2  exhibit  highly  significant 
and  reproducible  differences  in  immunological  and  clinical  scores.  Interestingly,  these  are  not 
correlated  with  parasitological  score,  suggesting  that  underlying  genetic  differences  in  the 
inflammatory  response  to  Lm  are  responsible  for  differential  clinical  outcomes.  We  have 
established  an  optimal  dosing  regime  for  obtaining  reproducible  Lm  infection  courses,  suitable 
for  analysis  of  BXD  mice  in  independent  cohorts  over  a  period  of  time.  Further,  we  have 
selected  and  optimized  an  informative  set  of  in  vivo  and  ex  vivo  phenotypic  readouts  to  study 
the  immunological,  parasitological  and  clinical  aspects  of  Lm  infection  in  mice.  Finally,  we  have 
selected  optimal  time  points  for  phenotypic  assessments.  With  these  preliminary  studies  now 
completed,  we  are  well  positioned  to  begin  phenotypic  scoring  of  the  BXD  mice.  At  this  point, 
we  request  additional  funding  to  support  expression  profiling  of  RNA  collected  from  draining 
popliteal  lymph  nodes  of  uninfected  and  infected  mice.  Additional  funding  would  cover  the  cost 
of  preparing  RNA  from  frozen  tissue  saved  in  our  archive,  preparing  reverse  translated  DNA 
libraries  suitable  for  lllumina  deep  sequencing  and  lllumina  sequencing  itself.  The  additional 
data  we  seek  will  permit  detection  of  eQTLs  that  regulate  the  host  response  to  Lm.  Further  by 
the  application  of  Bayesian  analysis,  it  may  be  possible  to  discern  cause/effect  relationships 
between  specific  DNA  variants,  expression  and  phenotype  QTLs. 
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4.  Severe  Acute  Respiratory  Syndrome  (SARS  Co-V)  project 

We  have  fully  executed  the  Material  Transfer  Agreements  with  CDC  and  NIH  for  transferring  the 
wild-type  (Urbani  strain)  and  mouse  adapted  (MA-15)  SARS  coronaviruses  to  UTHSC  and  have 
obtained  both  viruses.  The  personnel  engaged  in  the  SARS  project  have  gone  through  DOJ 
clearance  and  have  been  approved  for  access  to  select  agents  and/or  toxins.  They  have  also 
been  enrolled  in  the  UTHSC  Animal  Care  and  Use  Occupational  Health  &  Safety  Program  and 
fit-tested. 

In  a  related  project,  we  have  further  characterized  the  interferon  (IFN)  antagonism  mediated  by 
the  papain-like  proteases  (PLPs)  of  two  important  human  coronaviruses,  SARS-CoV  and  NL63, 
in  collaboration  with  S.  Baker  at  Loyola  University  -  Chicago.  The  results  show  that  the  PLP2  of 
NL63,  like  what  we  have  previously  demonstrated  for  the  PLP  of  SARS-CoV,  blocks  IFN 
induction  by  inhibiting  the  activation  of  IFN  regulatory  factor-3,  a  latent  transcription  factor 
pivotal  for  type  I  IFN  synthesis.  Furthermore,  the  PLP  IFN  antagonism  is  enhanced  by,  but  is  not 
strictly  dependent  on,  the  catalytic  activity  of  the  PLP  enzyme.  A  manuscript  has  been  published 
to  report  these  findings  (see  below). 

Reportable  outcomes 

Publications 

1.  Clementz  MA,  Chen  Z,  Banach  BS,  Wang  Y,  Sun  L,  Ratia  K,  Baez-Santos  YM,  Wang  J, 
Takayama  J,  Ghosh  AK,  Li  K,  Mesecar  AD,  Baker  SC.  Deubiquitinating  and  interferon 
antagonism  activities  of  coronavirus  papain-like  proteases.  Journal  of  Virology  2010;  84: 
4619-4629. 
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5.  Systems  Genetics  using  a  prototype  model  of  Chlamydia  psittaci  genetic 
susceptibility.  - 


5.1  Overview 

In  a  previous  study,  we  determined  that  polymorphisms  in  two  genes  (Iigp2  and 
IrgblO)  between  B6  and  D2  mice  account  for  differential  susceptibility 
(Difference  in  100%  Lethal  dose  >100,000  fold)  to  infection  with  C.  psittaci. 

This  model  was  investigated  further  to  serve  as  a  proof  of  concept  system  to 
perform  systems  genetic  studies  using  the  BXD  mice  in  infectious  diseases. 

We  generated  immune  and  disease  phenotype  data  for  59  BXD  strains  infected 
with  Chlamydia  and  microarray  based  transcriptional  analysis  in  40  BXD 
strains.  These  data  were  used  to  perform  QTL  analysis  for  each  parameter  and 
the  correlation  of  each  parameter  was  also  investigated. 

5.2  Multi-scale  QTL  analysis 

BXD  mice  were  infected  and  were  assessed  for  pathogen  load,  immune 
parameters  (32  cytokine  and  immune  cell  population),  and  disease  outcome. 
Results  indicate  that  while  Chrll  locus  influences  variations  in  many 
parameters  several  other  unique  genetic  loci  control  individual  parameters 
(Figure  1).  In  particular,  we  found  that  some  Th1/Th2  cytokines  (IFNy,  IL10) 
were  under  the  control  of  Chromosome  1  and  8.  The  chromosome  8  locus 
reached  a  statistically  significant  level  with  an  LRS  >20  and  mapped  to 
candidate  genes  Nod2  and  Cyld  that  are  known  to  be  negative  regulators  of 
NFkB  (Figure  2). 

Figure  1  (LEFT).  QTL  analysis  of  disease  phenotypes  (weight  loss  day  6)  and  immune 
phenotypes  (pathogen  load,  macrophage  activation  status,  neutrophil  recruitment,  20 
cytokine  levels:).  X-axis  represents  various  traits  and  Y-axis  represents  position  of  the 
genetic  locus  on  the  mouse  genome  (chromosome  1-19,  X).  Statistically  significant 
association  of  a  genetic  locus  with  a  given  trait  is  given  a  color.  Cold  colors  indicate  that  the 
B6  genotype  increases  trait  value  and  warm  colors  indicate  that  D2  genotype  increases 
value.  Results  indicate  that  the  previously  defined  QTL  on  Chr  11  dictates  immune 
responses  observed.  Multiple  other  genetic  loci  (Chr  1, 8,  9) . 


Figure  2.  QTL  analysis  of  IFNg  level  induced  by  Chlamydia.  The  X  axis  represents  the 
position  on  the  Chromosome  8  and  blue  line  indicates  level  of  association  of  IFNy  level 
with  the  genetic  marker  on  the  chromosome.  The  QTL  at  position  91MB  includes  20 
genes  including  Nod2  and  Cyld. 
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Figure  1 . 


5.3  Modeling  Immune  phenotypes 

In  order  to  generate  a  model  of  how  individual  immune  phenotypes  influence  each  other,  a 
correlation  analysis  was  performed  for  all  immune  phenotypes  using  Pearson  and  Spearman 
correlations.  Strongest  positive  correlation  with  the  disease  phenotype  (weight)  was  found  with 
macrophage  activation  status,  rather  than  Chlamydia  load  (Figure  3).  The  macrophage 
activation  status  itself  had  a  strong  negative  correlation  with  Chlamydia  load  as  well  as 
percentage  of  neutrophils  in  the  inflammatory  cell  population.  The  level  of  cytokines  such  as 
IFNy  and  TNFa  were  not  associated  with  the  weight  changes  observed  in  this  model. 


In  order  to  determine  whether  the 
phenotype  was  a  direct 
consequence  of  the  pathogen 
burden  we  examined  the 
correlation  of  pathogen  load  with 
weights  of  individual  mice 
according  to  the  genotype  at 
chromosome  11  encoding  the 
IRGs.  Results  demonstrated  that 
mice  that  have  the  B6  genotype  at 
the  locus  are  tolerant  to  increases 
in  the  pathogen  load  (Data  not 
shown).  Also,  there  was  a 
significant  overlap  in  the  pathogen 
load  between  two  genotypes. 
These  findings  suggest  the  IRGs 
encoded  in  chromosome  11  locus 
affect  disease  outcome  by  a 
mechanism  of  tolerance  (ability  to 
withstand  increasing  pathogen 
burden)  rather  than  resistance 
(ability  to  restrict  pathogen  load). 


Figure  3:  Graphical  representation  of  the  association  among 
phenotypes  .A  line  between  the  phenotypes  indicate  significant 
correlation.  Warm  colored  lines  represents  a  positive  correlation. 
Cold  colors  represents  a  negative  correlation.  Thicker  lines  (R  value 
>0.7).  Thin  lines  (R  value  0.5-0. 7).  The  strongest  positive  correlation 
with  the  disease  phenotype  (weight)  was  found  with  macrophage 
activation  status,  rather  than  Chlamydia  load.  The  macrophage 
activation  status  itself  had  a  strong  negative  impact  on  Chlamydia 
load  as  well  as  percentage  of  neutrophils  in  the  inflammatory  cell 
population.  The  level  of  cytokines  such  as  IFNy  and  TNFa  were  not 
associated  with  the  weight  changes  observed  in  this  model 
suggesting  cellular  phenotypes  drive  pathogenesis. 
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5.4  Transcriptional  analysis  and  eQTL  mapping. 

We  have  performed  a  microarray  based  transcriptional  analysis  of  peritoneal  lavage  samples 
obtained  from  40  BXD  strains  infected  intraperitoneally  with  Chlamydia  and  an  eQTL  analysis 
using  this  data  set.  Consistent  with  other  studies,  eQTLs  with  the  highest  associations  were 
primarily  cis-acting  QTLs,  however  the  analysis  also  reproduced  trans-acting  QTLs  that  have 
been  reported  by  others  in  the  past  (Figure  4). 


Marker  (MB) 

Figure  4:  eQTL  analysis  of  microarray  transcriptional  profiling  from 
Chlamydia  infected  BXD  strains.  Position  of  the  transcript  on  the 
Chromosome  is  indicated  on  the  Y  axis  and  position  of  the  QTL  that 
controls  transcriptional  I  levels  is  indicated  on  X  axis.  Each  dot 
represents  transcript/QTL  with  LRS>24  with  false  discovery  rate  of 
5%.  Numerous  cis-QTLs  as  well  as  trans-QTLs  are  represented. 


5.5  Reportable  outcomes 

Publications 

1.  Miyairi  I  I,  Laxton  J,  Ziebarth  J,  Williams  R,  Lu  L,  Byrne  G,  Cui  Y.  Systems  Genetics 
Analysis  of  Chlamydia  psittaci  Infection.  2010.  Proc.  12th  Int.  Symposium  on  Human  Chlamydial 
Infections,  Schachter  et  al  eds,  pp  197-200,  June  20-25,  Hof  bei  Salzburg,  Austria. 

(see  Appendices) 
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5.5  Model  construction 

Data  of  high  order  traits  and  transcriptional  analysis  was  combined  to  construct  hierarchical 
model  of  how  genotype  at  a  particular  QTL  influences  gene  transcription  and  outcome  of  high 
order  traits.  Over  10  models  were  constructed.  An  example  is  shown  in  Figure  5.  This  model 
suggests  that  Chr8  locus  (described  in  Figure  2)  controls  IFNy  levels  through  Prdx2.  This  is 
biologically  plausible  association  given  the  presence  of  two  candidate  genes  on  Chr  8  (Nod2 
and  Cyld),  and  because  Prdx2  is  known  to  be  a  negative  regulator  of  NFkB  and  control  Thl 
cytokine  levels.  This  will  require  further  validation. 


Figure  5.  Model  of  how  a  genetic  locus  controls 
downstream  gene  and  phenotype.  This  model 
indicates  that  a  gene/or  genes  in  the  Chr  8 
interval  (described  in  Figure  2)  controls  IFN 
gamma  levels  through  the  regulation  of  Farsa 
or  Prdx2. 
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6. 


Mouse  genomics  core 


6.1  Gene  expression  data  collection  from  BXD  spleen 

We  plan  to  collect  spleen  gene  expression  data  from  50-70  BXD  strains  and  20-30  inbred 
mouse  strains  (2  chips/strain)  for  genetic  analysis  of  DoD  priority  pathogens.  So  far  we  have  run 
the  Affy  mouse  GENE  1 .0  array  for  about  110  samples  including  52  BXD  strains  and  27  inbred 
strains  (see  table  below).  The  RNA  has  been  extracted  for  the  remaining  samples.  Most  of  them 
have  passed  QC  and  are  in  the  queue  to  be  run  on  the  chip. 


BXD  strain 

chip/strain 

BXD  strain 

chip/strain 

Inbred  strain 

chip/strain 

BXD2 

2 

BXD55 

1 

129P3/J 

2 

BXD5 

1 

BXD60 

2 

129S1/SvlmJ 

2 

BXD6 

1 

BXD61 

1 

129X1/SvJ 

2 

BXD8 

1 

BXD62 

1 

B6D2F1 

1 

BXD9 

3 

BXD64 

1 

BALB/cByJ 

2 

BXD  11 

1 

BXD66 

1 

BALB/cJ 

2 

BXD  12 

1 

BXD69 

2 

BTBR  T+  tf/J 

1 

BXD  13 

2 

BXD70 

1 

C57BL/6ByJ 

2 

BXD14 

2 

BXD71 

1 

C57BL/6J 

2 

BXD  16 

1 

BXD73 

2 

C57BL/10J 

1 

BXD  18 

2 

BXD75 

2 

C57BLKS/J 

2 

BXD21 

2 

BXD79 

1 

CBA/CaJ 

1 

BXD22 

1 

BXD80 

2 

DBA/2J 

1 

BXD24 

1 

BXD83 

1 

FVB/NJ 

1 

BXD25 

1 

BXD84 

1 

ILS 

2 

BXD27 

2 

BXD86 

2 

ISS 

1 

BXD32 

1 

BXD87 

1 

KK/HIJ 

1 

BXD33 

1 

BXD89 

3 

LG/J 

2 

BXD34 

2 

BXD90 

2 

LP/J 

2 

BXD39 

2 

BXD92 

1 

NOD/LtJ 

1 

BXD40 

2 

BXD95 

1 

NZB/BiNJ 

1 

BXD42 

1 

BXD96 

1 

NZO/HILtJ 

1 

BXD43 

1 

BXD97 

1 

NZW/LacJ 

2 

BXD44 

2 

BXD98 

1 

PL/J 

1 

BXD48 

2 

BXD99 

2 

PWD/PhJ 

1 

BXD50 

3 

BXD101 

1 

PWK/PHJ 

SJL/J 

WSB/EiJ 

1 

1 

1 

6.2  BXD  mice  colony  for  DoD  priority  pathogens 

We  maintain  more  than  500  cages  for  DoD  including  most  of  Jax  BXD  strains  and  all  of  UT  BXD 
strains,  and  BXD  parental  strains.  During  2009-2010  year,  we  have  provided  about  1000  mice 
for  DoD  projects.  Moreover,  we  are  developing  50  new  BXD  strains  in  order  to  increase  the 
power  of  genetic  analysis.  About  25  new  strains  have  passed  generation  4,  and  20  strains  are  in 
generation  6-8. 
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6.3  Sequence  of  DBA/2J 

The  parental  strains  of  BXD  Rl  mice  are  C57BL/6J  (B6)  and  DBA/2J  (D2).  The  B6  sequence  was 
completed  by  the  NIH  several  years  ago,  but  the  D2  genome-wide  sequence  data  is  still  missing. 
Though  the  BXD  panel  is  widely  used  to  find  the  gene  loci  that  control  many  diseases  and 
phenotypes,  cloning  the  genes  for  these  QTLs  is  quite  challenging  work,  creating  a  bottleneck  in 
the  research  process.  The  sequence  polymorphisms  between  B6  and  D2  are  key  information  for 
QTL  cloning,  especially  eQTL  cloning.  As  such,  we  sequenced  the  D2  genome  by  co-operating 
with  our  colleagues. 

Genomic  DNA  from  littermate  DBA/2J  females  (generation  224)  was  originally  prepared  from  liver 
tissue.  These  fully  inbred  animals  were  born  at  UTHSC,  their  parents  being  DBA/2J  F223 
foundation  stock  obtained  directly  from  the  Jackson  Laboratory  in  mid  2008  (date  of  birth 
07/27/2008). 

We  sequenced  the  genome  of  DBA/2J  using  next-generation  sequencing  technology:  ABI  SOLiD 
in  our  campus  and  lllumina  GA2  in  UCLA.  We  produced  10.04  billion  (6.66  billion  from  ABI  SOLiD 
and  0.38  billion  from  lllumina  GA2)  of  mate-paired  reads  that  align  to  the  B6  reference  mouse 
genome,  of  which  3.41  billion  (143.77  Gb  nucleotides)  have  been  mapped  to  the  genome  by  MAQ 
program  for  reads  from  lllumina  GA2,  and  by  ABI  Corona  Lite  pipeline  for  reads  from  ABI  SOLiD. 
This  indicates  an  average  depth-coverage  of  ~55x  (Table  1).  However,  the  depth-coverage 
increases  to  ~79x  when  we  use  only  MAQ  program  to  map  the  short  reads  from  both  lllumina  GA2 
and  ABI  SOLiD. 

Table  1 :  Summary  of  sequencing  data  for  the  D2  genome 


Sequencing 

method 

Clone 

insert 

size 

#  of  total 

reads 

(Billion) 

#  of  total 

nucleotides 

(Gb) 

#  of 

mappable 

reads 

(Billion)t 

#  of 

mappable 

nucleotides 

(Gb)t 

Sequence 
depth  (X)* 

ABI  SOLiD 

1-2  kb 

1.29 

32.27 

0.80/1.23 

19.95/30.83 

7.67/11.86 

2-3  kb 

4.32 

163.89 

1.92/3.09 

72.39/114.87 

27.84/44.18 

3-4  kb 

1.05 

46.06 

0.34/0.52 

16.28/24.98 

6.26/9.61 

Subtotal 

6.66 

242.22 

3.06/4.84 

108.62/170.23 

41.77/65.47 

lllumina 

GA2 

200  bp 

0.38 

41.87 

0.35 

35.15 

13.52 

Total 

10.04 

284.09 

3.41/5.19 

143.77/205.38 

55.29/78.99 

*Assumed  B6  reference  genome  size  is  2.6  Gb.  Read  depth  was  determined  through  the 
calculation  of  all  mappable  nucleotides  divided  by  the  size  of  B6  reference  genome, 
t  The  number  of  mappable  reads  and  nucleotides  for  data  from  ABI  SOLiD  using  ABI  Corona 
Lite  (top)  and  MAQ  (bottom)  mapping  tools. 

We  detected  approximately  4.6  million  SNPs  variants  between  the  DBA/2J  and  C57BL/6J 
reference  genome,  including  4,235,516  and  3,147,000  SNPs  from  ABI  SOLiD  by  ABI  Corona  Lite 
pipeline  and  lllumina  GA2  by  MAQ  program,  respectively  (Figure  1).  Compared  SNPs  between 
ABI  SOLiD  and  lllumina  GA2,  2.72  million  SNPs  are  shared  between  two  platforms.  Around  0.34 
million  SNPs  were  exclusively  presented  in  GA2  and  1 .3  million  SNPs  were  exclusively 
presented  in  SOLiD  platform  (figure  2).  The  number  of  SNPs  varied  among  chromosomes,  with 
the  highest  rate  on  chromosome  1  and  lowest  rate  on  chromosome  X  (Figure  1).  We  observe 
long  segments  of  either  extremely  high  (40  SNPs  per  10  kb)  or  extremely  low  (0.5  SNPs  per  10 
kb)  polymorphism  rates.  According  to  the  position  of  SNPs,  we  mapped  all  SNPs  to  the  structure 
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of  known  reference  genes.  As  a  result,  we  identified  1,667,877  (36.4%)  SNPs  were  located 
within  gene  regions  (Figure  1).  Of  1,667,877  SNPs,  the  vast  majority  (95.8%)  of  SNPs  was 
located  within  intron  regions,  followed  by  exon  (1.97%),  3’UTR  (1.84%)  and  5’  UTR  (0.33%). 
There  are  average  68  SNPs  per  gene  according  to  the  24,341  reference  genes  from  UCSC. 
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Chromosomes 

Figure  1.  SNPs  distribution  across  each  chromosome  from  SOLiD  and  GA2. 


To  further  assess  our  ability  to  call  genotypes  correctly  by  sequencing,  we  compared  our  data 
with  1 .03  million  SNPs  from  the  public  dbSNP  database.  Approximately,  0.93  million  SNPs  from 
public  dbSNPs  are  shared  with  our  SNPs.  In  addition,  to  confirm  the  new  SNPs  that  we  found, 
we  randomly  selected  about  500  SNPs  that  are  exclusively  presented  in  SOLiD  and  GA2 
platform  and  novel  as  indicated  by  an  absence  in  public  dbSNPs  for  experimental  validation 
using  Sanger  capillary  sequencing  method.  More  than  90%  of  SNP  have  been  validated  in  both 
SOLiD  and  GA2  platform. 
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Figure  2.  Comparison  of  SNPs  between  GA2,  SOLiD  and  dbSNPs. 
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7.  Construction  of  gene  network  models 


7.1  Overview 

During  the  first  year  of  the  grant  support,  we  have  developed  and  integrated  algorithms  (eQTL 
mapping,  Bayesian  network,  SNP  analysis  etc.)  to  model  the  genetic  pathways  underlying  the 
host  responses  to  pathogen  infections.  We  used  BXD  mice  infected  by  Chlamydia  psittaci  to 
test  our  algorithm  because  the  datasets  for  other  pathogens  are  not  yet  available.  The  analysis 
of  data  from  Chlamydia  infected  BXD  mice  will  provide  a  model  for  the  analysis  of  the  analysis 
of  data  from  BXD  mice  infected  by  other  pathogens.  B6  and  D2  mice  have  significant 
differences  in  response  to  infection  with  Chlamydia  [1],  Our  goal  is  to  determine  the  genetic 
variations  between  B6  and  D2  mice  that  result  in  these  differences  to  Chlamydia  susceptibility 
and  form  genetic  networks  that  link  genotype  to  gene  expression  to  phenotype.  These  networks 
provide  testable  hypotheses,  which  can  be  used  to  validate  the  networks  and  confirm  the 
importance  of  the  selected  genomic  and  gene  expression  differences  to  disease  outcomes. 

7.2  Adjustment  of  microarray  batch  effects. 

One  common  problem  in  evaluating  microarray  data  is  that  all  samples  in  a  study  cannot  be 
analyzed  simultaneously  because  of  practical  limits  to  the  number  of  samples  that  can  be 
prepared  or  analyzed  at  one  time.  The  microarray  data  obtained  will  then  contain  non-biological 
experimental  errors,  termed  batch  effects,  which  are  correlated  with  the  time  the  sample  was 
analyzed.  Microarray  data  from  44  BXD  strains  were  observed  to  contain  batch  effects,  as 
clustering  of  the  microarray  data  showed  clusters  were  highly  correlated  with  the  batch.  Initially, 
the  microarray  data  was  adjusted  so  that  the  averages  of  the  different  batches  for  each 
transcript  were  identical;  however,  clustering  analysis  showed  that  significant  batch  effects 
persisted.  ComBat  [2],  an  implementation  of  an  empirical  Bayes  method,  was  then  used  to 
adjust  the  data,  as  it  has  been  developed  to  reduce  batch  effects  for  data  where  the  number  of 
samples  is  much  fewer  than  the  number  of  measurements.  ComBat  has  been  successfully 
used  to  correct  for  batch  effects  in  several  previous  microarray  studies  [3-4].  After  adjustment 
with  ComBat,  clustering  of  the  microarray  data  based  on  batch  was  greatly  reduced. 

7.3  eQTL  mapping. 

The  eQTL  mapping  was  performed  with  QTL  Reaper,  and  1000  permutations  of  the  strain  labels 
were  used  to  adjust  for  multiple  testing  and  provide  a  false  discovery  rate  (FDR).  Fig.  1  shows 
the  eQTLs  that  were  found  for  the  infected  mice  with  a  likelihood-ratio  statistic  (LRS)  >  24, 
corresponding  to  a  FDR  or  4%.  The  eQTLs  fall  into  two  groups:  cis-e QTLs,  which  are  located 
on  a  diagonal  line  through  the  origin  with  a  slope  of  1  and  occur  when  the  marker  and  gene  are 
located  less  than  5  Mb  apart,  and  trans-e QTLs,  which  are  scattered  throughout  the  remainder  of 
the  plot  and  occur  when  the  gene  and  marker  are  not  at  the  same  location.  A  total  of  6262  cis- 
eQTLs  and  1343  trans-e  QTLs  were  found  for  a  FDR  of  4%.  Permutation  of  the  strain  labels 
also  allows  for  an  assessment  of  the  FDR  for  eQTL  hotspots  [5],  which  can  be  seen  in  Fig.  1  by 
vertical  bands  of  trans-e  QTLs,  and  are  locations  on  the  genome  where  a  marker  controls 
several  trans- genes.  Using  a  FDR  of  10%,  3  genomic  locations  were  found  to  be  eQTL 
hotspots:  chromosome  13  near  21  Mb,  chromosome  14  near  45  Mb,  and  chromosome  17  near 
20  Mb.  The  next  step  in  the  analysis  of  the  hotspots  is  determining  candidate  genes  near  the 
hotspot  markers  that  result  in  the  downstream  gene  expression  changes. 

7.4.  Bayesian  network  modeling. 

Small,  causal  networks  have  also  been  created  for  BXD  mice  infected  with  Chlamydia  using 
Bayesian  network  (BN)  modeling.  The  genetic  network  models  link  genotype  to  gene 
expression  to  higher  order  phenotypes,  such  as  pathogen  load  and  survival.  Briefly,  the  BNs 
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are  created  by  grouping  genes  and  higher  order  phenotypes  that  map  to  the  same  genomic 
marker  with  standard  QTL  mapping  procedures.  Then,  the  hypothetical  causal  pathways,  which 
link  the  QTL,  gene,  and  phenotype  are  scored,  and  those  groups  in  which  the  model  where  the 
QTL  is  causal  to  the  gene  and  the  gene  is  causal  to  the  phenotype  are  selected.  Models  1  and 
2  in  Fig.  2  show  simple  models  in  which  a  single  gene  influences  the  phenotype.  Additional 
genes  can  be  included  in  the  models,  as  shown  in  Models  3  and  4  of  Fig.  3,  if  they  improve  the 
accuracy  of  the  phenotype  prediction.  Several  of  the  genes  included  in  the  models  of  Fig.  2 
have  been  associated  with  GO  annotations  [6]  which  are  likely  involved  in  response  to 
Chlamydia  infection,  including  Prdx2,  which  is  involved  in  immune  response,  and  Gps2  and 
A/c/e/7,  which  are  associated  with  GTPase  activity,  which  has  been  previously  linked  to  response 
to  Chlamydia  infection  [1],  Fig.  3  shows  a  further  investigation  of  Model  1.  To  create  the  BNs, 
gene  expression  and  phenotype  data  are  discretized  into  high  and  low  values,  and  Fig. 3  shows 
the  probability  distributions  of  B6  or  D2  genotype  and  high  and  low  gene  expression  and 
phenotype  values  if  the  QTL  and  gene  are  fixed  to  certain  values.  Strains  with  the  B6  genotype 
at  the  rsl  3481 068  marker  on  chromosome  1 1  have  low  expression  of  /Vc/e/7  and  low  weight  5 
days  after  infection  (Day5wt),  while  the  D2  genotype  at  this  marker  results  in  higher  expression 
of  A/c/e/7  and  higher  Day5wt.  Additionally,  Fig.  3  shows  that  fixing  the  value  of  Ndell  to  being 
either  low  or  high  improves  the  prediction  accuracy  of  the  phenotype,  in  comparison  to  fixing  the 
genotype  at  the  rsl  3481 068  marker.  These  genetic  network  models  provide  testable 
hypothesis  that  can  be  used  to  verify  the  networks.  For  example,  if  the  Ndell  gene  is  silenced, 
Model  1  predicts  that  Day5wt  should  decrease  regardless  of  the  phenotype  at  the  rsl  3481068 
marker. 

Figure  1.  eQTL  mapping  of  BXD  mice  infected  with  Chlamydia  psittaci.  The  eQTLs 
shown  have  an  LRS  >  24,  corresponding  to  a  false  discovery  rate  of  4%.  The  vertical 
lines  indicate  the  separation  of  chromosomes,  which  are  indicated  with  the  numbers  at 
the  top  of  the  figure. 
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Figure  2.  Causal  pathways  linking  genotype  SNPs  to  gene  expression  levels  to 
phenotypes  in  BXD  mice  infected  with  Chlamydia  psittaci.  Nodes  outlined  in  red  are 
genotype  marker,  purple  nodes  are  genes,  and  blue  nodes  are  higher  order  phenotypes. 
Day5wt  is  the  weight  at  5  days  after  infection,  IFNg  is  the  amount  of  the  the  cytokine 
interferon-gamma  6  days  after  infection,  and  IFU  is  the  pathogen  load  6  days  after 
infection. 


Figure  3.  Further  investigation  into  Model  1  of  Fig.  2.  The  nodes  shaded  in  yellow  show 
the  probability  distributions  of  the  variables  after  fixing  the  node  shaded  in  grey  to  a 
certain  value. 
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7.5.  Key  Research  Accomplishments: 

•  Investigated  methods  to  deal  with  batch  effects  in  microarray  data  and  used  ComBat  to 
reduce  batch  effects  in  microarray  data  of  BXD  mice  infected  with  Chlamydia 

•  Performed  eQTL  mapping  of  Chlamydia  infected  BXD  strains  and  found  3  QTL  hotspots 

•  Constructed  small  Bayesian  network  models  which  link  genotype  to  gene  expression  to 
higher  order  phenotypes  associated  with  Chlamydia  infection 


29 


7.6  Reportable  Outcomes 

Publications 

Ziebarth,  J.,  B.  Li,  I.  Miyairi,  and  Y.  Cui,  Linking  genotype  to  phenotype  with  Bayesian  network 
modeling  of  Chlamydia  infection. 

Poster  presentation  at  UT-ORNL-KBRIN  Bioinformatics  Summit  2010,  March  20,  2010. 

The  Abstract  (see  appendices  2.1)  will  be  published  in  BMC  Bioinformatics. 

7.7  Conclusion: 

The  analysis  of  the  BXD  strains  infected  with  Chlamydia  has  currently  provided  two  results  that 
merit  future  investigation.  eQTL  mapping  of  microarray  data  showed  that  there  are  three 
potential  QTL  hotspots,  where  a  single  genome  location  effects  many  genes  on  several 
chromosomes.  The  next  step  in  this  analysis  is  to  investigate  the  genes  near  the  hotspot 
regions  to  select  the  genomic  and  gene  expression  differences  at  the  hotspot  that  cause  the 
downstream  effects.  Bayesian  network  modeling  casual  pathways  from  genotype  to  gene 
expression  data  to  phenotype  has  yielded  several  interesting  gene  network  models  that  require 
experimental  validation.  Additionally,  the  methods  used  to  perform  the  eQTL  mapping  and 
Bayesian  network  modeling  can  be  easily  adapted  for  use  with  mice  infected  with  other 
pathogens.  We  have  made  significant  progress  with  respect  to  Objective  3  of  the  Statement  of 
Work. 
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Appendices: 


1.  Administrative  and  logistical  matters 

1.1  List  of  personnel  with  percent  effort  on  the  project 


Personnel 

%  Effort 

Gerald  1.  Byrne,  PI 

10 

James  E.  Bina,  Co-1 

10 

Yan  Cui,  Co-1 

10 

Kui  Li,  Co-1 

10 

Lu  Lu,  Co-1 

10 

Mark  A.  Miller,  Co-1 

10 

Isao  Miyairi,  Co-1 

10 

Fabio  C.  Re,  Co-1 

10 

Robert  Williams,  Co-1 

5 

Olaimatu  Mahdi,  Project  Manager 

25 

Xiaowen  Renee  Bina,  Research  Faculty 

75 

Jonathan  Laxton,  Research  Assistant 

100 

Baoming  Lin,  Postdoc 

100 

Jesse  Dylan  Ziebarth,  Postdoc 

100 

Himangi  R.  Jayakar,  Postdoc 

100 

Jyothi  Pavathareddy,  Research  Technician 

100 

Li  Lui,  Postdoc 

100 

Sarah  Hasty,  Research  Assistant 

25 

Lishi  Wang,  Postdoc 

100 

1.2  Request  for  changes  in  personnel  and/or  personnel  efforts 

The  table  above  gives  the  list  of  all  UTHSC  personnel  supported  by  the  grant  in  year  1 .  Some  of 
the  personnel  listed  were  not  supported  for  the  whole  year  (June  2009-May  2010)  depending  on 
when  they  started/stopped  working  on  the  grant.  Personnel  no  longer  supported  by  the  grant 
are  in  bold. 

We  are  requesting  permission  to  make  the  following  changes  in  personnel: 

1 .  To  hire  two  undergraduate  students  during  the  summer  (July  and  August)  to  continue 

the  work  of  Dr.  Mark  Miller’s  former  postdoc  Himangi  Jayakar. 

2.  To  increase  the  percentage  support  of  the  Project  Manager  from  25%  to  50%,  reflecting 

the  increasing  responsibility  associated  with  managing  the  grant. 
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2. 


Publications 


2.1  Linking  genotype  to  phenotype  with  Bayesian  network  modeling  of  Chlamydia 
infection 

Jesse  D  Ziebarth1*,  Bao  Li1,  Isao  Miyairi1,2,  Yan  Cui1,3 

department  of  Molecular  Sciences  and  2Pediatrics,  3Center  of  Integrative  and  Translational 
Genomics,  University  of  Tennessee  Health  Science  Center,  Memphis  TN,  38163,  USA 
*jziebart@uthsc.edu 

Background 

Understanding  the  causal  pathways  that  link  genotypes  with  gene  expression  and  higher  order 
phenotypes,  such  as  disease  susceptibility,  has  been  a  recent  goal  of  systems  genetic  studies. 
Pairwise  correlations  can  select  genes  and  phenotypes  that  are  associated  with  a  common 
genetic  locus;  however,  the  directions  of  these  links  are  more  difficult  to  determine.  We  propose 
a  method  in  which  Bayesian  networks  are  used  to  establish  causal  pathways  from  genotype  to 
gene  expression  to  phenotype  and  apply  it  to  mice  infected  with  Chlamydia  psittaci. 

Results 

Recombinant  inbred  BXD  mice  strains  have  been  previously  used  to  study  the  genetic 
differences  that  cause  a  much  greater  resistance  to  Chlamydia  psittaci  infection  in  C57BL/6J 
mice  compared  with  DBA/2J  mice  [1],  Gene  expression  levels  and  higher  order  phenotypes, 
such  as  pathogen  load  and  survival,  were  collected  for  44  BXD  strains.  Directed  graphs  linking 
SNPs  [2]  in  the  BXD  genomes  to  this  data  were  created  with  QTL  mapping  and  the  Bayesian 
network  method.  Figure  1  shows  a  selection  of  the  resulting  causal  pathways;  many  of  the 
genes  in  these  pathways  are  knows  to  be  associated  with  immune  response  or  GTPase  activity, 
which  has  been  linked  to  Chlamydia  psittaci  infection  [1], 


Figure  1.  Causal  pathways  linking  genotype  SNPs  to  gene  expression  levels  to  phenotypes  in 
BXD  mice  infected  with  Chlamydia  psittaci 
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SYSTEMS  GENETICS  ANALYSIS  OF  CHLAMYDIA  PSITTACI 
INFECTION 


I.  Miyairi1, 2,  J.  Laxton1,  J.  Ziebarth1,  R.  Williams3,  L.  Lu3,  G.  Byrne1,  Y.  Cui1. 

1  2 

Department  of  Molecular  Sciences  ,  Pediatrics  ,  and  Anatomy  and 
Neurobiology  ,  University  of  Tennessee  Health  Science  Center,  Memphis,  TN, 
USA. 

Introduction 

Individual  differences  in  susceptibility  to  infectious  diseases  have  been 
mapped  to  numerous  disease  defining  genetic  polymorphisms  or  quantitative 
trait  loci  (QTL)  by  forward  genetics  (1-3).  However,  such  genotype-phenotype 
associations  do  not  provide  information  on  the  causal  pathways  of  how  disease 
occurs.  Inbred  mouse  models  have  been  widely  utilized  for  identification  of 
susceptibility  genes  and  systems  genetics.  Recombinant  inbred  (RI)  strains 
contain  unique,  approximately  equal  proportions  of  genetic  contributions  from 
two  progenitor  inbred  strains.  A  panel  of  78  established  BXD  RI  strains 
derived  from  C57BL/6J  (B6)  and  DBA/2J  (D2)  strains  have  been  genotyped 
over  13,377  SNP  markers  and  are  available  for  repeated  testing.  Combined 
with  a  web-based  analytical  tool  (GeneNetwork.org)  this  allows  high  precision 
Quantitative  trait  loci  (QTL)  mapping.  In  a  previous  analysis,  we  determined 
that  a  QTL  on  chromosome  11  containing  two  polymorphic  innate  immune 
genes  (Irgm3  and  IrgblO)  in  the  family  of  immunity-related  GTPases  (IRG) 
were  responsible  for  the  innate  difference  in  susceptibility  to  a  systemic 
infection  to  C.psittaci  between  the  B6  and  D2  mice  (3).  Characterization  of  the 
immunological  differences  between  the  B6  and  D2  mice  revealed  significant 
differences  in  Chlamydia  load,  inflammatory  responses,  and  cytokine  profiles. 
While  the  IRGs  were  shown  to  control  Chlamydia  load  (2),  alternative 
immunomodulatory  functions  of  this  gene  have  been  implicated  (4-5).  In  our 
model  it  is  unclear  whether  the  various  immunological  differences  are  a 
function  of  pathogen  load,  immunomodulatory  functions  or  yet  unrelated  to 
the  IRGs. 

We  hypothesized  that  a  multiscale  analysis  of  BXD  genotype,  pathogen  load, 
and  immune  parameters  can  be  combined  to  define  causative  immune 
pathways  induced  by  Chlamydia  psittaci  infection. 


Methods 

1 .  Infection  and  sample  collection. 

Chlamydia  infection:  8-10  week  old  male  mice  (C57BL/6J,  DBA/2J,  and  57 
BXD  strains)  were  infected  in  groups  of  2-5/strain  for  a  total  of  372  mice.  C. 
psittaci  6BC  strains  were  propagated  in  L  cells,  titrated  and  stored  at  -80C°. 
Intraperitoneal  infection  with  C. psittaci  6BC  (104  IFU)  was  performed  using 
the  same  stock  source  to  minimize  variations  across  experiments. 

Monitoring  and  sample  collection:  Infected  mice  were  monitored  daily  for 
weight  changes.  On  day  6-post  infection,  mice  were  euthanized  to  obtain 
peritoneal  lavage  samples  for  pathogen  load,  flow  cytometry,  and  cytokine 
analysis. 

2.  Assessment  of  immune  phenotypes. 

Chlamydia  load:  Titration  was  performed  by  a  cell  culture  based  IFU  assay. 
Flow  cytometry:  Standard  methods  were  used.  Briefly,  murine  peritoneal 
exudates  were  blocked  with  Fc  block  and  incubated  with  fluorochrome- 
conjugated  Abs.  The  following  Abs  was  used:  (Macrophage)  F4/80-APC, 
(Neutrophil)  Ly6G-PE,  (Class  II)  IA/AE-Alexa488,  (Myeloid)  and  CD  11b- 
PE-Cy7.  Data  was  expressed  as  percent  of  macrophages  or  neutrophils  in  the 
CD  lib  expressing  myeloid  cells.  Class  II  expression  was  used  as  a  marker  for 
macrophage  activation  status  and  data  was  expressed  as  percent  of  F4/80 
expression  cells  that  also  expressed  class  II. 

Cytokine  analysis:  Peritoneal  lavage  supernatants  were  stored  in  -80°C  till 
assessment.  We  used  the  Mouse  32-plex  kit  to  analyze  levels  of  32  cytokines. 

3.  Data  analysis. 

QTL  analysis:  We  used  interval  mapping  method  available  on  the 
GeneNetwork.org  interface  to  identify  genetic  loci  regulating 
infection/immune  phenotypes.  Interval  Mapping  is  a  statistical  test  of 
association  between  trait  values  and  the  genotypes  of  marker  loci  through  the 
genome.  A  significant  association  is  interpreted  as  indicating  the  presence  of  a 
QTL  linked  to  the  marker  that  shows  the  association. 

Correlation  analysis:  The  correlation  analysis  was  performed  using  Pearson 
and  Spearman  correlations  and  a  network  graph  was  drawn  using  a  built  in 
function  in  the  GeneNetwork.  This  function  enables  graphical  association 
among  large  groups  of  phenotypes. 


Results 

1 .  Immune  phenotype  data 

Mice  exhibited  a  range  of  disease  manifestation  ranging  from  30  %  weight 
reduction  to  10  %  weight  increase  in  6  days.  Chlamydia  load  ranged  from  4  to 
6.5  log  IFU/mouse.  Neutrophils  ranged  2-40%  of  the  myeloid  population.  Of 
the  macrophage  population  2-80%  expressed  Class  II.  A  measurable 
difference  was  observed  in  at  22  of  the  32  cytokines  examined. 

2.  QTL  analysis 

The  majority  of  the  variability  observed  in  weight  change,  macrophage 
activation  status,  Chlamydia  load,  and  neutrophil  recruitment  was  attributed  to 
the  previously  reported  chromosome  1 1  locus.  While  many  other  traits  were 
influenced  by  the  chromosome  1 1  locus,  subset  of  chemokines  were  controlled 
by  a  combination  of  minor  QTLs  on  chromosomes  4,  5,  9,  and  16  whereas 
IFNy,  TNFa,  IL10  levels  were  under  the  primary  control  of  a  QTL  on 
chromosome  8  (Table). 

3.  Correlation  analysis 

The  strongest  positive  correlation  with  the  disease  phenotype  (weight)  was 
found  with  macrophage  activation  status,  rather  than  Chlamydia  load  (Table). 
The  macrophage  activation  status  itself  had  a  strong  negative  correlation  with 
Chlamydia  load  as  well  as  percentage  of  neutrophils  in  the  inflammatory  cell 
population.  The  level  of  cytokines  such  as  IFNy  and  TNFa  were  not  associated 
with  the  weight  changes  observed  in  this  model. 

In  order  to  determine  whether  the  phenotype  was  a  direct  consequence  of  the 
pathogen  burden  we  examined  the  correlation  of  pathogen  load  with  weights 
of  individual  mice  according  to  the  genotype  at  chromosome  1 1  encoding  the 
IRGs.  Results  demonstrated  that  mice  that  have  the  B6  genotype  at  the  locus 
are  tolerant  to  increases  in  the  pathogen  load.  Also,  there  was  a  significant 
overlap  in  the  pathogen  load  between  two  genotypes.  These  findings  suggest 
the  IRGs  encoded  in  chromosome  11  locus  affect  disease  outcome  by  a 
mechanism  of  tolerance  (ability  to  withstand  increasing  pathogen  burden) 
rather  than  resistance  (ability  to  restrict  pathogen  load). 

Discussion 

Results  indicate  the  involvement  of  multiple  QTL  in  determining  the  immune 
phenotypes  associated  with  a  systemic  Chlamydia  psittaci  infection.  A 
systems  approach  utilized  in  this  study  allows  us  to  infer  immune  responses 
and  its  hierarchy  that  contribute  to  disease  process.  Specifically,  the  disease 
status  controlled  by  the  chromosome  11,  which  is  most  likely  the  IRGs, 


appeared  to  be  a  function  of  the  macrophage  activation  status  rather  than  direct 
Chlamydia  load.  In  this  model  the  macrophage  activation  status  influences 
Chlamydia  load  as  well  as  the  inflammatory  cell  population.  While,  this  will 
require  further  validation,  this  likely  represents  the  multiple  functions  of  the 
IRG.  Cytokine  levels  that  intuitively  would  be  associated  with  disease  severity 
such  as  IFNy  and  TNFa  were  found  to  fluctuate  independent  of  disease  status 
and  were  controlled  by  other  genetic  loci.  Candidate  genes  on  the  chromosome 
8  interval  are  currently  under  investigation. 


Table:  QTL  mapping  results  and  Pearson  correlation  coefficient. 


Phenotype 

QTL 

(Chromosome) 

Macro 

Neutro 

Weight 

Chlamydia  TNFa 

Mip2  IFNg 

Mac  activation  status 

Chrll 

1 

Neutrophils 

Chrll 

-0.726 

1 

Weight 

Chrll 

0.738 

-0.673 

1 

Chlamydia  load 

Chrll,  5 

-0.736 

0.556 

-0.69 

1 

TNFa 

Chr  8,  11 

Chr  4,  5.  9,  11, 

-0.307 

0.404 

-0.25 

0.374 

1 

Mip2 

16 

-0.531 

0.282 

-0.354 

0.269 

0.149 

1 

IFNg 

Chr  1,  8 

-0.261 

0.356 

-0.246 

0.242 

0.712 

0.106  1 

High  positive  correlations  0.7  to  1,  Moderate  positive  correlation  0.5  to  0.7,  High  negative  correlation 
-0.7  to  -1.0,  Moderate  negative  correlation  -0.5  to  -0.7. 
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Abstract 

Background:  Francisella  tularensis  (FT)  is  a  gram-negative  facultative  intracellular  coccobacillus  and  is  the  causal 
agent  of  a  life-threatening  zoonotic  disease  known  as  tularemia.  Although  FT  preferentially  infects  phagocytic  cells 
of  the  host,  recent  evidence  suggests  that  a  significant  number  of  bacteria  can  be  found  extracellularly  in  the 
plasma  fraction  of  the  blood  during  active  infection.  This  observation  suggests  that  the  interaction  between  FT  and 
host  plasma  components  may  play  an  important  role  in  survival  and  dissemination  of  the  bacterium  during  the 
course  of  infection.  Plasminogen  (PLG)  is  a  protein  zymogen  that  is  found  in  abundance  in  the  blood  of 
mammalian  hosts.  A  number  of  both  gram-positive  and  gram-negative  bacterial  pathogens  have  the  ability  to  bind 
to  PLG,  giving  them  a  survival  advantage  by  increasing  their  ability  to  penetrate  extracellular  matrices  and  cross 
tissue  barriers. 

Results:  We  show  that  PLG  binds  to  the  surface  of  FT  and  that  surface-bound  PLG  can  be  activated  to  plasmin  in 
the  presence  of  tissue  PLG  activator  in  vitro.  In  addition,  using  Far-Western  blotting  assays  coupled  with  proteomic 
analyses  of  FT  outer  membrane  preparations,  we  have  identified  several  putative  PLG-binding  proteins  of  FT. 

Conclusions:  The  ability  of  FT  to  acquire  surface  bound  PLG  that  can  be  activated  on  its  surface  may  be  an 
important  virulence  mechanism  that  results  in  an  increase  in  initial  infectivity,  survival,  and/or  dissemination  of  this 
bacterium  in  vivo. 


Background 

Francisella  tularensis  (FT)  is  a  Gram-negative  intracellu¬ 
lar  pathogen  that  is  the  etiological  agent  of  a  multi-syn¬ 
dromic  disease  with  a  high  morbidity/mortality  that  is 
referred  to  as  tularemia.  The  pneumonic  form  of  tulare¬ 
mia  is  of  particular  concern  because  of  the  high  mortal¬ 
ity  rate  (up  to  60%)  following  inhalation  of  as  few  as  ten 
organisms  [1-4].  Francisella  species  are  found  through¬ 
out  the  Northern  Hemisphere  and  infect  a  variety  of 
vertebrate  and  invertebrate  hosts  [5,6].  Infections  with 
FT  can  be  contracted  from  blood  sucking  insects,  such 
as  the  deer  fly  [5,7],  mosquitoes  [8,9],  and  ticks  [5,7,10], 
and  by  open-wound  contact  with  infected  animal  tissue 
[5,11,12], 

Upon  entry  into  a  susceptible  vertebrate  host,  FT  is 
readily  phagocytized  by  resident  macrophages  and  den¬ 
dritic  cells  and  quickly  escapes  into  the  cytoplasm 
[13,14]  where  it  multiplies.  Late  in  its  replicative  cycle, 
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FT  induces  apoptotic  death  of  the  host  phagocyte, 
resulting  in  release  of  progeny  bacteria  that  can  infect 
new  host  cells.  Recent  studies  have  shown  that  signifi¬ 
cant  numbers  of  FT  are  found  in  the  acellular  plasma 
fraction  of  mice  infected  intradermally  or  intranasally 
with  either  FT  Live  Vaccine  Strain  (LVS)  (Type  B)  or 
FT  Schu  S4  (Type  A)  [15],  and  intranasally  with  FT 
novicida  [16].  These  findings  suggest  that,  in  addition  to 
utilizing  the  intracellular  cytoplasmic  niche  for  replica¬ 
tion  and  protection  from  humoral  immunity,  FT  may 
also  have  a  significant  extracellular  phase.  Several  stu¬ 
dies  have  shown  that  deposition  of  host  complement 
component  C3  on  the  surface  of  FT  is  required  for 
opsonophagocytosis  by  activating  CR3  and  CR4- 
mediated  phagocytosis  by  macrophages  and  dendritic 
cells  [14,17,18].  It  is  also  known  that  FT  is  relatively 
resistant  to  complement-mediated  lysis  [19].  A  recent 
report  suggested  that  resistance  of  FT  to  membrane 
attack  complex-mediated  lysis  may  be  due  (at  least  in 
part)  to  its  ability  to  bind  to  factor  H  from  host  plasma 
[20].  It  is  possible  that  the  ability  of  FT  to  bind  to  factor 
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H  and  potentially  to  other  host  plasma  components 
plays  a  significant  role  in  its  pathogenesis. 

It  has  been  long  established  that  a  broad  spectrum  of 
both  gram-positive  and  gram-negative  bacterial  patho¬ 
gens  gain  a  survival  advantage  by  interacting  with  com¬ 
ponents  of  the  host  coagulation/fibrinolytic  system  in 
humans  [21-24].  For  instance,  the  ability  to  acquire  sur¬ 
face-associated  plasmin  has  been  documented  as  an 
important  virulence  mechanism  in  Group  A  streptococci 
[25],  Borrelia  burgdorferi  [26],  and  Yersinia  pestis  [27] 
by  aiding  in  the  organism’s  ability  to  penetrate  the 
extracellular  matrix  and  to  disseminate  to  distal  sites  in 
the  host.  Plasminogen  (PLG)  is  a  92-kDa  glycoprotein 
zymogen  that  is  involved  in  fibrinolysis.  This  precursor 
protein  is  converted  to  an  active  serine  protease  (plas¬ 
min)  by  cleavage  of  the  peptide  bond  between  residues 
R’60and  V561  in  vivo  via  urokinase-type  (uPA)  and/or 
tissue-type  (tPA)  PLG  activators.  Plasmin  has  an  impor¬ 
tant  role  in  blood  clot  resolution  because  of  its  role  in 
the  degradation  of  fibrin  polymers.  Because  plasmin  has 
other  substrates  that  include  pro-collagenases,  pro- 
metalloproteinases,  and  extracellular  matrix  proteins, 
such  as  fibronectin,  laminin,  and  vitronectin,  the  ability 
of  a  bacterium  to  acquire  surface-associated  plasmin  can 
result  in  an  enhanced  ability  of  the  pathogen  to  pene¬ 
trate  the  extracellular  matrix  and  to  disseminate  to  dis¬ 
tal  sites  in  the  host  [21,23,25].  In  this  report  we  show 
that  PLG  binds  to  the  surface  of  FT  in  vitro  and  that 
surface-bound  PLG  can  be  converted  to  the  active  plas¬ 
min  form.  In  addition,  using  a  combination  of  Far-Wes¬ 
tern  blotting  analyses  coupled  with  proteomic 
methodologies,  we  have  identified  several  FT  proteins 
that  can  bind  to  human  PLG  in  vitro. 

Results 

Binding  of  PLG  from  fresh  human  plasma  to  the  surface 
of  FTLVS 

We  used  an  ELISA  assay  to  determine  that  PLG  in  fresh 
frozen  plasma  (FFP)  binds  to  FTLVS  grown  to  mid-log 
phase  in  BHI  (Figure  1).  Binding  was  inhibited  when  e- 
aminocaproic  acid  (eACA),  known  to  inhibit  binding  of 
PLG  to  lysine  groups  in  proteins,  was  included  in  the 
incubation  mixture.  To  help  eliminate  the  possibility  of 
non-specific  binding  of  PLG  due  to  its  high  concentra¬ 
tions  in  human  plasma  and  also  to  rule  out  the  contri¬ 
butions  of  other  plasma  proteins,  we  used  purified 
human  Glu-PLG  (huPLG)  and  noted  similar  results  to 
those  observed  when  FFP  was  used  (Figure  2A).  We 
also  found  that  huPLG  binds  to  the  highly  virulent  Schu 
S4  strain  of  FT  at  moderately  higher  levels  than 
observed  with  FTLVS  (Figure  2B).  We  confirmed  that 
binding  of  huPLG  to  FT  is  a  lysine-dependent  interac¬ 
tion  by  showing  that  increasing  concentrations  of  eACA 
can  inhibit  binding  of  huPLG  to  FTLVS  in  a  dose- 
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Figure  1  FT  binds  to  PLG  from  human  plasma  FTLVS  cultured  to 
mid-log  phase  in  BHI  broth  were  bound  to  wells  of  microtiter  plates 
and  then  incubated  for  1  hour  with  fresh  frozen  human  plasma 
(FFP)  in  the  presence  or  absence  of  100  mM  e-amino  caproic  acid 
(eACA),  a  PLG-binding  inhibitor.  A  modified  ELISA  was  performed  to 
measure  FTLVS-bound  PLG.  The  results  shown  are  representative  of 
3  experiments  of  similar  design.  Bars  indicate  +/-  SEM  in  triplicate. 
Statistical  analysis  was  performed  via  one-way  ANOVA  using  a 
Dunnett's  Multiple  Comparison  post-test  (***  P  <  .001). 


dependent  fashion  (Figure  3).  When  similar  concentra¬ 
tions  of  glycine  were  used  as  an  inhibitor  control,  no 
inhibition  of  huPLG  binding  was  observed  (data  not 
shown).  Confocal  microscopic  analyses  suggested  that 
huPLG  binds  to  the  surface  of  FT  (Figure  4);  however,  it 
is  possible  that  some  of  the  staining  observed  was  the 
result  of  huPLG  penetration  into  the  outer  envelope  of 
FT.  Although  is  has  been  reported  that  culture  media 
composition  can  have  a  significant  impact  of  the  surface 
properties  and  virulence  characteristics  of  FTLVS  [28], 
we  observed  no  differences  in  the  ability  of  PLG  to  bind 
to  the  surface  of  FTLVS  grown  in  modified  Mueller- 
Hinton  medium  vs.  brain-heart  infusion  broth  (data  not 
shown). 

Plasmin  activation  on  the  surface  of  FT  LVS  in  vitro  by  a 
PLG  activator 

In  other  bacterial  systems,  surface-bound  PLG  can  be 
converted  to  its  proteolytically  active  plasmin  form  that 
contributes  to  the  organism’s  virulence  [21-24].  To  test 
whether  huPLG  bound  to  FTLVS  can  be  converted 
to  plasmin,  we  used  a  chromogenic  plasmin  substrate 
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Figure  2  Purified  huPLG  binds  to  FTLVS  and  FTSchuS4.  FTLVS 
(Panel  A)  and  FTSchuS4  (Panel  B)  were  bound  to  microtiter  wells 
and  incubated  for  2  hours  with  purified  huPLG  (3  pg/ml)  in  the 
presence  or  absence  of  1 0  mM  eACA).  A  modified  ELISA  was 
performed  to  measure  FTLVS-bound  huPLG.  The  results  shown  are 
representative  of  four  (Panel  A)  and  one  (Panel  B)  experiments, 
respectively,  of  similar  design.  Bars  indicate  +/-  SEM  in  triplicate. 
Statistical  analysis  was  performed  via  one-way  ANOVA  using  a 
Dunnett's  Multiple  Comparison  post-test  (***  P  <  .001). 


PLG:  -+  +  +  +  +  + 
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Figure  3  eACA  inhibits  huPLG  binding  to  FT  in  a  dose- 
dependent  fashion.  FTLVS  was  coated  onto  microtiter  plate  wells 
and  incubated  for  2  hours  with  purified  huPLG  (3  pg/mL)  in  the 
presence  or  absence  of  titrated  concentrations  of  eACA.  The  results 
shown  are  representative  of  3  experiments  of  similar  design.  Bars 
indicate  +/-  SEM  in  triplicate.  Statistical  analysis  performed  via  one¬ 
way  ANOVA  using  a  Kruskal-Wallis  test  determined  a  p-value  of 
<  0.0001. 


(H-D-Val-Leu-Lys-pNA)  to  detect  proteolytic  activity 
following  the  addition  of  tissue  PLG  activator  (tPA) 
(Figure  5).  We  also  found  that  plasmin  on  the  surface  of 
FT  can  break  down  fibronectin  (Figure  6),  suggesting 
that  FT-bound  plasmin  can  potentially  participate  in  the 
degradation  of  extracellular  matrices. 

Identification  of  putative  Francisella  PLG-binding 
membrane  proteins 

We  used  Far-Western  Blot  methodology  and  mass  spec¬ 
trometry  to  identify  potential  PLG  receptors  in  the  Sar- 
kosyl  soluble  and  insoluble  FTLVS  membrane  fractions 
(Figure  7).  Sarkosyl  is  a  weak  anionic  detergent  in  which 
many  outer  membrane  proteins  of  Gram-negative  bac¬ 
teria  are  insoluble  [29].  We  transferred  the  Sarkosyl- 
treated  proteins  to  a  PVDF  membrane  and  incubated 
the  membrane  with  PLG  and  identified  bound  PLG  by 
reaction  with  anti-PLG  mAbs  (Figure  7a).  We  used  the 
relative  migration  rates  of  the  reactive  bands  to  identify 
the  reactive  proteins  on  a  duplicate  Coomassie-stained 
polyacrylamide  gel  (Figure  7b),  which  were  then  excised 
for  proteomic  analysis  by  mass  spectrometry.  Several 
prominent  PLG-binding  proteins  were  noted  in  the  total 
membrane  fraction  of  FTLVS,  all  but  one  of  which  was 
found  in  the  Sarkosyl  insoluble  fraction  (Figure  7b).  The 
identity  of  the  prominent  proteins  from  this  assay  (Fig¬ 
ure  7c)  are  the  products  of  the  following  genes: 
FTL  1328  (outer  membrane  associated  protein,  fopAl), 
FTL  1042  (FKBP-type  peptidyl-prolyl  cis-trans  isomer- 
ase  family  protein),  FTL_0336  (peptidoglycan-associated 
lipoprotein),  FTL_0421  (hypothetical  lipoprotein,  lpn-A), 
and  FTL_0645  (hypothetical  lipoprotein). 

Discussion 

Until  recently  FT  has  been  considered  an  intracellular 
pathogen  whose  dissemination  to  tissues  distal  to  the 
site  of  initial  infection  was  highly  dependent  on  its  abil¬ 
ity  survive  within  host  macrophages.  The  observation 
that  FT  can  be  found  in  relatively  high  numbers  in  the 
acellular  plasma  fraction  of  its  mammalian  host  [15,16] 
suggested  that  FT  may  have  a  significant  extracellular 
component  to  its  life  cycle  and  that  interactions  between 
FT  and  one  or  more  plasma  proteins  could  contribute 
to  its  ability  to  disseminate  within  the  host.  There  are  a 
number  of  examples  of  bacterial  pathogens  that  utilize 
interactions  with  host  plasma  components  to  enhance 
their  ability  to  colonize  and  to  penetrate  the  extracellu¬ 
lar  matrices  of  host  cells/tissues.  A  wide  range  of  bacter¬ 
ial  pathogens  (including  Francisella)  subvert  the 
destructive  mechanisms  of  the  complement  cascade  by 
acquiring  surface-bound  complement  control  proteins 
[20,30-34].  Moreover,  a  number  of  Gram-positive  bac¬ 
terial  pathogens  including  streptococcal  spp.  [35,36],  sta¬ 
phylococcal  spp.  [37-40],  and  Bacillus  anthracis  [41,42], 
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Figure  4  PLG  binds  to  the  outer  envelope  of  FT.  Laser  scanning  confocal  microscopy  of  PLG-associated  FTLVS  was  performed  as  described  in 
"Materials  and  Methods".  Bound  huPLG  ligand  was  detected  using  sheep  anti-human  PLG  antibody  followed  by  incubation  with  Dylight-488 
conjugated  donkey,  anti-sheep/goat  IgG  secondary  antibody.  Samples  were  visualized  using  a  Zeiss  LSM  510  confocal  microscope. 


PLG:  +  + 

tPA:  +  + 

Figure  5  FT  surface-bound  huPLG  can  be  converted  to 
plasmin.  FTLVS  was  incubated  with  huPLG  at  a  concentration  of  96 
pg/mL.  After  removal  of  unbound  huPLG,  a  chromogenic  plasmin 
substrate  (D-VLK-pNA),  tissue  PLG  activator  (tPA),  or  both  were  then 
added  to  test  the  proteolytic  ability  of  each  sample  preparation. 
Conversion  of  the  chromogenic  substrate  was  measured  by 
comparison  of  A405  nm.  The  results  shown  are  representative  of  3 
experiments  of  similar  design.  Bars  indicate  +/-  SEM  in  triplicate. 
Statistical  analysis  was  performed  via  one-way  ANOVA  using  a 
Dunnett's  Multiple  Comparison  post-test  (***  P  <  .001). 


?i8=»  *  • 


100— 


50 

37 


25— 

FTA  alone:  P  P  t 

FT/PLG:  t  t  P 

tPA :  t  P  P 

Figure  6  Fibronectin  is  a  substrate  for  plasmin  bound  to  FT 

FTLVS  (109  CFU)  were  incubated  with  100  gg  of  huPLG  and  0.5  pg 
tissue  tPA  for  1  hour  at  37°C.  After  removal  of  unbound  huPLG  and 
tPA,  3  pg  fibronectin  was  added  and  allowed  to  incubate  for  24 
hours  at  37°C.  Supernatant  from  each  preparation  were  separated 
by  SDS-PAGE  and  transferred  to  PVDF  membrane.  Degradation  of 
fibronectin  was  detected  by  Western  blot  analysis  as  described  in 
"Materials  and  Methods”. 
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Figure  7  Identification  of  putative  PLG-binding  proteins  of  FT.  Sarkosyl-soluble  and  insoluble  protein  fractions  of  FTLVS  were  separated  by 
SDS-PAGE  and  transferred  to  PVDF  membrane.  Membranes  were  then  blotted  with  huPLG  (3  ug/mL)  followed  by  anti-PLG  antibody  and  HRP- 
conjugated  secondary  antibody  to  detect  PLG-binding  proteins  (Panel  A).  Protein  bands  on  an  identical  Coomassie  Blue-stained  SDS-PAGE  gel 
corresponding  to  those  identified  via  blotting  (Panel  B)  were  excised  and  identified  using  proteomic  methodologies  (Panel  C). 


as  well  as  Gram-negative  bacteria  such  as  Pseudomonas 
aeruginosa  [43]  have  been  shown  to  augment  their  inva¬ 
sive  capacity  by  interacting  with  fibrinogen,  fibronectin, 
and/or  PLG.  Yersinia  pestis  is  probably  the  best-charac¬ 
terized  example  of  a  pathogen  that  exploits  the  host 
fibrinolytic  system  to  penetrate  host  tissues.  Yersinia 
expresses  a  surface  serine  protease  (designated  Pla) 
whose  substrates  include  several  complement  compo¬ 
nents,  PLG,  and  alpha2-antiplasmin  (the  primary  circu¬ 
lating  inhibitor  of  plasmin).  Pla  also  has  adhesin  activity 
and  binds  to  laminin  (a  glycoprotein  of  mammalian 
basement  membranes).  Because  Pla  upregulates  plasmin 
activity,  and  because  laminin  is  a  substrate  of  plasmin, 
Yersinia  can  very  efficiently  penetrate  basement  mem¬ 
branes  of  host  tissues  [for  review,  see  Suomalainen  et. 
al.  [44]].  Clearly,  interaction  with  plasma  components  is 
a  strategy  that  is  used  by  many  bacterial  pathogens  to 
gain  a  survival  advantage  within  their  hosts. 

The  goal  of  the  studies  described  here  was  to  deter¬ 
mine  whether  FT  has  the  potential  to  use  the  host  fibri¬ 
nolytic  system  (specifically  PLG)  to  enhance  its  ability  to 
penetrate/disseminate  following  infection  of  a  mamma¬ 
lian  host.  Our  results  indicate  that  both  FTLVS  and 
FTSchuS4  are  able  to  acquire  surface  bound  PLG  in 
vitro  and  that  this  zymogen  can  be  converted  by  a  host- 
derived  PLG  activator  into  its  active  serine  protease 
form  (plasmin)  while  bound  to  FTLVS.  The  ability  of 
PLG  to  bind  its  ligands  typically  involves  its  lysine-bind¬ 
ing  kringle  domains.  This  specific  interaction  between 
PLG  and  exposed  lysine  residues  can  be  inhibited  with 
the  lysine-analogue  eACA  and,  to  a  lesser  extent,  with 
free  lysine.  Our  findings  revealed  that  binding  of  PLG  to 
the  surface  of  FTLVS  could  be  inhibited  by  eACA  in  a 
dose-dependent  fashion.  Moreover,  we  showed  that 
plasmin  bound  to  the  surface  of  FT  could  degrade 


fibronectin.  This  finding  supports  our  hypothesis  that 
the  ability  of  FT  to  bind  to  serum  plasmin  may  enhance 
its  ability  to  penetrate  extracellular  matrices,  enhancing 
its  ability  to  disseminate  in  vivo. 

Using  a  ligand-blotting  technique  coupled  with  pro¬ 
teomic  methodologies  we  identified  five  FTLVS  proteins 
that  were  able  to  bind  to  PLG,  each  of  which  are  highly 
conserved  among  the  various  FT  type  A  and  B  strains. 
Three  of  these  proteins  are  lipoproteins  (gene  products 
of  FTL  0336,  FTL_0421,  and  FTL_0645).  Two  of  the 
lipoproteins  are  unique  to  FT,  while  the  third,  peptido- 
glycan-associated  lipoprotein  (PAL),  is  highly  conserved 
among  gram-negative  bacteria.  The  specific  use  of  sur¬ 
face-exposed  lipoproteins  as  receptors  for  host  PLG  is 
not  unusual  and  has  been  well  documented  in  other 
human  bacterial  pathogens,  such  as  some  members  of 
the  genus  Borrelia  and  Treponema.  Several  members  of 
the  genus  Borrelia  use  complement  regulator-acquiring 
surface  proteins  (CRASP)  to  bind  both  PLG  and  com¬ 
plement  factor  H  to  aid  in  the  ability  of  the  organism  to 
both  disseminate  and  to  resist  innate  immunity  [45-50]. 
An  additional  example  of  a  PLG-binding  lipoprotein  is 
OppA  of  Treponema  denticola,  which  has  been  sug¬ 
gested  to  play  a  role  in  periodontal  disease  in  humans 
[51].  With  this  in  mind,  there  lies  the  possibility  that 
lipoproteins  of  Francisella  species  may  have  the  capacity 
to  bind  multiple  host-derived  proteins  in  addition  to 
PLG. 

Here  we  have  shown  that  FT  can  bind  to  PLG  and 
that  surface-bound  PLG  can  be  activated  by  tPA  to  its 
proteolytic  form  (plasmin).  The  binding  of  PLG  on  the 
surface  of  FT  could  play  a  role  in  several  phases  of 
tularemia,  including  the  initial  entry  into  the  host 
through  insect  bites  and/or  broken  skin  where  active 
fibrinolytic  processes  would  provide  an  early 
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opportunity  for  FT  to  acquire  proteolytic  activity  that 
might  augment  the  establishment  or  dissemination  of 
infection.  During  later  phases  of  tularemia  the  acquisi¬ 
tion  of  plasmin  on  the  cell  surface  may  contribute  to  its 
pathogenicity  by  degrading  host  innate  effector  mole¬ 
cules  and  extracellular  matrix  components.  Based  on  the 
new  report  that  FT-bound  plasmin  can  degrade  immu¬ 
noglobulins  [52],  as  well  as  the  established  ability  of  FT 
to  acquire  surface-bound  factor  H  [20],  it  also  appears 
likely  that  FT  uses  plasma  components  to  interfere  with 
host  humoral  immune  mechanisms  throughout  the 
course  of  FT  infection.  Future  studies  to  identify  addi¬ 
tional  plasma  components  that  can  be  surface  acquired 
by  FT  may  uncover  additional  virulence  mechanisms 
used  by  this  pathogen  during  its  extracellular  life  cycle. 

Conclusions 

FT  interacts  with  at  least  two  serum  components  (plas¬ 
min,  and  complement  factor  H),  and  it  seems  likely  that 
FT  also  uses  interactions  with  additional  host  serum 
components  to  gain  a  survival  advantage.  Our  lab  is 
examining  FT  interactions  with  additional  targets, 
including  fibrinogen  and  fibronectin,  both  of  which  are 
substrates  for  plasmin  and  are  host  components  that  are 
known  to  be  exploited  by  numerous  pathogens  for  adhe¬ 
sion  to  and  penetration  of  extracellular  matrix  layers. 
The  interaction  of  FT  with  host  serum  components  may 
play  a  significant  role  in  the  survival  and  dissemination 
of  this  highly  pathogenic  bacterium.  Gaining  a  better 
understanding  of  these  interactions  could  be  a  critical 
step  in  the  development  of  therapeutic  and  prophylactic 
interventions  for  tularemic  disease. 

Methods 

Bacterial  strains  and  culture 

F.  tularensis  Live  Vaccine  Strain  (FTLVS)  was  a  kind  gift 
of  Dr.  Karen  Elkins  (FDA,  Bethesda,  MD).  FT  Schu  S4 
was  obtained  from  the  CDC.  All  bacterial  cultures  were 
grown  overnight  in  Brain-Heart  Infusion  broth  (37  g/L, 
pH  6.8)  from  frozen  stocks  at  37°C  with  shaking  to  mid¬ 
log  phase  (OD600  =  ~0.7)  before  use. 

Reagents 

Human  fresh  frozen  plasma  (FFP)  was  purchased  from 
Lifeblood  Mid-South  Regional  Blood  Center  (Memphis, 
TN).  Purified  human  Glu-PLG  (huPLG),  human  single¬ 
chain  tissue  PLG  activator  (tPA),  and  the  plasmin 
colorimetric  substrate  (H-D-Val-Leu-Lys-pNA)  were 
purchased  from  Molecular  Innovations  (Novi,  MI). 
Bovine  serum  albumin  (fraction  V)  was  purchased  from 
Thermo-Fisher  Scientific  (Pittsburgh,  PA).  Polyclonal 
sheep  anti-human  PLG,  anti-human  fibronectin,  and 
donkey  anti-sheep/goat  IgG:Dylight-488  antibody  pre¬ 
parations  purchased  from  AbD  Serotec  (Raleigh,  NC). 


Monoclonal  anti-goat/sheep  IgG-horseradish  peroxidase 
conjugated  secondary  antibody  (clone  GT-34)  and 
E-aminocaproic  acid  (A7824)  were  purchased  from 
Sigma- Aldrich  (St.  Louis,  MO).  Ninety-six  well  MAXI- 
SORP  ELISA  plates  were  purchased  from  Nunc  (Roche¬ 
ster,  NY). 

PLG  binding  ELISA  assays 

FTLVS  was  cultured  overnight  to  mid-log  phase,  pelleted 
at  6,400  x  g  for  30  minutes,  washed  twice  with  phos¬ 
phate-buffered  saline  (PBS),  and  resuspended  in  PBS 
with  0.1%  Na  azide  to  an  OD600  =  0.1.  The  resulting  bac¬ 
terial  suspension  was  added  to  microtiter  plates  (100  pL/ 
well;  approximately  2.5  x  108  bacterial  cells)  before  being 
incubated  overnight  at  4°C  to  facilitate  binding.  The  wells 
were  then  washed  twice  with  200  pL  of  Tris-buffered  sal¬ 
ine  (TBS)  pH  7.45  containing  0.05%  Tween-20  (TBST)  to 
remove  unbound  bacteria  and  then  pre-blocked  with  200 
pL  of  TBST  containing  1%  bovine  serum  albumin  (1% 
BSA-TBST)  for  1  hour  at  RT°  to  prevent  non-specific 
protein  binding.  After  removal  of  the  blocking  solution, 
90%  citrated  human  plasma  or  3  pg/mL  huPLG  in  1% 
BSA-TBST  was  added  to  each  well  (100  pL),  with  or 
without  the  indicated  concentrations  of  e-amino  caproic 
acid  (eACA),  and  incubated  for  1-2  hours  at  37°C  with 
gentle  rocking.  Wells  were  washed  three  times  with 
TBST  and  then  sheep  anti-human  PLG-specific  antibody 
(1:2,000  dilution  in  1%  BSA-TBST)  was  added  (100  pL/ 
well)  and  allowed  to  incubate  for  1  hour  at  37°C. 
Unbound  primary  antibodies  were  removed  by  washing 
three  times  with  TBST,  followed  by  the  addition  of  HRP- 
conjugated  anti-sheep/goat  IgG  mAb  (GT-34,  1:5,000 
dilution  in  1%  BSA-TBST;  100  pL/well)  and  incubation 
for  1  hour  at  37°C.  Unbound  secondary  antibodies  were 
removed  by  washing  four  times  with  TBST,  and  OptEIA 
TMB  colorimetric  substrate  solution  (Becton-Dickenson, 
Franklin  Lakes,  NJ)  was  added  to  each  well  (100  pL/well) 
and  incubated  at  37°C  for  20  min.  to  allow  color  develop¬ 
ment.  Absorbance  at  450  nm  was  determined  using  a 
SpectraMAX  340  plate  reader  (Molecular  Devices,  Sun¬ 
nyvale,  CA). 

Indirect  immunofluorescence  assays 

FTLVS  was  cultured  and  washed  as  described  above. 
After  diluting  the  washed  bacteria  to  OD60o  =  0.1,  1  mL 
aliquots  were  incubated  with  a  total  of  40  pgs  of  PLG  or 
PBS  (negative  control)  for  30  minutes  at  37°C  with  gen¬ 
tle  rotation.  Bacteria  were  then  washed  three  times  with 
PBS  by  centrifugation,  resuspended  in  100  pL  of  PBS, 
followed  by  spotting  20  pL  of  each  sample  onto  glass 
coverslips.  The  samples  were  then  air-dried  overnight  at 
37°C.  After  methanol  fixation,  the  coverslips  were 
blocked  with  1%  BSA-PBS  at  room  temperature  before 
adding  sheep  anti-human  PLG  (1:100  diluted  in  1% 
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BSA-PBS)  for  30  minutes  at  room  temperature.  The 
coverslips  were  gently  washed  with  PBS  before  adding 
donkey  anti-sheep/goat  IgG:Dylight-488  (1:100  diluted 
in  1%  BSA-PBS),  followed  by  incubation  for  30  minutes 
at  room  temperature.  After  washing  again  with  PBS, 
coverslips  were  mounted  onto  glass  slides  using  100% 
glycerol  containing  0.1  M  n-propyl  gallate  and  images 
were  collected  on  a  Zeiss  LSM  510  confocal  microscope 
with  an  Axiovert  100  M  base  with  a  lOOx  Plan  Apoc- 
hromat  1.4  NA  oil  DIC  objective  using  the  argon  laser 
for  488  nm  excitation  and  505-530  nm  bandpass  emis¬ 
sion  filter  for  imaging  Dylight488  fluorescence  and  the 
HeNel  543  nm  laser  for  illumination  of  the  DIC  images. 
Both  images  were  collected  using  identical  detector  gain 
and  amplifier  offset  settings,  and  the  images  shown  are 
1.0  pm  optical  slices.  Digital  images  were  visualized 
using  Zeiss  AxioVision  LE  software. 

Chromogenic  plasmin  activation  assay 

FTLVS  was  cultured  overnight  to  mid-log  phase,  washed 
twice  with  TBS  and  then  resuspended  in  TBS  to  an 
ODgoo  of  0.7.  Aliquots  of  the  bacterial  suspension 
(50  pL)  was  added  to  50  pL  of  TBS  alone  or  TBS  con¬ 
taining  huPLG  (192  pg/mL)  and  incubated  for  1  hour  at 
37°C.  The  cells  were  washed  3x  with  TBST  containing 
0.1%  BSA,  and  pellets  were  resuspended  in  200  pL  of 
TBS  and  then  split  into  two  100  pL  aliquots.  50  pL  of 
50  mM  Tris-HCl  (pH  7.45)  with  or  without  333  pM  of 
the  chromogenic  plasmin  substrate  (H-D-Val-Leu-Lys- 
pNA)  and  50  pL  1.2  pg  of  tPA  or  TBS  alone  was  added 
to  each  sample  and  incubated  at  37°C  for  3  h.  Bacteria 
were  pelleted  via  centrifugation  and  150  pL  of  each 
supernatant  was  pipetted  into  a  96-well  plate  and  absor¬ 
bance  at  405  nm  was  determined  as  a  measure  of  plas¬ 
min  activity. 

Membrane  protein  fractionation 

Outer  membrane  enriched  fractions  were  isolated  by  a 
procedure  adapted  from  de  Bruin,  et  al  [53].  FTLVS 
were  grown  in  BHI  broth  (500  ml)  to  mid-log  phase  and 
then  were  pelleted  via  centrifugation  at  6,400  x  g  for  30 
minutes.  Cells  were  resuspended  in  cold  PBS  and  then 
lysed  by  sonication.  Unlysed  bacterial  cells  were  sepa¬ 
rated  from  the  whole-cell  lysate  by  centrifugation  at 
10,000  x  g  for  20  minutes  at  4°C.  The  insoluble  mem¬ 
brane  fraction  was  then  isolated  by  ultracentrifugation 
for  1  hour  at  100,000  x  g  at  4°C.  After  removal  of  the 
soluble  protein  fraction,  the  pelleted  total  membrane 
fraction  was  resuspended  in  1%  sarkosyl  with  vortexing 
and  subjected  to  a  second  round  of  ultracentrifugation 
for  1  hour  at  100,000  x  g  at  4°C.  The  Sarkosyl-insoluble 
pellet  was  resuspended  in  50  mM  Tris  pH  8.  The  pro¬ 
tein  concentration  of  both  the  Sarkosyl-soluble  and  Sar¬ 
kosyl-insoluble  fractions  was  determined  using  the  DC 


protein  assay  (Bio-Rad,  Hercules,  CA)  according  to 
manufacturer  directions.  Samples  were  stored  at  -20°C 
until  use. 

Fibronectin  degradation  assay 

Overnight  cultures  of  FTLVS  were  washed  three  times 
with  PBS,  109  CFU  were  pipetted  into  1.5  mL  tubes,  and 
bacteria  were  pelleted  via  centrifugation  at  18,900  x  g 
for  10  minutes.  Bacterial  pellets  were  then  resuspended 
in  50  pi  of  PBS  with  or  without  PLG  (2  mg/ml),  fol¬ 
lowed  by  the  addition  of  50  pi  of  tPA  (10  pg/mL)  and 
incubation  at  37°C  with  gentle  shaking  for  1  hour.  The 
bacterial  suspensions  were  pelleted  via  centrifugation  at 
18,900  x  g,  washed  3x  with  PBS  and  resuspended  with 
100  pi  of  50  mM  Tris,  100  mM  NaCl,  5  mM  CaCl2) 
with  3  pg  fibronectin  (BD  Biosciences)  and  incubation 
at  37°C  with  gentle  shaking  for  24  hours.  After  the  incu¬ 
bation  was  complete,  bacteria  were  pelleted  via  centrifu¬ 
gation  at  18,900  x  g  and  the  supernatants  were 
solublized  by  boiling  in  2x  SDS-PAGE  sample  buffer 
containing  2-mercaptoethanol.  Samples  were  subjected 
to  10%  SDS-PAGE  and  then  electrophoretically  trans¬ 
ferred  to  a  PVDF  membrane  (Immobilon-P,  Millipore). 
The  PVDF  membrane  was  pre-blocked  with  1%  BSA- 
TBST  for  1  hour  at  RT  to  minimize  non-specific  protein 
binding,  and  was  then  incubated  with  sheep  anti-human 
fibronectin-specific  antibody  (diluted  1:2000  in  1%  BSA- 
TBST)  for  1  hour  at  RT  with  gentle  rocking.  The  PVDF 
membrane  was  washed  three  times  with  TBST  to 
remove  unbound  primary  antibody.  The  membrane  was 
then  incubated  in  a  solution  of  anti-sheep/goat  IgG 
monoclonal  antibody  (GT-34,  diluted  1:5000  in  1%BSA- 
TBST)  with  rocking  for  1  hr  at  RT.  The  PVDF  mem¬ 
branes  were  washed  3  times  with  TBST  to  remove 
unbound  secondary  antibody.  The  blot  was  developed 
using  Pierce  PicoWest  chemiluminescence  reagents  and 
images  were  captured  using  a  Bio-Rad  ChemiDoc  XRS 
system. 

Far-Western  blotting  analysis 

Approximately  100  pg  of  each  protein  fraction  was  pre¬ 
cipitated  using  ice-cold  acetone,  pelleted  via  centrifuga¬ 
tion  at  18,900  x  g  for  15  minutes,  and  air-dried  at  room 
temperature.  The  samples  were  then  solublized  by  boil¬ 
ing  in  lx  SDS-PAGE  sample  buffer  containing  2-mer- 
captoethanol.  Duplicate  20  pL  aliquots  of  each  sample 
were  subjected  to  15%  SDS-PAGE  to  separate  the  pro¬ 
teins  based  on  their  size.  One  set  of  the  samples  was 
then  electrophoretically  transferred  to  a  PVDF  mem¬ 
brane  (Immobilon-Psq,  Millipore).  The  PVDF  mem¬ 
brane  was  pre-blocked  with  1%  BSA-TBST  for  1  hour  at 
room  temperature  to  minimize  non-specific  protein 
binding  and  was  then  incubated  in  a  solution  of  huPLG 
(3  ug/mL  in  1%  BSA-TBST)  for  one  hour  with  rocking 
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at  37°C.  Unbound  PLG  was  removed  by  washing  three 
times  with  TBST.  Sheep  anti-human  PLG-specific  anti¬ 
body  (diluted  1:2,000  in  1%  BSA-TBST)  was  added 
(100  pL/well)  and  allowed  to  incubate  for  1  hour  at  RT° 
with  rocking.  The  PVDF  membrane  was  washed  three 
times  with  TBST  to  remove  unbound  primary  antibody. 
The  membrane  was  then  incubated  in  a  solution  of 
anti-sheep/goat  IgG  monoclonal  antibody  (GT-34, 
diluted  1:5,000  in  1%BSA-TBST)  with  rocking  for  1  hr 
at  room  temperature.  The  PVDF  membranes  were 
washed  three  times  with  TBST  to  remove  unbound  sec¬ 
ondary  antibody.  The  blot  was  developed  using  Pierce 
PicoWest  chemiluminescence  reagents  and  imaged 
using  a  Bio-Rad  ChemiDoc  XRS  system. 

Proteomic  identification  of  PLG-binding  FT  proteins 

Protein  bands  were  excised  from  Coomassie-stained 
SDS-PAGE  gels,  cut  into  small  pieces,  incubated  in  50% 
acetonitrile/ 100  mM  ammonium  bicarbonate  until  col¬ 
orless,  and  dried  via  vacuum  centrifugation.  The  protein 
was  digested  by  adding  20  pi  of  a  20  ng/pl  trypsin  solu¬ 
tion  and  incubating  overnight  at  37°C.  Peptides  were 
extracted  from  the  gel  slices  via  sonication  in  50  pi  60% 
acetonitrile/5%TFA,  dried  via  vacuum  centrifugation, 
and  reconstituted  in  15  pi  0.1%  TFA.  Tryptic  peptides 
were  desalted/enriched  using  a  C18  ZipTip  column 
(Millipore,  Billerica,  MA)  according  to  manufacturer’s 
instructions  and  the  eluant  was  spotted  on  a  MALDI 
plate  and  dried.  Samples  were  analyzed  using  a 
MALDI-LTQ  mass  spectrometer  (ThermoFinnigan,  San 
Jose,  CA).  A  full  MS  scan  in  high-mass  range  (m/z  600- 
4000,  5  microscans)  was  performed.  The  50  most 
intense  peaks  in  the  full  MS  spectrum  were  selected, 
and  MSMS  scans  were  performed  for  those  ions  in 
high-mass  range  (m/z  50-4000,  5  microscans),  the  nor¬ 
malized  collision  energy  for  MSMS  was  35.  Xcalibur 
software  was  used  to  process  the  mass  spectrometric 
data,  and  the  NCBInr  database  and  the  Bioworks  3.2 
search  engine  software  were  used  for  database 
searching. 
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Abstract 


Our  ability  to  rapidly  respond  to  emerging  influenza  is  hampered  somewhat  by  the 
lack  of  a  susceptible  small  animal  model.  To  develop  a  more  sensitive  model  we 
5  pathotyped  eighteen  low-pathogenic  non-mouse-adapted  influenza  A  viruses  of  human 
and  avian  origin  in  DBA/2  and  C57BL/6  mice.  The  majority  of  the  isolates  (13/18) 
induced  severe  morbidity  and  mortality  in  DBA/2  mice  upon  intranasal  challenge  with 
one  million  infectious  doses.  Also,  at  a  100-fold  lower  dose,  more  than  50%  of  the 
viruses  induced  severe  weight  loss  and  mice  succumbed  to  the  infection.  In  contrast,  only 
10  two  virus  strains  were  pathogenic  for  C57BL/6  mice  upon  high  dose  inoculation. 

Therefore  DBA/2  mice  are  a  suitable  model  to  validate  influenza  A  virus  vaccines  and 
antiviral  therapies  without  the  need  for  extensive  viral  adaptation.  Correspondingly,  we 
used  the  DBA/2  model  to  assess  the  level  of  protection  afforded  by  pre-existing 
pandemic  H1N1  2009  virus  (HINlpdm)  cross-reactive  human  antibodies  detected  by 
15  hemagglutination  inhibition  assay.  Passive  transfer  of  these  antibodies  prior  to  infection 
protected  mice  from  HINlpdm-induced  pathogenicity,  demonstrating  the  effectiveness 
of  these  cross-reactive  neutralizing  antibodies  in  vivo. 

Keywords:  Pandemic,  Influenza,  Animal  Model,  DBA/2,  Pandemic  H1N1  2009  virus 
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Introduction 


Respiratory  tract  infections  are  the  third  leading  cause  of  mortality  in  the  world  2S. 
Influenza,  a  disease  of  the  airways  caused  by  influenza  viruses,  is  responsible  for 
5  approximately  half  a  million  deaths  and  3-5  million  hospitalizations  per  year 29.  In 

addition  to  the  annual  disease  burden,  influenza  A  virus  is  more  notoriously  known  for  its 
ability  to  cause  pandemics.  Three  pandemics  have  been  reported  in  the  twentieth  century: 
the  first  that  occurred  in  1918  (Spanish  influenza)  killed  20-50  million  individuals  15;  the 
other  2  in  1957  and  1968,  although  less  lethal,  killed  millions  due  to  the  lack  of 
10  preexisting  immunity.  In  April  2009,  2  cases  of  febrile  illness  were  confirmed  to  be 
caused  by  swine-origin  influenza  A  virus  (H1N1) 3'7.  Continuous  spread  within  North 
America  and  other  parts  of  the  world  has  signaled  the  first  influenza  pandemic  of  this 
century. 

To  study  the  pathogenicity  of  influenza  A  viruses,  including  the  current  pandemic  A 
15  (H1N1)  2009  virus  (HINlpdm),  in  mammalian  hosts  and  to  determine  the  effectiveness 

of  pharmaceutical  interventions,  it  is  essential  to  have  a  sensitive  animal  model.  Although 
influenza  has  some  important  differences  in  mice  and  humans,  a  murine  model  is  the  only 
animal  model  thus  far  described  that  allows  for  relatively  high  group  numbers  and  any 
relatively  high  throughput.  Unfortunately,  only  few  strains  of  influenza  A  virus  -  almost 
20  exclusively  belonging  to  the  highly  pathogenic  avian  influenza  virus  isolates  of  the  H5 
and  H7  subtype  -  are  pathogenic  in  most  commonly  used  mouse  strains  without 
adaptation  through  serial  passaging.  The  hemagglutinin  proteins  of  these  H5  and  H7 
viruses  contain  a  basic  amino-acid  cleavage  site,  allowing  them  to  spread  systemically 
’  .  Most  other  subtypes  of  influenza  virus,  including  H1N1  and  H3N2,  either  do  not 
25  infect  or  cause  very  mild  disease  in  mice.  The  requirement  for  adaptation  of  a  pandemic 
virus  to  commonly  used  mouse  strains  can  lead  to  a  delay  in  the  gathering  of  important 
data  to  help  guide  public  health  control  strategies.  As  such,  the  lack  of  a  sensitive  small 
animal  model  to  study  the  infection  dynamics  of  various  subtypes  of  avian  influenza 
viruses  severely  hampers  the  rapid  and  effective  response  required  during  a  pandemic  or 
30  prepandemic  situation. 
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This  study  was  designed  to  demonstrate  the  utility  of  DBA/2  mice,  previously 
reported  to  be  susceptible  to  highly  pathogenic  influenza  viruses  ',  to  study  infections 
caused  by  several  influenza  A  virus  subtypes  isolated  from  birds  or  humans  without  the 
need  for  prior  adaptation.  To  assess  the  utility  of  the  model  to  respond  to  emerging 
5  strains,  we  used  DBA/2  mice  to  examine  the  functional  activity  of  sera  from  individuals 
previously  shown  to  have  pre-existing  cross-reactive  HINlpdm  antibodies.  It  is 
hypothesized  that  these  individuals  may  be  partially  protected  from  infection  because  of 
the  presence  of  cross-reactive  neutralizing  antibodies  produced  after  infection  with  a 
different  but  related  H1N1  virus.  This  hypothesis  is  supported  by  in  vitro 
10  microneutralization  and  hemagglutination  inhibition  assays  10'18;  however,  it  is  not  yet 
known  whether  these  antibodies  are  also  functional  in  vivo. 
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Material  and  Methods 


Mice  and  viruses 

Six-  to  ten-week-old  female  C57BL/6  and  DBA/2  mice  were  purchased  from  Jackson 
5  Laboratories  (Bar  Harbor,  ME)  and  housed  in  the  Animal  Resource  Center  at  St.  Jude 
Children’s  Research  Hospital  (St.  Jude).  The  mice  received  food  and  water  ad  libitum, 
and  all  experiments  were  conducted  in  accordance  with  rules  of  the  Institutional  Animal 
Control  and  Use  Committee  of  St.  Jude. 

Twenty-three  influenza  A  viruses  (Table  1)  from  9  different  hemagglutinin  subtypes 
10  (H1-H7,  H9  and  H10)  were  propagated  in  the  chorioallantoic  cavity  of  10-day-old 

embryonated  chicken  eggs.  The  allantoic  fluid  containing  infectious  particles  was 
harvested  48  h  after  inoculation  and  the  infectious  virus  titer  (egg  infectious  dose  50, 
EID50)  was  determined.  All  virus  stocks  had  a  minimum  titer  of  10s  5  EIDso/ml. 

15  Inoculation  of  mice  with  influenza  A  virus 

C57BL/6  and  DBA/2  mice  were  inoculated  with  influenza  A  viruses  intranasally  in 
30  pi  of  sterile  PBS  after  sedation  with  avertin  (2,2,2-tribromoethanol,  Sigma-Aldrich, 
MO,  USA).  The  mouse  lethal  dose  50  (MLD50)  was  determined  after  infecting  mice  with 
10-fold  serial  dilutions  of  the  viruses  from  106E1D5o  to  101  EID50.  Morbidity  and 
20  mortality  were  monitored  for  21  days  and  the  MLD50  values  were  calculated  by  the  Reed- 
Munch  method  21 .  Groups  of  5  mice  per  inoculum  size  per  isolate  were  tested  with  the 
exception  of  seasonal  H1N1  (106  EID50,  n  =  6  and  104  EID50,  n  =  8),  HINlpdm  (104 
EID50,  n  =  9),  H2N3  (106  EID50,  n  =  3  and  104  EID50,  n  =  4),  H2N4  (104  EID50,  n  =  4), 
H4N6  (106  EID50,  n  =  3),  H5N9  (10s  EID50,  n  =  6  and  104  EID50,  n  =  4),  H5N7  (10s 
25  EID50,  n  =  10),  H7N3  (106  EID50,  n  =  9  and  104  EID50,  n  =  4),  H7N9  (104  EID50,  n  =  4), 
H9N2/Y280  (10s  EID50,  n  =  10),  H9N5  (10s  and  104  EID50,  n  =  4),  H10N5  (10s  EID50,  n 
=  6  and  104  EID50,  n  =  8)  and  H10N7  (104  EID50,  n  =  4)  for  DBA/2  mice  and  H5N7  (106 
EID50,  n  =  8),  H6N5  (10s  EID50,  n  =  6),  H7N3  (10s  EID50,  n  =  10),  H7N9  (10s  EID50,  n  = 
4)  and  H9N2  (106  EID50,  n  =  4)  for  C57BL/6. 

30 
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Lung  viral  titers 

Lungs  were  collected  on  days  2  and  7  postinoculation  with  104EID5o  of  influenza  A 
virus  and  stored  at  -80°C.  They  were  homogenized  in  1.0  ml  of  minimal  essential 
medium,  and  homogenates  were  spun  for  5  min  at  1 000#  to  remove  cellular  debris.  The 
5  supernatant  was  used  to  quantify  the  amount  of  infectious  virus  present  in  the  lungs. 
Depending  on  the  virus  isolate,  viral  titers  were  determined  in  eggs  or  Madin-Darby 
Canine  Kidney  (MDCK)  cells  as  described  previously  '. 

Hemagglutination  inhibition  and  virus  neutralization  assays 
10  Influenza  A  virus  neutralizing  activity  in  serum  was  quantified  by  hemagglutination 

inhibition  (HI)  and  virus  microneutralization  (VN)  assay.  Sera  were  first  treated  with 
receptor  destroying  enzyme  (RDE)  (RDE  (II)  “Seiken”,  Denka  Seiken  UK  Ltd,  UK)  for 
18  h  at  37°C  followed  by  30  minutes  inactivation  at  56°C.  HI  assays  were  done  with  4 
hemagglutination  units  of  the  virus  and  0.5%  turkey  red  blood  cells  (HINlpdm)  or  0.5% 
15  chicken  red  blood  cells  (avian  virus  isolates)  as  described  previously  9.  For  a  VN  assay 
the  sera  were  diluted  two-fold  starting  at  1:10  dilution  in  PBS  and  incubated  for  lh  at 
37°C  with  100  TCTDso  of  A/California/4/09  virus.  Next,  100pl  of  the  mixture  of  virus  and 
sera  were  added  to  MDCK  cells  for  1  h  at  37°C.  Following  the  aspiration  of  the 
supernatant,  cells  were  washed  with  PBS  and  fresh  200pl  of  minimal  essential  medium 
20  supplemented  with  0.1%  bovine  serum  albumin  (A8412,  Sigma- Aldrich),  antibiotics 
(Invitrogen,  NY,  USA),  vitamins  (Invitrogen)  and  lpg/ml  TPCK  trypsin  (Worthington, 
NJ,  USA)  was  added.  After  3-4  days  at  37°C,  the  assay  was  developed  by  HA  assay 
using  turkey  red  blood  cells.  The  average  HI  and  VN  titer  was  calculated  following  log2 
transformation  of  the  highest  serum  dilution  able  to  inhibit  hemagglutination  or  virus 
25  replication  respectively. 

Passive  antibody  transfer 

Human  sera  were  collected  as  part  of  a  clinical  trial  conducted  during  the  2007-08 
and  2008-09  influenza  seasons  in  the  Greater  Vancouver  Area  of  British  Columbia, 

30  Canada  or  in  the  vicinity  of  the  Greater  Hartford  Area  of  Connecticut.  All  participants 
received  the  standard  dose  of  the  licensed  trivalent  split-virus  influenza  vaccine 


6 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


containing  A/Solomon  Islands/3/2006-like  (H1N1),  AAVisconsin/67/2005-like  (H3N2), 
and  B/Malaysia/2506/2004-like  viruses  in  2007-08  or  A/Brisbane/59/2007  (HlNl)-like, 
A/Brisbane/10/2007  (H3N2)-like,  and  B/Florida/4/2006-like  viruses  in  2008-09.  Sera 
were  collected  before  vaccination  and  4  weeks  after  vaccination.  Post-vaccination  sera 
5  from  individuals,  aged  65  years  and  older,  with  a  detectable  HI  and  VN  titer  toward 
HINlpdm  (A/California/4/2009)  were  pooled  and  heat-inactivated  for  30  min  at  56°C. 

To  study  the  effect  of  neutralizing  antibodies  we  used  age-matched  pooled  human  sera 
without  a  detectable  HI  and  VN  titer  to  HINlpdm  (A/Califomia/4/2009),  seasonal  H1N1 
(A/Brisbane/59/2007)  and  H7N3  (A/shorebird/Delaware/22/2006)  vims.  Ferret 
10  polyclonal  serum  obtained  from  ferrets  14  days  after  inoculation  with  the  HINlpdm  vims 
(HI  titer  of  2560,  VN  titer  of  320)  or  PBS  were  used  as  a  positive  and  negative  control 
respectively.  Four  hundred  micro  liters  of  pooled  human  sera,  diluted  1 : 1  in  PBS,  were 
injected  intraperitoneally  into  10  mice  24  h  prior  to  inoculation  with  a  lethal  dose  of 
virus.  The  positive  and  negative  controls  were  also  injected  into  10  mice  each  for  the 
15  HINlpdm  experiment,  while  5  PBS  control  mice  were  included  in  the  H7N3  and 
seasonal  H1N1  follow-up  experiment. 

Cytokine  analysis 

Lungs  were  collected  on  days  2  and  7  postinoculation  with  104EID5o  of  influenza  A 
20  virus  and  concentrations  of  CCL2,  CCL5,  interleukin-6  (IL-6),  tumor  necrosis  factor- 

alpha  (TNF-a)  and  interferon-gamma  (IFN-y)  were  determined  as  described  previously  . 
ELISA  was  performed  according  to  manufacturer’ s  instructions  (Quantikine  kits,  R&D 
Systems,  Minneapolis,  MN).  At  least  4  animals  infected  with  a  particular  strain  of 
influenza  A  virus  were  tested  for  1  cytokine  at  a  given  time  point. 

25 

Statistical  analysis 

Statistical  analyses  of  differences  in  mortality  were  determined  by  using  the  log-rank 
test.  The  Student’s  T-test  was  used  to  analyze  differences  in  lung  virus  titers  between  the 
different  strains  of  mice  following  ln-transformation  of  the  data  as  well  as  to  determine 
30  statistical  significance  in  cytokine  and  chemokine  production  and  weight  loss  after 
influenza  A  virus  infection. 
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Results 


Increased  susceptibility  of  DBA/2  mice  to  influenza  A  virus  isolates 

To  assess  the  utility  of  DBA/2  mice  as  a  more  universal  small  animal  model  for 
5  influenza,  we  tested  a  range  of  different  viral  subtypes  for  their  ability  to  induce 

morbidity  and  mortality  in  this  host.  A  set  of  twenty-three  viruses  belonging  to  9  different 
hemagglutinin  subtypes  was  selected  and  used  to  inoculate  DBA/2  mice.  At  a  dose  of  106 
EID50,  18  or  78%  of  the  viruses  were  pathogenic  and  mice  succumbed  to  infection  5-12 
days  postinoculation,  depending  on  the  virus  isolate  (Table  1).  Inoculation  with  a  lower 
10  dose  of  virus  (104  EID50)  caused  severe  weight  loss  and  death  of  DBA/2  mice  in  14  of  22 
(64%)  virus  isolates  tested  (Table  1).  These  isolates  included  a  seasonal  human  H1N1 
virus  from  2008  and  both  H10  viruses.  The  ability  to  induce  severe  disease,  as  measured 
by  weight  loss  and  mortality,  in  DBA/2  mice  was  not  limited  to  certain  subtypes  of 
influenza;  however,  virus  isolates  of  the  H2  and  H4  subtype  were  only  mildly  pathogenic. 
15  Few  isolates  were  pathogenic  at  102  EID50,  and  these  included  the  2  mouse-adapted 

influenza  A  viruses  (A/Puerto  Rico/8/1934  (100%  mortality)  and  X31  (60%  mortality)),  2 
highly  pathogenic  H5N1  influenza  A  viruses  (A/Hong  Kong/213/2003  (84%  mortality) 
and  A/Vietnam/ 1203/2004  (100%  mortality)),  a  low  pathogenic  H7N3 
(A/shorebird/Delaware/22/2006,  45%  mortality)  virus,  and  the  HINlpdm  (100% 

20  mortality).  Eighteen  isolates  were  also  tested  in  C57BL/6  mice  at  a  dose  of  106  EID50. 

Only  the  mouse-adapted  virus  A/Puerto  Rico/8/1934  virus,  3  highly  pathogenic  viruses,  a 
H6N1  virus  (A/teal/Hong  Kong/W312/1997),  and  the  HINlpdm  virus  caused  severe 
disease  and  mortality  (Table  1). 

25  Nonadapted  avian  influenza  A  viruses  can  replicate  to  high  titers  in  C57BL/6  and 
DBA/2  mice 

We  hypothesize  that  the  large  difference  in  pathogenicity  of  avian  influenza  A  viruses 
between  DBA/2  and  C57BL/6  mice  may  be  due  to  increased  replication  efficiency.  To 
test  this,  we  measured  day  2  and  7  postinoculation  pulmonary  viral  loads  in  mice  infected 
30  with  104  EID50  of  H7N3  or  H10N5  virus  (Figure  1).  These  avian  virus  isolates  were 

selected  for  their  exceptionally  large  difference  in  pathology  score  between  DBA/2  and 
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C57BL/6  mice.  High  dose  inoculation  (106  EID50)  with  H7N3  or  H10N5  virus  induced 
7%  and  4%  maximum  weight  loss  in  C57BL/6  mice,  while  DBA/2  mice  succumbed  to 
infection  with  102  or  104  EIDgo  respectively.  On  day  2  postinoculation,  lung  viral  titers  of 
DBA/2  and  C57BL/6  mice  were  similar  for  both  H7N3  (106  25  vs.  106  75  EID50/ml, 

5  respectively)  and  H10N5  virus  (105  3  vs.  106  0  EIDgo/ml,  respectively,  Figure  1). 

Interestingly,  titers  in  H7N3-infected  lungs  were  approximately  1-logio  higher  than  those 
of  HI 0N5 -infected  lungs.  Day  7  postinoculation,  H7N3-infected  lungs  of  DBA/2  mice 
had  higher  virus  loads  than  that  of  C57BL/6  mice  (P  <  0.05,  Figure  1).  Lungs  of  H10N5- 
infected  DBA/2  mice  also  contained  more  virus  than  those  of  C57BL/6  mice,  however 
10  this  difference  was  not  significant  (P  =  0.07). 

To  assess  the  ability  of  other  non-  or  low-pathogenic  avian  influenza  A  virus  isolates 
to  infect  mice,  we  measured  the  serological  response  against  the  challenge  virus  in 
convalescent  sera  as  a  surrogate  marker  for  viral  replication.  In  C57BL/6  mice,  an  HI  titer 
was  detected  postinoculation  with  5  of  the  8  virus  isolates  studied,  suggesting  that  the 
15  majority  of  the  isolates  replicate  in  the  respiratory  tract  of  mice  (Table  2).  Inoculation  of 
3  viruses,  an  H4N1,  H6N5  and  an  H5N7  virus  did  not  cause  seroconversion  in  C57BL/6 
mice.  Because  DBA/2  mice  are  generally  more  susceptible  to  influenza  virus  infection 
and  fewer  isolates  were  non-  or  low-pathogenic,  we  tested  convalescent  sera  following 
inoculation  with  only  5  virus  isolates.  Of  the  5  convalescent  sera  tested,  3  contained  a 
20  detectable  HI  titer  whereas  2  (an  H4N6  and  an  H9N6  virus)  did  not  (Table  2).  This  data 
suggests  that  most  influenza  virus  isolates  are  capable  of  replicating  in  the  respiratory 
tract  of  mice,  but  the  outcome  after  infection  depends  entirely  on  the  mouse  strain,  virus 
strain,  or  a  combination  of  both. 

25  Pandemic  H1N1  2009  virus  A/California/4/2009  is  highly  pathogenic  in  DBA/2  mice 
Based  on  the  results  presented  above,  we  next  looked  at  the  replication  of  the 
HINlpdm  viruses  in  DBA/2  and  C57BL/6  mice.  Inoculation  of  mice  with  106— 102  EID50 
of  A/California/4/09,  a  representative  HINlpdm  virus,  resulted  in  100%  mortality  after 
8-12  days,  and  that  with  101  EID50  caused  significant  weight  loss  in  100%  and  mortality 
30  in  50%  of  DBA/2  mice  (Table  1).  In  contrast,  C57BL/6  mice  lost  a  significant  amount  of 
weight  (14%,  P<0.05)  by  day  7  when  inoculated  with  106  EID50  and  60%  of  the  mice 
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died.  Inoculation  with  105  to  104  EID50  did  not  cause  death  of  C57BL/6  mice.  Therefore, 
the  MLD50  for  A/California/4/2009  was  105-fold  lower  in  DBA/2  than  C57BL/6  mice,  a 
finding  consistent  with  other  viral  strains. 

Increased  pathogenicity  is  often  associated  with  higher  viral  loads  and  increased 
5  levels  of  proinflammatory  cytokines  such  as  CCL2,  IL6,  and  TNFa.  At  Day  2 
postinoculation  with  104  EID50  of  HINlpdm,  lung  viral  titers  of  DBA/2  mice  and 
C57BL/6  mice  were  similar  (Figure  la).  There  were  also  no  significant  differences  in  the 
levels  of  CCL2,  CCL5,  and  IL6  between  both  strains  (Table  3),  but  levels  of  TNFa  were 
significantly  higher  in  DBA/2  than  C57BL/6  mice  (P  <  0.01).  At  Day  7  postinfection, 

10  lung  homogenates  of  DBA/2  mice  had  higher  viral  loads  than  C57BL/6  mice  (P  <  0.01, 
Figure  2b)  as  well  as  significantly  higher  concentrations  of  CCL2  and  IL6  (P  <  0.01, 
Table  3). 

Passive  transfer  of  cross-reactive  human  neutralizing  antibodies  protects  against 
15  HINlpdm-induced  pathogenicity 

To  assess  the  functionality  of  human  cross-reactive  polyclonal  antibodies,  passive 
antibody  transfer  experiments  were  performed  using  the  highly  susceptible  DBA/2  mice. 
Human  sera  were  obtained  from  individuals  not  previously  exposed  to  the  HINlpdm 
virus  but  who  had  cross-reactive  neutralizing  antibodies  to  it  with  HI  titer  of  160  and  VN 
20  titer  of  95.  Twenty-four  hours  after  passive  transfer,  mice  were  inoculated  with  102  EID50 
(3-5  LD50)  of  A/California/4/2009  virus  and  morbidity  and  mortality  were  assessed.  As 
expected,  all  control  mice  treated  with  PBS  died  within  10  days  of  inoculation  (Figure 
2a).  Mice  treated  with  ferret  sera  containing  high  levels  of  HINlpdm-specific 
neutralizing  antibodies  did  not  succumb  to  infection  (P  <  0.001)  and  lost  significantly 
25  less  weight  on  days  7,  10,  13,  and  15  than  mice  in  all  other  groups  (P  <  0.01,  Figure  2). 
DBA/2  mice  injected  with  pooled  human  serum  containing  cross-reactive  antibodies  had 
a  higher  survival  rate  (75%)  than  mice  in  the  PBS  control  group  (0%,  P  <  0.01)  and  mice 
receiving  human  sera  without  detectable  cross-neutralizing  antibodies  (15%,  P  <  0.05). 
Increased  survival  of  the  mice  injected  with  cross-neutralizing  antibodies  was 
30  accompanied  by  a  decrease  in  percent  weight  loss  on  days  10(11%)  and  13  (20%) 
postinfection  (P  <  0.01,  Figure  2b).  These  data  indicate  that  HINlpdm-specific 
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neutralizing  antibodies  induced  after  infection  with  a  related  H1N1  virus  can  protect 
DBA/2  mice  from  a  lethal  challenge  and  are  likely  responsible  for  the  age  related  attack 
rates  seen  in  humans.  The  protective  effect  of  human  sera  without  any  neutralizing 
antibodies  was  validated  using  two  additional  challenge  models,  a  seasonal  H1N1 
5  (A/Memphis/33/2008)  and  an  avian  H7N3  (A/shorebird/Delaware/22/2006)  virus  isolate. 

The  survival  rates  increased  significantly  for  the  seasonal  H1N1  virus  (90%,  P  <  0.01) 
but  not  for  H7N3  (17%,  P  >  0.05)  suggesting  that  in  vitro  HI  or  MN  assays  underestimate 
the  levels  of  pre-existing  neutralizing  immunity  to  influenza  A  virus  strain  in  humans. 
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Discussion 

The  present  study  establishes  that  DBA/2  mice  are  very  susceptible  to  most  influenza 
A  virus  isolates  and  that  infection  often  results  in  debilitating  pneumonia  and  subsequent 
death.  This  sensitivity  of  DBA/2  mice  was  used  to  demonstrate  that  HINlpdm  is  more 
5  pathogenic  than  the  circulating  seasonal  H1N1  viruses  and  that  pre-existing  human  cross¬ 
reactive  neutralizing  antibodies  can  prevent  HINlpdm-induced  mortality  and  morbidity. 

Small  animal  models,  like  the  mouse,  have  frequently  been  used  for  influenza  virus 
research  including  areas  such  as  pathogenesis,  vaccine  efficacy,  and  antiviral  therapies. 
The  preferred  strains,  C57BL/6  and  Balb/c,  display  few  clinical  symptoms  upon  high 
10  dose  inoculation  with  most  influenza  A  virus  isolates  and  only  highly  pathogenic  or 

mouse-adapted  viruses  cause  severe  morbidity  and  mortality  at  low  doses.  As  such,  it  was 
generally  believed  that  mice  are  resistant  to  most  human  and  avian  influenza  A  viruses. 
The  current  study,  as  well  as  a  recent  report  by  Driksell  el  al 5,  provide  substantial 
evidence  to  suggest  that  Balb/c  and  C57BL/6  mice  are  susceptible  to  infection  with  many 
15  different  influenza  A  virus  isolates  (H1-H7,  H9-H1 1);  albeit  that  the  infection  does  not 
cause  significant  disease.  In  contrast,  DBA/2  mice  become  sick  and  often  succumb  to 
infection  with  the  majority  of  the  tested  isolates.  This  enhanced  susceptibility  of  DBA/2 
mice  was  previously  reported  for  highly  pathogenic  and  mouse-adapted  influenza  viruses 
124  and  now  includes  many  non-adapted  avian  influenza  isolates. 

20  The  mechanism  for  the  difference  in  susceptibility  between  C57BL/6  and  DBA/2  is 
not  yet  fully  understood  but  involves  multiple  genetic  differences  between  the  two  mouse 
lineages  affecting  several  pathways  and  processes  1 .  Certain  influenza  viruses  grow  to 
higher  titers  in  DBA/2  mice  (A/Hong  Kong/213/2003  (H5N1)  or  A/Memphis/33/2008 
(H1N1,  data  not  shown)),  while  others  do  not  (H7N3  and  H10N5,  this  study). 

25  Irrespective  of  the  difference  in  viral  load,  DBA/2  mice  respond  more  vigorously 

producing  larger  quantities  of  certain  pro-inflammatory  molecules  like  TNF-a  which  was 
shown  to  correlate  with  increased  morbidity  and  mortality  in  humans  4.  Previous  work 
has  shown  that  the  lack  of  a  functional  He  in  DBA/2  mice  may  result  in  a  less  effective 
adaptive  immune  response  increasing  the  viral  load  during  later  stages  of  the  infection 
30  1,n.  In  the  same  study,  histopathologic  evaluation  DBA/2  and  C57BL/6  mice  after  H5N1 

virus  infection  was  able  to  clearly  show  more  extensive  involvement  of  the  entire  lung  as 
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well  as  necrotic  epithelium  on  day  4  in  the  DBA/2  mice.  At  day  7  there  was  a  dramatic 
difference  in  the  number  of  cells  that  stain  positive  for  influenza  between  the  two  strains. 
In  the  C57BL/6  mice,  an  inflammatory  response  was  observed  in  the  parenchyma  with 
multifocal  macrophages  within  the  infiltrate  that  stain  positive  for  influenza.  In  contrast, 

5  the  DBA/2  mice  maintain  widespread  intensely  positive  staining  of  the  airway 

epithelium,  alveolar  macrophages.  Type  I  and  Type  II  cells  in  the  alveolar  wall.  Overall, 
the  combined  effect  of  poor  clearance,  excessive  inflammation,  and  elevated  virus  titers 
likely  creates  a  highly  pathogenic  environment  in  DBA/2  mice. 

The  exceptional  susceptibility  of  DBA/2  mice  for  influenza  virus  infections  was  used 
10  to  demonstrate  the  effectiveness  of  HINlpdm  cross-neutralizing  human  antibodies  in 
reducing  HINlpdm  driven  mortality.  Approximately  50%  of  the  US  population  over  75 
years  of  age  have  low  levels  of  pre-exposure  cross-reactive  antibodies  to  HINlpdm  13. 
The  age  of  these  individuals  suggests  that  these  antibodies  were  generated  between  1918 
and  1930  when  an  antigenically  similar  H1N1  virus  was  possibly  circulating  among 
15  humans.  We  show  in  this  study  that  these  cross-reactive  antibodies  are  also  fully 

functional  in  vivo ,  and  may  offer  some  degree  of  protection  to  an  otherwise  at-risk  elderly 
population  against  the  current  HINlpdm.  The  identity  of  the  neutralizing  epitopes  has  yet 
to  be  identified  but  could  be  located  within  the  antigenic  site  Sa  J°.  Interestingly,  there 
was  a  small  but  significant  level  of  protection  when  mice  were  injected  with  human 
20  serum  without  detectable  levels  of  virus  neutralizing  antibodies.  The  protection  was 

specific  for  H1N1  viruses  since  the  transfer  of  HI  negative  sera  did  not  affect  the  survival 
after  challenge  with  a  H7N3  virus.  The  protection  is  possibly  mediated  by  antibodies 
specific  for  the  M2  protein  19,25,26  or  to  non-neutralizing  epitopes  on  the  HA  or  NA 

16  22  2  14  23 

’  ’  '  ’  .  These  data  also  suggest  that  the  current  in  vitro  assays  (HI  and  VN)  are 
25  underestimating  the  level  of  pre-existing  protective  immunity  in  the  human  population. 
Addition  of  serum  factors,  like  Clq,  have  been  shown  to  increase  the  sensitivity  of  these 
assays  6'17.  The  limited  capacity  of  HI  or  VN  assays  to  detect  neutralizing  antibodies  was 
previously  noted  in  pre-clinical  H5N 1  ferret  vaccine  studies.  A  single  dose  of  inactivated 
H5N 1  vaccine  did  not  induce  a  detectable  HI  or  VN  titer,  however  the  animals  were 
30  protected  from  a  lethal  challenge  with  highly  pathogenic  H5N 1  virus  8. 
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Although  the  DBA/2  mouse  model  will  provide  a  useful  analytical  tool  to  study 
viruses  and  antiviral  therapies,  the  model  may  not  necessarily  reflect  the  natural  response 
to  influenza  viruses  in  humans.  With  the  exception  of  highly  pathogenic  H5N1  viruses, 
one  should  be  careful  interpreting  pathogenesis  data  upon  infection  with  these  avian 
5  isolates  as  described  in  the  current  study. 

To  summarize,  we  have  confirmed  that  the  DBA/2  mouse  model  is  a  suitable  and 
highly  susceptible  animal  model  to  study  infection  of  influenza  A  viruses  of  various 
subtypes,  including  those  previously  known  to  infect  humans.  Also,  this  model  will  allow 
us  to  define  the  requirements  of  viruses  of  avian  origin  to  infect  mammalian  hosts  and 
10  rapidly  evaluate  vaccines  or  antiviral  therapies  in  the  event  of  a  pandemic  emergency. 


14 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


Acknowledgements 


We  thank  David  Walker  for  isolating  and  propagating  influenza  A  virus  isolates.  We  also 
acknowledge  Drs.  Nancy  Cox,  Sasha  Klimov,  and  Ruben  Donis  (Centers  of  Disease 
5  Control),  Dr.  Ron  Fouchier  (Erasmus  Medical  Center),  Drs.  Malik  Peiris  and  Guan  Yi 
(University  of  Hong  Kong),  and  the  World  Health  Organization  Global  Influenza 
Surveillance  Network  as  a  source  of  influenza  A  viruses.  Finally,  we  thank  Drs.  M 
Ducatez  and  S.  Schultz-Cherry  for  critically  reviewing  the  paper.  This  project  was  funded 
in  part  by  grants  from  the  National  Institute  Allergy  and  Infectious  Diseases,  National 
10  Institutes  of  Health,  Department  of  Health  and  Human  Services  under  contract  no. 

HHSN266200700005C,  by  the  Department  of  Defense  award  W81XWH-09-1-0391,  by  the 
Centers  of  Infectious  Diseases  Control  at  St  Jude  Children’s  Research  Hospital  and  the 
American  Lebanese  Syrian  Associated  Charities  (ALSAC). 


15 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


Reference  List 


1 .  Boon,  A.  J.  Debeauchamp,  A.  Hollmann,  J.  Luke,  M.  Kotb,  S.  Rowe,  D. 
Finkelstein,  G.  Neale,  L.  Lu,  R.  W.  Williams,  and  R.  J.  Webby.  2009.  Host 
Genetic  Variation  Affects  Resistance  to  Infection  with  a  Highly  Pathogenic  H5N1 
Influenza  A  Virus  in  Mice.  J.  Virol.  83:10417-26 

2.  Chen,  Z.,  S.  Kadowaki,  Y.  Hagiwara,  T.  Yoshikawa,  K.  Matsuo,  T.  Kurata, 
and  S.  Tamura.  2000.  Cross-protection  against  a  lethal  influenza  virus  infection  by 
DNA  vaccine  to  neuraminidase.  Vaccine.  18:3214-3222. 

3.  Dawood,  F.  S.,  S.  Jain,  L.  Finelli,  M.  W.  Shaw,  S.  Lindstrom,  R.  J.  Garten,  L. 

V.  Gubareva,  X.  Xu,  C.  B.  Bridges,  and  T.  M.  Uyeki.  2009.  Emergence  of  a 
novel  swine-origin  influenza  A  (H1N1)  virus  in  humans.  N.  Engl.  J.  Med. 
360:2605-2615. 

4.  de  Jong,  M.  D.,  C.  P.  Simmons,  T.  T.  Thanh,  V.  M.  Hien,  G.  J.  Smith,  T.  N. 
Chau,  D.  M.  Hoang,  N.  V.  Chau,  T.  H.  Khanh,  V.  C.  Dong,  P.  T.  Qui,  B.  V. 
Cam,  d.  Q.  Ha,  Y.  Guan,  J.  S.  Peiris,  N.  T.  Chinh,  T.  T.  Hien,  and  J.  Farrar. 

2006.  Fatal  outcome  of  human  influenza  A  (H5N1)  is  associated  with  high  viral 
load  and  hypercytokinemia.  Nat.  Med.  12:1203-1207. 

5.  Driskell,  E.  A.,  C.  A.  Jones,  D.  E.  Stallknecht,  E.  W.  Howerth,  and  S.  M. 
Tompkins.  2010.  Avian  influenza  virus  isolates  from  wild  birds  replicate  and  cause 
disease  in  a  mouse  model  of  infection.  Virology.  399:280-289. 

6.  Feng,  J.  Q.,  K.  Mozdzanowska,  and  W.  Gerhard.  2002.  Complement  component 
Clq  enhances  the  biological  activity  of  influenza  virus  hemagglutinin-specific 
antibodies  depending  on  their  fine  antigen  specificity  and  heavy-chain  isotype.  J. 
Virol.  76:1369-1378. 

7.  Garten,  R.  J.,  C.  T.  Davis,  C.  A.  Russell,  B.  Shu,  S.  Lindstrom,  A.  Balish,  W. 

M.  Sessions,  X.  Xu,  E.  Skepner,  V.  Deyde,  M.  Okomo-Adhiambo,  L.  Gubareva, 
J.  Barnes,  C.  B.  Smith,  S.  L.  Emery,  M.  J.  Hillman,  P.  Rivailler,  J.  Smagala,  G. 
M.  de,  D.  F.  Burke,  R.  A.  Fouchier,  C.  Pappas,  C.  M.  puche- Aranda,  H.  Lopez- 
Gatell,  H.  Olivera,  I.  Lopez,  C.  A.  Myers,  D.  Faix,  P.  J.  Blair,  C.  Yu,  K.  M. 
Keene,  P.  D.  Dotson,  Jr.,  D.  Boxrud,  A.  R.  Sambol,  S.  H.  Abid,  G.  K.  St,  T. 
Bannerman,  A.  L.  Moore,  D.  J.  Stringer,  P.  Blevins,  G.  J.  mmler-Harrison,  M. 
Ginsberg,  P.  Kriner,  S.  Waterman,  S.  Smole,  H.  F.  Guevara,  E.  A.  Belongia,  P. 
A.  Clark,  S.  T.  Beatrice,  R.  Donis,  J.  Katz,  L.  Finelli,  C.  B.  Bridges,  M.  Shaw, 
D.  B.  Jernigan,  T.  M.  Uyeki,  D.  J.  Smith,  A.  I.  Klimov,  and  N.  J.  Cox.  2009. 
Antigenic  and  genetic  characteristics  of  swine-origin  2009  A(H1N1)  influenza 
viruses  circulating  in  humans.  Science.  325:197-201. 

8.  Govorkova,  E.  A.,  R.  J.  Webby,  J.  Humberd,  J.  P.  Seiler,  and  R.  G.  Webster. 

2006.  Immunization  with  reverse-genetics-produced  H5N1  influenza  vaccine 


16 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


protects  ferrets  against  homologous  and  heterologous  challenge.  J.  Infect.  Dis. 
194:159-167. 

9.  Hancock,  K.,  Y.  Yeguilla,  X.  Lu,  W.  Zhong,  E.  N.  Butler,  H.  Sun,  F.  Liu,  L. 
Dong,  J.  R.  Devos,  P.  M.  Gargiullo,  T.  L.  Brammer,  N.  J.  Cox,  T.  M.  Tumpey, 
and  J.  M.  Katz.  2009.  Cross-Reactive  Antibody  Responses  to  the  2009  Pandemic 
H1N1  Influenza  Virus.  N.  Engl.  J.  Med. 

10.  Hancock,  K.,  V.  Veguilla,  X.  Lu,  W.  Zhong,  E.  N.  Butler,  H.  Sun,  F.  Liu,  L. 
Dong,  J.  R.  Devos,  P.  M.  Gargiullo,  T.  L.  Brammer,  N.  J.  Cox,  T.  M.  Tumpey, 
and  J.  M.  Katz.  2009.  Cross-Reactive  Antibody  Responses  to  the  2009  Pandemic 
H1N1  Influenza  Virus.  N.  Engl.  J.  Med. 

1 1 .  Hicks,  J.  T.,  F.  A.  Ennis,  E.  Kim,  and  M.  Verbonitz.  1978.  The  importance  of  an 
intact  complement  pathway  in  recovery  from  a  primary  viral  infection:  influenza  in 
decomplemented  and  in  C5-deficient  mice.  J.  Immunol.  121:1437-1445. 

12.  Horimoto,  T.  and  Y.  Kawaoka.  1994.  Reverse  genetics  provides  direct  evidence 
for  a  correlation  of  hemagglutinin  cleavability  and  virulence  of  an  avian  influenza 
A  virus.  J.  Virol.  68:3120-3128. 

13.  Itoh,  Y.,  K.  Shinya,  M.  Kiso,  T.  Watanabe,  Y.  Sakoda,  M.  Hatta,  Y. 

Muramoto,  D.  Tamura,  Y.  Sakai-Tagawa,  T.  Noda,  S.  Sakabe,  M.  Imai,  Y. 
Hatta,  S.  Watanabe,  C.  Li,  S.  Yamada,  K.  Fujii,  S.  Murakami,  H.  Imai,  S. 
Kakugawa,  M.  Ito,  R.  Takano,  K.  Iwatsuki-Horimoto,  M.  Shimojima,  T. 
Horimoto,  H.  Goto,  K.  Takahashi,  A.  Makino,  H.  Ishigaki,  M.  Nakayama,  M. 
Okamatsu,  K.  Takahashi,  D.  Warshauer,  P.  A.  Shult,  R.  Saito,  H.  Suzuki,  Y. 
Furuta,  M.  Yamashita,  K.  Mitamura,  K.  Nakano,  M.  Nakamura,  R. 
Brockman-Schneider,  H.  Mitamura,  M.  Yamazaki,  N.  Sugaya,  M.  Suresh,  M. 
Ozawa,  G.  Neumann,  J.  Gern,  H.  Kida,  K.  Ogasawara,  and  Y.  Kawaoka.  2009. 
In  vitro  and  in  vivo  characterization  of  new  swine-origin  H1N1  influenza  viruses. 
Nature.  20;460:102 1-1025. 

14.  Johansson,  B.  E.,  D.  J.  Bucher,  and  E.  D.  Kilbourne.  1989.  Purified  influenza 
virus  hemagglutinin  and  neuraminidase  are  equivalent  in  stimulation  of  antibody 
response  but  induce  contrasting  types  of  immunity  to  infection.  J.  Virol.  63:1239- 
1246. 

15.  Johnson,  N.  P.  and  J.  Mueller.  2002.  Updating  the  accounts:  global  mortality  of 
the  1918-1920  "Spanish"  influenza  pandemic.  Bull.  Hist  Med.  76:105-115. 

16.  McLain,  L.  and  N.  J.  Dimmock.  1989.  Protection  of  mice  from  lethal  influenza  by 
adoptive  transfer  of  non-neutralizing  haemagglutination-inhibiting  IgG  obtained 
from  the  lungs  of  infected  animals  treated  with  defective  interfering  virus.  J.  Gen. 
Virol.  70:2615-2624. 

17.  Mehlhop,  E.,  S.  Nelson,  C.  A.  Jost,  S.  Gorlatov,  S.  Johnson,  D.  H.  Fremont,  M. 
S.  Diamond,  and  T.  C.  Pierson.  2009.  Complement  protein  Clq  reduces  the 


17 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


stoichiometric  threshold  for  antibody-mediated  neutralization  of  West  Nile  virus. 
Cell  Host.  Microbe.  6:381-391. 

18.  MMWR.  2009.  Serum  cross-reactive  antibody  response  to  a  novel  influenza  A 
(H1N1)  virus  after  vaccination  with  seasonal  influenza  vaccine.  MMWR  Morb. 
Mortal.  Wkly.  Rep.  58:521-524. 

19.  Mozdzanowska,  K.,  K.  Maiese,  M.  Furchner,  and  W.  Gerhard.  1999.  Treatment 
of  influenza  virus-infected  SCID  mice  with  nonneutralizing  antibodies  specific  for 
the  transmembrane  proteins  matrix  2  and  neuraminidase  reduces  the  pulmonary 
virus  titer  but  fails  to  clear  the  infection.  Virology.  254:138-146. 

20.  Ohuchi,  M.,  M.  Orlich,  R.  Ohuchi,  B.  E.  Simpson,  W.  Garten,  H.  D.  Klenk, 
and  R.  Rott.  1989.  Mutations  at  the  cleavage  site  of  the  hemagglutinin  after  the 
pathogenicity  of  influenza  virus  A/chick/Penn/83  (H5N2).  Virology.  168:274-280. 

21 .  Reed,  L.  J.  and  H.  Muench.  1938.  A  simple  method  for  estimating  fifty  percent 
endpoints.  Am.  J.  Hyg.  27:493-497. 

22.  Sambhara,  S.,  A.  Kurichh,  R.  Miranda,  T.  Tumpey,  T.  Rowe,  M.  Renshaw,  R. 
Arpino,  A.  Tamane,  A.  Kandil,  O.  James,  B.  Underdown,  M.  Klein,  J.  Katz, 
and  D.  Burt.  2001.  Heterosubtypic  immunity  against  human  influenza  A  viruses, 
including  recently  emerged  avian  H5  and  H9  viruses,  induced  by  FLU-ISCOM 
vaccine  in  mice  requires  both  cytotoxic  T-lymphocyte  and  macrophage  function. 
Cell  Immunol.  211:143-153. 

23.  Sandbulte,  M.  R.,  G.  S.  Jimenez,  A.  C.  Boon,  L.  R.  Smith,  J.  J.  Treanor,  and  R. 
J.  Webby.  2007.  Cross-reactive  neuraminidase  antibodies  afford  partial  protection 
against  H5N1  in  mice  and  are  present  in  unexposed  humans.  PLoS.  Med.  4:e59. 

24.  Srivastava,  B.,  P.  Blazejewska,  M.  Hessmann,  D.  Bruder,  R.  Geffers,  S.  Mauel, 
A.  D.  Gruber,  and  K.  Schughart.  2009.  Host  genetic  background  strongly 
influences  the  response  to  influenza  a  virus  infections.  PLoS.  One.  4:e4857. 

25.  Tompkins,  S.  M.,  Z.  S.  Zhao,  C.  Y.  Lo,  J.  A.  Misplon,  T.  Liu,  Z.  Ye,  R.  J. 
Hogan,  Z.  Wu,  K.  A.  Benton,  T.  M.  Tumpey,  and  S.  L.  Epstein.  2007.  Matrix 
protein  2  vaccination  and  protection  against  influenza  viruses,  including  subtype 
H5N1.  Emerg.  Infect.  Dis.  13:426-435. 

26.  Treanor,  J.  J.,  E.  L.  Tierney,  S.  L.  Zebedee,  R.  A.  Lamb,  and  B.  R.  Murphy. 

1990.  Passively  transferred  monoclonal  antibody  to  the  M2  protein  inhibits 
influenza  A  virus  replication  in  mice.  J.  Virol.  64:1375-1377. 

27.  Vey,  M.,  M.  Orlich,  S.  Adler,  H.  D.  Klenk,  R.  Rott,  and  W.  Garten.  1992. 
Hemagglutinin  activation  of  pathogenic  avian  influenza  viruses  of  serotype  H7 
requires  the  protease  recognition  motif  R-X-K/R-R.  Virology.  188:408-413. 


18 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


28.  World  Health  Organization.  The  Global  Burden  of  Disease:  2004  update.  2008. 
Ref  Type:  Report 

29.  World  Health  Organization.  Fact  sheet  Number  211.  Influenza  (Seasonal).  Fact 
sheet  Number  21 1. Influenza  (Seasonal) .  2009. 

Ref  Type:  Report 

30.  Xu,  R.,  D.  C.  Ekiert,  J.  C.  Krause,  R.  Hai,  J.  E.  Crowe,  Jr.,  and  I.  A.  Wilson. 

2010.  Structural  basis  of  preexisting  immunity  to  the  2009  H1N1  pandemic 
influenza  virus.  Science.  328:357-360. 


19 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


Figure  Legends 


Figure  1 

Virus  titer  in  lungs  of  C57BL/6  (•)  and  DBA/2  (V  )  mice  2  and  7  days  post  inoculation 
with  with  104  EID50  of  A/shorebird/Delaware/22/2006  (H7N3),  A/blue  winged- 
teal/Alberta/27 1/2007  (H10N5)  and  A/Califomia/4/2009  (HINlpdm).  *  P  <  0.01.  1 
Infectious  titer  in  HIDso/ml  for  H7N3  and  H10N5,  and  in  TCIDso/ml  for  HINlpdm 

Figure  2 

Human  cross-reactive  2009  pandemic  H1N1  neutralizing  antibodies  are  functional  in 
vivo.  Human  serum  pools  with  (A)  or  without  (O)  a  detectable  cross-reactive  2009 
pandemic  H1N1  (A/California/4/2009)  neutralizing  antibody  titer  were  injected 
intraperitoneally  24  h  prior  to  intranasal  lethal  challenge.  Control  mice  were  injected  with 
PBS  (■  )  or  convalescent  serum  obtained  from  2009  pandemic  H1N1  virus-infected 
ferrets  «».  Survival  was  monitored  for  21  days  (A)  and  weight  loss  for  16  days  (B).  Data 
shown  in  (A)  are  the  cumulative  results  of  2  experiments  and  those  in  (B)  represent  the 
average  weight  loss  of  1  indicative  experiment.  '  '  P  <  0.0 1 . 
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Tables 

Table  1  Percent  mortality  in  DBA/2  and  C57BL/6  mice  inoculated  with  various  isolates 
of  influenza  A  virus 


Influenza  A  virus  isolate 

Subtype  - 

DBA/2 

C57BL/6 

106  EIDso  10 

4  EIDso 

106  EIDso 

A/Puerto  Rico/8/ 1934 

H1N1 

100 

100 

100 

A/Memphis/3/2008  1 

H1N1 

100 

88 

0 

A/California/4/2009  2 

H1N1 

100 

100 

60 

A/mallard/ Albert  a/79/2003 

H2N3 

33 

0 

ND4 

A/mallard/ Albert  a/3  3/2004 

H2N4 

0 

0 

0 

X31  (A/Hong  Kong/1/1968) 3 

H3N2 

100 

100 

0 

A/pintail  duck/ Alberta/66/2005 

H4N1 

20 

0 

0 

A/mallard/ Albert  a/147/2007 

H4N6 

0 

ND 

ND 

A/Hong  Kong/213/2003  1 

H5N1 

100 

100 

76 

A/Vietnam/ 1203/2004  1 

H5N1 

100 

100 

100 

AJ  shorebird/Delaware/ 1 01/2004 

H5N7 

20 

0 

0 

A/ruddy  turnstone/Delaware/1 03/2007 

H5N9 

100 

25 

ND 

A/teal/Hong  Kong/W3 12/1997 

H6N1 

100 

100 

40 

A/mallard/ Albert  a/ 1 54/2003 

H6N5 

100 

40 

0 

A/  shorebird/Delaware/22/2006 

H7N3 

100 

100 

0 

A/Netherlands/33/2003  1 

H7N7 

100 

100 

100 

A/mallard/ Albert  a/177/2004 

H7N9 

0 

0 

0 

A/quail/Hong  Kong/Gl/1997 

H9N2 

100 

0 

0 

A/duck/Hong  Kong/Y280/1997 

H9N2 

100 

20 

0 

A/mallard/ Albert  a/162/2007 

H9N5 

0 

0 

ND 

A/mallard/ Albert  a/22 1/2006 

H9N6 

0 

0 

ND 

A/blue-winged  teal/ Alberta/27 1/2007 

H10N5 

100 

75 

0 

A/mallard/ Albert  a/56/2004 

H10N7 

100 

75 

0 

Human  influenza  A  virus  isolate. 

2  Human  2009  pandemic  H1N1  virus  isolate. 

3  X3 1  is  a  6  +  2  reassortant  virus  containing  the  HA  and  NA  gene  segments  from  A/Hong 
Kong/1/68  and  6  gene  segments  from  A/Puerto  Rico/8/1934  virus. 

4  ND,  not  done. 
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Table  2  Serum  antibody  responses  in  mice  infected  with  106  EID50  of  non-  or  low- 
pathogenic  influenza  A  viruses 


Influenza  A  virus  isolate 

Subtype 

Hemagglutination 

Inhibition  titer  1 

DBA/2  C57BL/6 

A/Memphis/3/2008 

H1N1 

ND2 

40 

A/mallard/ Alberta/33/2004 

H2N4 

40 

20 

A/pintail  duck/ Alberta/66/2005 

H4N1 

ND 

<10 

A/mallard/ Alberta/ 1 47/2007 

H4N6 

<10 

ND 

A/shorebird/Delaware/ 1 01/2004 

H5N7 

40 

<10 

A/mallard/ Alberta/ 1 54/2003 

H6N5 

ND 

<10 

A/mallard/ Alberta/ 177/2004 

H7N9 

40 

80 

A/mallard/ Alberta/22 1/2006 

H9N6 

<10 

ND 

A/blue-winged  teal/ Alberta/27 1/2007 

H10N5 

ND 

160 

A/mallard/ Alberta/5 6/2004 

H10N7 

ND 

40 

1  Average  geometric  mean  titer  as  measured  by  the  HI  assay. 

2  ND,  not  done.  HI  titer  <  10  is  undetectable. 


22 


Downloaded  from  jvi.asm.org  at  UT  Health  Sciences  Library  on  June  25,  2010 


Table  3.  Concentrations  of  proinflammatory  cytokines  in  lungs  of  C57BL/6  and  DBA/2 
mice  infected  with  104  EID50  of  pandemic  H1N1  2009  virus 


Cytokine 

Day  2 

Day  7 

(pg/ml) 

DBA/2 

C57BL/6 

DBA/2 

C57BL/6 

CCL2 

4013 

3387 

10145  ’ 

4520 

CCL5 

4212 

2924 

5703 

7060 

IL-6 

710 

950 

1057  * 

735 

TNFa 

236  * 

102 

96 

109 

IFNy 

46 

42 

1097 

900 

Significant  difference  (P<0.05)  in  cytokine  concentration  between  C57BL/6  and  DBA/2 
mice  on  day  2  or  7  after  infection. 
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Coronaviruses  encode  multifunctional  proteins  that  are  critical  for  viral  replication  and  for  blocking  the  innate 
immune  response  to  viral  infection.  One  such  multifunctional  domain  is  the  coronavirus  papain-like  protease 
(PLP),  which  processes  the  viral  replicase  polyprotein,  has  deubiquitinating  (DUB)  activity,  and  antagonizes  the 
induction  of  type  f  interferon  (IFN).  Here  we  characterized  the  DUB  and  IFN  antagonism  activities  of  the  PUP 
domains  of  human  coronavirus  NU63  and  severe  acute  respiratory  syndrome  (SARS)  coronavirus  to  determine  if 
DUB  activity  mediates  interferon  antagonism.  We  found  that  NU63  PUP2  deconjugated  ubiquitin  (Ub)  and  the 
Ub-line  molecule  ISG15  from  cellular  substrates  and  processed  both  lysine-48-  and  lysine-63-  linked  polyubiquitin 
chains.  This  PUP2  DUB  activity  was  dependent  on  an  intact  catalytic  cysteine  residue.  We  demonstrated  that  in 
contrast  to  PUP2  DUB  activity,  PDP2-mediated  interferon  antagonism  did  not  require  enzymatic  activity.  Further¬ 
more,  addition  of  an  inhibitor  that  blocks  coronavirus  protease/DUB  activity  did  not  abrogate  interferon  antago¬ 
nism.  These  results  indicated  that  a  component  of  coronavirus  PUP-mediated  interferon  antagonism  was  indepen¬ 
dent  of  protease  and  DUB  activity.  Overall,  these  results  demonstrate  the  multifunctional  nature  of  the  coronavirus 
PUP  domain  as  a  viral  protease,  DUB,  and  IFN  antagonist  and  suggest  that  these  independent  activities  may 
provide  multiple  targets  for  antiviral  therapies. 


The  front-line  defense  of  a  host  cell  against  virus  infection  is 
the  innate  immune  system,  which  utilizes  multiple  membrane 
and  cytoplasmic  sensors,  such  as  toll-like  receptors  (TLRs)  and 
RNA  helicases,  to  detect  pathogen-associated  molecular  pat¬ 
terns  like  viral  RNA  (3,  9,  31,  47,  54).  Activation  of  these 
sensors  by  viral  RNA  intermediates  sets  off  a  cascade  of  sig¬ 
naling  events  that  ultimately  turn  on  transcription  factors,  such 
as  NF-kB,  ATF2/c-Jun,  IRF-7,  and  IRF-3.  These  activated 
transcription  factors  translocate  to  the  nucleus  and  upregulate 
transcription  of  interferon  (IFN)  mRNAs.  The  translation  and 
subsequent  secretion  of  IFNs  activates  cells  to  upregulate  in¬ 
terferon-stimulated  genes  (ISGs)  to  establish  an  antiviral  state 
hostile  to  viral  replication.  Of  importance  for  this  study,  many 
of  the  signaling  events  that  link  the  sensors  to  the  transcription 
factors  are  mediated  by  the  activities  of  kinases  and  ubiquiti- 
nating  enzymes  that  modify  and  activate  critical  intermediates 
in  the  cascade  (7,  8,  20,  25).  For  example,  signaling  from  RIG-I 
can  proceed  through  MAVS/TRAF3/TANK  to  TBKI  and  in¬ 
ducible  IkB  kinase  (IKKi),  which  ultimately  phosphorylate 
IRF-3.  Recent  studies  indicate  that  both  TRAF3  and  TANK 
are  modified  by  lysine-63-linked  polyubiquitination  and  can  be 
inactivated  by  DUBA,  a  cellular  deubiquitinating  (DUB)  en- 
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zyme  (28).  Thus,  ubiquitinating  enzymes  and  DUBs  are  critical 
players  in  modulating  the  innate  immune  response. 

For  positive-strand  RNA  viruses  that  replicate  in  the  cyto¬ 
plasm  of  the  cell  using  double-stranded  RNA  intermediates, 
the  cytoplasmic  innate  immune  sensors  and  subsequent  signal¬ 
ing  cascades  represent  a  minefield  that  must  be  either  neutral¬ 
ized,  navigated  by  stealth,  or  both.  Recent  studies  have  re¬ 
vealed  that  viral  proteases  can  act  as  “multitaskers”  during 
viral  replication  by  not  only  processing  viral  polyproteins  but 
also  cleaving/inactivating  key  players  in  the  innate  immune 
response.  For  example,  the  hepatitis  C  virus  NS3-4A  protease 
cleaves  the  viral  replicase  polyprotein  and  inactivates  the  key 
signaling  proteins  TRIF  and  MAVS  (36,  38,  42,  46).  Picorna- 
virus  3C  protease  is  essential  for  processing  the  replicase 
polyprotein  and  inactivating  NF-kB  and  RIG-I  (2,  49).  In  these 
studies,  the  catalytic  function  of  the  viral  proteases  was  essen¬ 
tial  for  the  inactivation  of  host  factors  involved  in  signaling  the 
innate  immune  response.  Recently  the  coronavirus  papain-like 
protease  domains  have  also  been  identified  as  modulators  of 
the  innate  immune  response;  however,  the  mechanisms  of  in¬ 
hibition  are  not  entirely  clear  (13,  16,  70). 

Coronaviruses  are  enveloped  viruses  with  large  RNA  ge¬ 
nomes  (28  to  32  kb)  that  cause  disease  in  humans  ranging  from 
common  colds  (human  coronavirus  [CoV]  229E  [FICoV-229E] 
and  OC43)  to  croup  and  pneumonia,  seen  mostly  in  very  young 
and  old  populations  (HCoV-NL63  and  -HKU1),  to  severe 
acute  respiratory  syndrome  (SARS)  coronavirus  (SARS-CoV) 
with  10%  mortality  (53).  Upon  entry,  coronavirus  genomic 
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RNA  is  translated  to  produce  two  large  polyproteins,  ppla  and 
pplab.  These  polyproteins  are  processed  by  viral  cysteine 
proteases,  papain-like  (PLPs/PLpro)  and  picornavirus  SC- 
like  (3CLpro),  to  generate  mature  nonstructural  proteins 
(nsp’s)  that  assemble  with  host  cell  membranes  to  form  double 
membrane  vesicles  (DMVs)  (18,  19,  61).  These  DMVs  are  the 
site  of  viral  RNA  synthesis  producing  double-stranded  inter¬ 
mediates  and  genomic/subgenomic  mRNAs.  Interestingly,  ro¬ 
bust  replication  of  SARS-CoV  was  shown  to  trigger  low  but 
detectable  levels  of  beta  interferon  (IFN-(3)  (13,  62,  63).  The 
low-level  IFN  response  to  a  vigorously  replicating  RNA  virus 
suggests  that  SARS-CoV  either  evades  or  inactivates  the  in¬ 
nate  immune  response.  Indeed,  recent  studies  indicate  that 
SARS-CoV  encodes  multiple  proteins,  such  as  nspl,  protein 
3b,  protein  6,  and  the  nucleocapsid  protein  that  modulate 
multiple  pathways  of  the  innate  immune  response  (17,  27,  30, 
48,  68).  In  addition,  we  showed  that  the  SARS-CoV  papain¬ 
like  protease  (PLpro)  domain  acts  as  an  interferon  antagonist 
that  blocks  the  phosphorylation  and  subsequent  nuclear  trans¬ 
location  of  IRF-3  (13).  We  also  showed  via  X-ray  structural 
studies  that  the  SARS-CoV  PLpro  domain  is  similar  to  cellular 
deubiquitinating  enzymes  (57),  and  we  and  others  demon¬ 
strated  that  PLpro  is  both  a  protease  and  a  DUB  (4,  39,  40). 
Initially,  we  hypothesized  that  either  the  protease  or  DUB 
activity  would  be  required  for  modulating  the  innate  immune 
response,  but  we  found  that  inactivation  of  the  catalytic  cys¬ 
teine  residue  of  PLpro,  which  ablates  both  proteolysis  and 
deubiquitinating  activity,  decreased  but  did  not  abrogate 
PLpro’s  ability  to  block  activation  of  interferon  (13).  These  results 
are  consistent  with  a  protease/DUB-independent  mechanism 
that  contributes  to  interferon  antagonism.  Frieman  and  co¬ 
workers  also  showed  that  catalytic  mutants  of  PLpro  retained 
interferon  antagonism,  and  they  reported  that  deletion  of  the 
PLP  ubiquitin-like  (Ubl)  domain  upstream  of  the  catalytic  site 
resulted  in  a  loss  of  antagonism  (16).  Studies  by  Zheng  et  al.  of 
the  PLP  domain  of  murine  hepatitis  virus  (MHV)  suggested 
that  viral  DUB  activity  may  be  required  for  interferon  antag¬ 
onism,  although  they  reported  that  MHV  PLP2  catalytic  cys¬ 
teine  mutants  became  less  efficient  at,  but  did  not  eliminate, 
blockade  of  the  interferon  response  (70).  Therefore,  further 
studies  are  required  to  clarify  the  role  of  coronavirus  protease/ 
DUB  activity  in  PLP-mediated  interferon  antagonism. 

In  this  study,  we  analyzed  the  DUB  and  IFN  antagonism 
profiles  of  the  papain-like  proteases  of  human  coronavirus 
NL63  and  SARS-CoV  (group  1  and  group  2  coronaviruses, 
respectively).  We  show  that  HCoV-NL63  PLP2  is  a  deubiq¬ 
uitinating  and  delSGylating  (ISG15-removing)  enzyme. 
HCoV-NL63  infection,  like  that  of  SARS-CoV,  triggers  a  weak 
type  I  IFN  response  in  human  airway  epithelial  cell  cultures. 
We  also  evaluated  the  role  of  PLP2  and  PLpro  enzymatic 
activity  in  interferon  antagonism.  By  generating  dose-response 
profiles  of  IFN  antagonism,  we  found  that  the  papain-like 
proteases  do  not  require  enzymatic  activity  to  inhibit  type  I 
IFN  induction.  Using  a  pharmacological  approach,  we  found 
that  the  inhibition  of  PLpro  did  not  alter  the  antagonism  of 
IRF-3-dependent  reporters  but  did  affect  an  NF-kB-  depen¬ 
dent  reporter.  Overall,  we  show  that  multifunctional  corona¬ 
virus  PLPs  target  the  activities  of  key  transcription  factors 
involved  in  the  induction  of  type  I  interferons  and  thereby 
hinder  the  activation  of  the  innate  immune  system. 


MATERIALS  AND  METHODS 

Cells  and  HCoV-NL63.  HeLa  cells,  HEK293  cells,  and  HEK293-TLR3  (stable 
expression  of  human  TLR3  receptor)  cells  were  cultured  using  Dulbecco’s  mod¬ 
ified  Eagle’s  medium  containing  10%  (vol/vol)  fetal  calf  serum,  supplemented 
with  penicillin  (100  U/ml)  and  streptomycin  (100  pg/ml).  The  HCoV-NL63  (P8) 
virus  and  LLC-MK2  cells  were  kindly  provided  by  Lia  van  der  Hoek  (University 
of  Amsterdam,  Amsterdam,  Netherlands)  and  propagated  as  described  previ¬ 
ously  (12).  A  plaque-purified  isolate  of  HCoV-NL63  was  kindly  provided  by 
Christian  Drosten  and  propagated  in  CaCo2-TC7  cells  (21).  This  virus  stock  was 
used  to  infect  human  airway  epithelial  cells  as  described  previously  (1). 

Plasmid  DNA.  pcDNA3.1-Flag-Ub  was  kindly  provided  by  Adriano  Marchese 
(Loyola  University  Medical  Center).  IFN-(3-Luc  was  kindly  provided  by  John 
Hiscott  (Jewish  General  Hospital,  Montreal,  Canada).  pISRE-Luc  has  been 
previously  described  (35).  pRL-TK  Renilla  luciferase  reporter  was  purchased 
from  Promega.  N-RIG-Flag,  NF-kB-Luc,  and  nsp2/3-GFP  were  kindly  provided 
by  Ralph  Baric  (University  of  North  Carolina).  pcDNA3-myc6-mISG15  was 
kindly  provided  by  Min-Jung  Kim  (Pohang  University  of  Science  and  Technol¬ 
ogy,  Pohang,  Republic  of  Korea).  pcDNA3-UbelL  and  pcDNA3-UbcH8  were 
kind  gifts  from  Robert  M.  Krug  (University  of  Texas). 

PLP1  and  PLP2  core  domain  synthesis,  cloning,  and  site-directed  mutagen¬ 
esis.  To  obtain  high  expression  in  eukaryote  cells,  the  codon  usage  of  the 
HCoV-NL63  PLP1  core  domain  (amino  acids  1018  to  1277  of  HCoV-NL63)  and 
PLP2  core  domain  (amino  acids  1570  to  1884)  were  optimized  based  on  human 
codon  usage  frequency,  and  the  potential  splicing  sites  and  poly(A)  signal  se¬ 
quences  were  removed  and  cloned  into  pcDNA3.1-V5/HisB  at  the  BamHI  and 
EcoRI  sites  as  an  in-frame  fusion  with  the  V5  peptide.  The  native  viral  sequence 
for  the  remainder  of  nsp3  (including  the  transmembrane  domain  downstream  of 
PLP2)  was  cloned  into  pcDNA3.1-PLP2(N)  using  the  existing  EcoRI  site  and 
Xhol  to  generate  transmembrane  (TM)-containing  PLP2  (PLP2-TM)  in  frame 
with  the  V5  peptide.  To  generate  specific  mutations  in  the  catalytic  residues 
(C1062  and  H1212  in  PLP1  and  C1678  and  H1836  in  PLP2),  mutagenic  primers 
were  incorporated  into  newly  synthesized  DNA  using  the  QuikChange  II  XL 
site-directed  mutagenesis  protocol  (Stratagene,  La  Jolla,  CA)  according  to  the 
manufacturer’s  instruction.  Mutated  nucleotides  are  indicated  in  bold.  PLP1 
C1062A  was  amplified  using  the  forward  primer  5'  AAC  AAC  GCC  TGG  ATC 
AGC  ACC  ACC  CTG  GTG  CAA  CTG  3'  and  reverse  primer  5'  GAT  CCA 
GGC  GTT  GTT  GTC  GCT  CTG  GTC  CAG  CAC  CCG  3'.  PLP1  H1212A  was 
amplified  using  the  forward  primer  5'  AGC  GGC  GCC  TAC  CAG  ACC  AAC 
CTG  TAC  AGC  TTC  AAC  3'  and  reverse  primer  5'  CTG  GTA  GGC  GCC 
GCT  GCC  CTT  CAC  GCC  CAG  GTA  CAC  3'.  PLP2  C1678A  was  amplified 
using  the  forward  primer  5'  AAC  AAC  GCC  TGG  GTG  AAC  GCC  ACC  TGC 
ATC  ATC  CTG  3'  and  reverse  primer  5'  CAC  CCA  GGC  GTT  GTT  GTC  GGT 
GGT  GCC  CAG  CAC  CCG  3'.  PLP2  H1836A  was  amplified  using  the  forward 
primer  5'  AAC  GGC  GCC  TAC  GTG  GTG  TAC  GAC  GCC  GCC  AAC  AAC 
3'  and  reverse  primer  5'  CAC  GTA  GGC  GCC  GTT  GTC  GAA  GCT  GCC 
GCT  GAA  GG  3'.  The  primers  used  for  mutagenesis  were  designed  according 
to  the  modified  methods  of  Zheng  et  al.  (71).  All  introduced  mutations  were 
confirmed  by  DNA  sequencing. 

Assay  of  deubiquitinating  activity  in  cultured  cells.  The  effect  of  HCoV-NL63 
PLP1  and  PLP2  on  ubiquitinated  proteins  in  cultured  cells  was  assessed  as 
described  previously  (14).  Briefly,  HeLa  cells  cultured  in  60-mm  dishes  were 
cotransfected  with  0.4  pg  of  pcDNA3.1-Flag-Ub  plus  appropriate  amounts  of 
constructs  containing  PLP1,  PLP2,  or  the  corresponding  catalytic  mutants. 
Transfection  was  performed  with  Lipofectamine  2000  per  the  manufacturer’s 
instructions.  The  empty  vector  pcDNA3.1/V5-HisB  was  used  to  standardize  the 
total  amount  of  DNA  used  for  transfection.  After  48  h,  cells  were  harvested  by 
adding  250  pi  2X  LB  A,  containing  20  mM  A-ethylmaleimide  (NEM)  (Calbio- 
chem)  and  20  mM  iodoacetamine  (Sigma).  Cell  lysates  were  then  analyzed  for 
ubiquitin  (Ub)-conjugated  proteins  by  Western  blotting  with  anti-Flag  M2  anti¬ 
body  (1:10,000)  (Sigma).  Proteins  were  separated  via  SDS-PAGE,  followed  by 
transfer  to  a  polyvinylidene  difluoride  (PVDF)  membrane  in  transfer  buffer 
(0.025  M  Tris,  0.192  M  glycine,  and  20%  methanol)  for  2  h  at  4°C.  The  mem¬ 
brane  was  blocked  using  5%  dried  skim  milk  in  Tris-buffered  saline  (TBS)  (0.9% 
NaCl,  10  mM  Tris-HCl,  pH  7.5)  plus  0.1%  Tween  20  (TBST)  for  2  h  at  room 
temperature.  The  blot  was  probed  with  the  indicated  antibody  overnight  at  4°C. 
The  membrane  was  washed  in  TBST  three  times  for  20  min  (each).  Following  the 
washes,  the  membrane  was  incubated  with  peroxidase-conjugated  secondary 
antibody  (donkey  anti-rabbit  or  goat  anti-mouse  IgG  horseradish  peroxidase 
[HRP])  (Amersham)  at  a  dilution  of  1:10,000  for  2  h  at  room  temperature.  The 
membrane  was  then  washed  three  times  with  TBST  and  detected  with  Western 
Lightning  Chemoluminescence  Reagent  Plus  (PerkinElmer  LAS  Inc.).  To  con¬ 
firm  the  expression  levels  of  PLP1,  PLP2,  and  the  mutants,  anti-V5  antibody 
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(Invitrogen)  was  used  to  detect  the  V5-tagged  proteins.  Calnexin  was  detected 
with  anti-calnexin  monoclonal  antibody  (MAb)  (BD  Transduction  Lab)  as  a 
protein  loading  control. 

NL63  PLP2  cleavage  of  K48-  and  K63-linked  ubiquitin  chains.  The  NL63 
PLP2  wild-type  protein  was  purified  as  previously  described  (4, 12),  and  ubiquitin 
chains  were  purchased  from  Boston  Biochem  (K48-Ub6  [catalog  no.  UC-217] 
and  K63-Ub6  [catalog  no.  UC-317]).  Proteolytic  cleavage  of  K48-linked  or  K63- 
linked  ubiquitin  chains  was  carried  out  under  the  following  conditions:  0.01  p,g  of 
purified  PLP2  was  incubated  with  2.5  to  5  pg  of  K48-Ub6  or  K63-Ub6  at  25°C  in 
a  10-  to  20-pl  volume  containing  50  mM  HEPES,  pH  7.5,  0.1  mg/ml  bovine 
serum  albumin  (BSA),  100  mM  NaCl,  and  2  mM  dithiothreitol  (DTT).  A  control 
reaction  was  incubated  under  identical  conditions  with  the  exclusion  of  enzyme. 
At  specified  time  points,  the  reactions  were  quenched  with  the  addition  of 
SDS-PAGE  sample  loading  dye  to  a  1 X  concentration  (25  mM  Tris,  pH  6.8,  280 
mM  (3-mercaptoethanol,  4%  glycerol,  0.8%  SDS,  0.02%  bromophenol  blue)  and 
heat  treated  at  95°C  for  5  min.  The  samples  were  analyzed  by  electrophoresis  on 
a  15%  SDS-PAGE  gel  and  stained  with  Coomassie  dye. 

Assay  of  delSGylating  activity  in  cultured  cells.  The  effect  of  HCoV-NL63 
PLP2  on  ISGylated  proteins  in  cultured  cells  was  assessed  as  described  previ¬ 
ously  (69).  Briefly,  HEK293  cells  cultured  in  60-mm  dishes  were  cotransfected 
with  0.5  pg  of  pcDNA3-myc6-mISG15,  0.25  pg  of  pcDNA3-UbelL,  and  0.25  pg 
of  pcDNA3-UbcH8  plus  appropriate  amounts  of  constructs  containing  PLP2  or 
the  corresponding  catalytic  mutant  with  a  total  of  2  pg  of  plasmid  DNA  for  each 
transfection.  The  empty  vector  pcDNA3.1/V5-HisB  was  used  to  standardize  the 
total  amount  of  DNA  used  for  transfection.  Transfection  was  performed  with 
Lipofectamine  2000  per  the  manufacturer’s  instructions.  After  30  h,  cells  were 
harvested  by  adding  250  pi  2X  LB  A  containing  20  mM  /V-ethylmaleimide  (Cal- 
biochem)  and  20  mM  iodoacetamine  (Sigma).  Cell  lysates  were  then  analyzed  for 
ISG-conjugated  proteins  by  Western  blotting  with  monoclonal  anti-myc  antibody 
(1:2,000;  MBL  Companies,  Japan)  as  described  above.  To  confirm  the  expression 
levels  of  PLP2  and  the  mutants,  anti-V5  antibody  (1:5,000;  Invitrogen)  was  used 
to  detect  the  V5-tagged  proteins.  Actin  was  detected  with  antiactin  MAb  (Beyo- 
time  Institute  of  Biotechnology,  China)  as  a  protein  loading  control. 

Enzyme-linked  immunosorbent  assay  for  IFN-|J  secretion  in  HAE  culture 
supernatants.  Human  airway  epithelial  (HAE)  cultures  were  generated  as  pre¬ 
viously  described  (1).  Cultures  were  inoculated  with  100  pi  of  105  PFU/ml 
HCoV-NL63,  2,000  hemagglutinating  units  (HAU)/ml  Sendai  virus  (SeV),  or 
medium  for  1  h  at  37°C,  after  which  the  inoculum  was  removed  and  apical  washes 
with  200  pi  of  F12  medium  (Gibco)  were  performed  at  indicated  times.  The 
IFN-(3  concentration  was  determined  by  a  commercial  enzyme-linked  immu¬ 
nosorbent  assay  (ELISA)  (PBL  Biomedical  Laboratories)  according  to  the  man¬ 
ufacturer’s  instructions.  Data  were  analyzed  and  plotted  using  the  GraphPad 
Prism  5.0  software  program. 

Luciferase  reporter  gene  assay.  To  determine  if  PLpro  and/or  PLP2  modulates 
IFN  expression  in  host  cells,  the  IFN-(3-Luc  reporter,  consisting  of  the  IFN-p 
promoter  upstream  of  firefly  luciferase,  was  transfected  into  HEK293  cells  along 
with  PLpro,  AUbl-PLpro-Sol,  AUbl-PLpro-TM,  PLP2,  PLP2-TM,  or  the  C1678A 
and  H1836A  PLP2  catalytic  mutants  in  the  soluble  or  transmembrane  versions. 
pRL-TK,  encoding  Renilla  luciferase  under  the  control  of  the  herpes  simplex 
thymidine  kinase  promoter  (low  to  moderate  Renilla  expression),  was  used  for 
normalization  of  transfection  efficiency.  HEK293  cells  were  transfected  by  Li¬ 
pofectamine  2000  or  LT1  transfection  reagent  (MirusBio)  according  to  the 
manufacturer’s  instructions  and  incubated  for  24  h.  Cells  were  then  mock  in¬ 
fected  or  infected  with  Sendai  virus  (Cantell  strain;  Charles  River  Laboratories) 
at  the  dose  of  100  HAU/ml  for  16  h  or  transfected  with  N-RIG-Flag  for  24  h  to 
activate  the  RIG-I-dependent  IFN  pathway.  To  detect  TLR3-dependent  IFN 
expression,  HEK293-TLR3  cells  were  transfected  with  IFN-(3-Luc  and  PLP2  for 
24  h.  Cells  were  then  treated  with  poly(IC)  for  6  h  as  described  previously  (13). 
Firefly  luciferase  and  Renilla  luciferase  activities  were  assayed  using  the  Dual- 
Luciferase  reporter  assay  kit  (Promega)  and  a  luminometer. 

Drug  inhibition  of  SARS-CoV  PLpro.  HEK293  cells  were  transfected  with 
plasmids  encoding  PLpro-TM  (13),  IFN-p-Luc,  ISRE-Luc,  NF-kB-Luc,  pRL- 
TK,  nsp2/3-GFP  substrate,  and/or  N-RIG-Flag.  The  pcDNA3.1/V5-HisB  vector 
was  used  to  standardize  the  DNA  concentration  for  transfection.  Dimethyl 
sulfoxide  (DMSO)  vehicle  control  or  100  (jlM  GRL-0617S  (56)  was  added  at  the 
time  of  transfection.  Tumor  necrosis  factor  alpha  (TNF-a)  (10  ng/ml)  from 
Roche  was  used  to  stimulate  the  NF-kB-Luc  reporter  (6  h  of  stimulation).  Cell 
lysates  were  prepared  per  the  manufacturer’s  instruction  using  IX  passive  lysis 
buffer  (Promega),  and  luciferase  activity  was  measured  using  the  Dual-Lucif- 
erase  reporter  assay  kit  (Promega)  and  a  luminometer.  Cell  lysates  were  also 
mixed  1:1  with  2x  sample  buffer  and  subjected  to  Western  blotting  as  described 
above.  PLpro-TM  was  detected  with  anti-V5  (Invitrogen),  and  nsp2/3-GFP  and 


nsp3-GFP  were  detected  with  rabbit  anti-green  fluorescent  protein  (anti-GFP) 
(Invitrogen). 

RESULTS 

NL63  PLP2  but  not  PLP1  has  deubiquitinating  activity. 

Previously  we  showed  that  HCoV-NL63  replicase  nonstruc- 
tural  protein  3  (nsp3)  encodes  two  papain-like  proteases,  PLP1 
and  PLP2,  that  process  the  amino-terminal  end  of  the  replicase 
polyprotein  (Fig.  1A).  In  addition,  we  have  shown  that  PLP2 
can  process  K48-linked  polyubiquitin  chains,  suggesting  that 
this  protease  has  deubiquitinating  activity  (12).  Both  the 
polyprotein  cleavage  and  K48-linked  polyubiquitin  chain  pro¬ 
cessing  are  dependent  on  a  cysteine  residue  in  the  catalytic 
triad  of  this  cysteine  protease  (12).  To  determine  if  one  or  both 
of  these  NL63  PLPs  can  remove  ubiquitin  conjugated  to  cel¬ 
lular  proteins,  HeLa  cells  were  transfected  with  plasmid  DNA 
encoding  PLP1  or  PLP2  along  with  pcDNA3.1-3XFlag-Ub, 
and  the  effect  of  expression  of  PLP1  and  PLP2  on  the  extent  of 
ubiquitinated  cellular  proteins  was  assessed  via  Western  blot¬ 
ting  with  anti-Flag  antibodies.  We  found  that  increasing  ex¬ 
pression  of  PLP2  resulted  in  a  dramatic  reduction  in  the  level 
of  Ub-conjugated  proteins  (Fig.  IB).  We  noted  that  PLP2 
appears  to  have  global  DUB  activity,  since  no  particular  cel¬ 
lular  substrates  were  spared.  As  expected,  this  PLP2  DUB 
activity  is  dependent  on  an  intact  catalytic  cysteine  residue 
1678,  and  mutation  of  this  residue  to  alanine  resulted  in  the 
loss  of  DUB  activity  (Fig.  1C).  In  contrast,  PLP1  did  not  show 
any  significant  reduction  of  Flag-Ub  conjugates  at  the  concen¬ 
trations  tested  (Fig.  ID).  These  results  indicate  that  NL63 
PLP2,  like  SARS-CoV  PLpro,  has  potent  DUB  activity  that 
can  remove  ubiquitin  conjugates  from  many  cellular  substrates. 

NL63  PLP2  processes  lysine-63-linked  in  addition  to  lysine- 
48-linked  polyubiquitin.  Cellular  proteins  can  be  covalently 
modified  with  ubiquitin  at  one  or  multiple  lysines  through  an 
isopeptide  bond  that  links  the  carboxy  terminus  of  ubiquitin  to 
a  lysine  on  the  target  protein.  Ubiquitin  itself  contains  seven 
lysines  that  can  be  further  conjugated  to  the  C  terminus  of 
another  ubiquitin  molecule,  forming  different  types  of  poly- 
ubiquitin-linked  chains  on  the  targeted  protein.  The  two  most 
common  types  of  polyubiquitin  chains  are  linked  through  ubiq¬ 
uitin  lysine  48  (K48)  and  lysine  63  (K63).  These  modifications 
play  key  regulatory  roles  in  protein  degradation  and  pathway 
signaling  and  have  been  associated  with  controlling  several 
pathways  of  innate  and  adaptive  immunity  (7).  Previous  studies 
indicated  that  SARS-CoV  PLpro  processes  both  K48-  and 
K63-linked  ubiquitin  (39,  40).  To  assess  if  HCoV-NL63  PLP2 
has  isopeptidase  activity  that  will  deconjugate  K63-linked  ubiq¬ 
uitin  in  addition  to  K48-linked  isopeptidase  activity,  purified 
PLP2  enzyme  was  incubated  with  hexameric  K48-linked  and 
K63-linked  polyubiquitin  chains.  SDS-PAGE  analysis  of  the 
cleavage  products  shows  that  PLP2  can  cleave  the  substrates  in 
vitro  (Fig.  2),  and  with  prolonged  incubation  times,  both  chains 
can  be  reduced  to  monoubiquitin  (data  not  shown).  These  data 
show  that  both  major  forms  of  polyubiquitin  can  be  recognized 
and  degraded  by  HCoV-NL63  PLP2. 

PLP2  possesses  delSGylating  activity.  Several  viral  DUBs, 
including  SARS-CoV  PLpro,  can  also  deconjugate  Ub-like  moi¬ 
eties  such  as  ISG15  (39,  40).  Conjugation  of  ISG15  has  been 
shown  to  be  important  for  protection  against  viral  infection 
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FIG.  1.  HCoV-NL63  PLP2,  but  not  PLP1,  has  a  dose-dependent  global  deubiquitinating  activity  in  cultured  cells.  (A)  Schematic  diagram  of 
the  NL63  genomic  RNA  and  the  resulting  polyprotein  lab,  which  contains  three  viral  proteases.  PLP1  and  PLP2  cleavage  sites  are  indicated,  as 
are  the  resulting  nonstructural  proteins.  The  V5-tagged  constructs  of  PLP2  used  in  this  study  are  listed,  and  the  catalytic  residues  numbering  from 
ORF  la  are  shown.  DNA  encoding  HCoV-NL63  PLP2  (B),  PLP2  C1678A  (C),  or  PLP1  (D)  was  transfected  into  HeLa  cells  along  with 
pcDNA3.1-3xFlag-Ub.  Cell  lysates  were  prepared  at  24  h  posttransfection  and  analyzed  for  Flag-Ub-conjugated  proteins  by  Western  blotting  with 
an  anti-Flag  antibody.  Mouse  anti-V5  was  used  to  confirm  the  expression  of  PLP1  and  PLP2,  and  anticalnexin  antibody  was  used  to  detect  calnexin, 
which  serves  as  a  protein  loading  control.  Molecular  weight  markers  shown  on  the  left  of  each  gel  are  in  thousands. 


(15,  33,  34).  We  assessed  whether  HCoV-NL63  PLP2  can 
delSGylate  cellular  c-myc-tagged  ISG15  (c-myc-ISG15)  conju¬ 
gates.  HeLa  cells  were  transfected  with  increasing  amounts  of 
plasmid  DNA  encoding  PLP2  along  with  c-myc-ISG15  and  the 
ISG15  conjugation  machinery  UbelL  and  UbcH8  to  enhance 
ISGylation  of  host  cell  proteins.  The  ability  of  PLP2  to 
delSGylate  cellular  proteins  was  then  assayed  via  Western  blot¬ 
ting  with  anti-myc  antibody.  We  found  that  expression  of  in¬ 
creasing  amounts  of  PLP2  was  associated  with  a  dramatic  re¬ 
duction  in  the  levels  of  ISGylated  cellular  proteins  (Fig.  3A),  in 
agreement  with  a  previous  report  (50).  The  delSGylating  ac¬ 
tivity  was  dependent  on  intact  catalytic  residues  C1678  and 
H1836,  since  mutation  of  these  residues  to  alanine  resulted 
in  the  loss  of  delSGylating  activity  (Fig.  3B).  Thus,  HCoV- 


NL63  PLP2  is  a  potent  DUB/delSGylating  enzyme  that  acts 
on  many  modified  cellular  substrates. 

IFN-P  release  from  human  airway  epithelial  (HAE)  cell 
cultures.  FIAE  cultures  are  a  useful  model  system  for  studying 
human  respiratory  viruses,  including  HCoV-NL63,  since  they 
mimic  human  bronchial  epithelium,  which  is  the  primary  site  of 
infection  (1,  60).  HAE  cell  cultures  were  infected  with  HCoV- 
NL63  or  Sendai  virus,  and  apical  wash  samples  were  collected 
at  24,  48,  72,  96,  and  144  h  postinfection.  The  presence  of 
secreted  IFN-p  in  the  apical  wash  was  measured  by  ELISA. 
Mock-infected  cultures  released  no  detectable  IFN-p,  whereas 
cultures  inoculated  with  a  potent  IFN  inducer,  Sendai  virus, 
released  more  than  400  pg/ml  of  IFN-p.  In  contrast,  inocula¬ 
tion  of  HAE  cultures  with  HCoV-NL63  stimulated  low  but 
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FIG.  2.  Processing  of  K48-  and  K63-linked  ubiquitin  chains  by 
PLP2.  NL63  PLP2  was  incubated  with  K48-linked  (left)  or  K63- 
linked  (right)  Ub6  chains  for  the  indicated  time  points  before  being 
analyzed  by  SDS-PAGE.  Uncleaved  Ub6  is  run  in  the  first  lane  of 
each  gel.  The  expected  sizes  of  the  Ub  species  are  indicated  to  the 
left  of  all  gels.  Molecular  weight  (MW)  markers  include  250,000-, 
100,000-,  75,000-,  50,000-,  37,000-,  25,000-,  20,000-,  15,000-,  and 
10,000-molecular-weight  bands. 


detectable  levels  of  IFN-p  from  24  to  144  h  postinfection  (Fig. 
4).  These  results  are  reminiscent  of  the  reports  of  low  but 
detectable  levels  of  IFN-p  produced  from  SARS-CoV-infected 
cells  (13,  63)  and  indicate  that  either  HCoV-NL63  is  a  weak 
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FIG.  3.  NL63  PLP2  has  global  delSGylating  activity  in  cultured 
cells.  HEK293  cells  were  transfected  with  pcDNA3-myc6-mISG15, 
pcDNA3-UbelL,  and  pcDNA3-UbcH8  plus  indicated  amounts  of  the 
PLP2  expression  construct  (A)  or  PLP2  expression  construct  and  the 
corresponding  catalytic  mutants  (B).  To  ensure  that  the  total  amount 
(2  p,g/transfection)  of  plasmids  for  transfection  was  equal  under  every 
condition,  empty  vector  pcDNA3.1W5-HisB  (v)  was  used  to  standard¬ 
ize  the  total  amount  of  DNA.  After  30  h,  cells  were  harvested,  and  cell 
lysates  were  analyzed  for  ISG-conjugated  proteins  by  Western  blotting 
with  monoclonal  anti-myc  antibody.  Expression  levels  of  V5-tagged 
PLP2  and  mutant  enzymes  were  detected  with  anti-V5  antibody.  Actin 
was  detected  with  antiactin  MAb  antibody  as  a  protein  loading  control. 


Mock  SeV  MCoV  NL63 

FIG.  4.  Evaluating  interferon-p  secretion  from  human  airway  epi¬ 
thelial  (HAE)  cell  cultures  following  HCoV-NL63  infection.  Apical 
washes  were  collected  from  HAE  cultures  at  24,  48,  72,  96,  and  120  h 
postinfection,  and  secreted  IFN-P  was  measured  by  ELISA.  Superna¬ 
tant  from  HAE  cells  infected  with  Sendai  virus  (SeV)  for  24  h  was  used 
as  a  positive  control. 

inducer  of  the  IFN-p  response  or,  like  SARS-CoV,  HCoV- 
NL63  encodes  potent  IFN  antagonists. 

PLP2  inhibits  both  RIG-I-  and  TLR3-dependent  IFN-P  ex¬ 
pression.  To  determine  if  PLP2  is  an  IFN  antagonist,  we  trans¬ 
fected  HEK293  cells  with  plasmids  encoding  HCoV-NL63 
PLP2  or  SARS-CoV  PLpro  along  with  IFN-p-luciferase  and 
Renilla  luciferase  reporters  for  24  h.  Then,  the  RIG-I-depen- 
dent  pathway  leading  to  IFN-p  expression  was  activated  by 
Sendai  virus  infection  for  16  h  or  by  a  dominant  active  N- 
terminal  portion  of  RIG-I  (N-RIG).  We  found  that  activation 
of  the  IFN-p  promoter  by  Sendai  virus  (Fig.  5A)  or  N-RIG 
(Fig.  5B)  was  inhibited  in  the  presence  of  either  NL63  PLP2  or 
SARS-CoV  PLpro.  To  determine  if  HCoV-NL63  PLP2  inhib¬ 
its  TLR3-mediated  activation  of  IFN-p  production,  PLP2  and 
the  reporters  were  transfected  into  HEK293-TLR3  cells,  and 
then  the  TLR3-mediated  pathway  was  activated  by  addition  of 
poly(IC)  to  the  cell  culture  medium.  We  found  that  activation 
of  the  IFN-p  promoter  by  the  TLR3-mediated  pathway  was 
significantly  inhibited  by  HCoV-NL63  PLP2  and  SARS-CoV 
PLpro  (Fig.  5C).  These  results  demonstrate  that  the  IFN  an¬ 
tagonism  mediated  by  coronavirus  PLPs  is  conserved  in  two 
distinct  viruses,  although  there  is  only  19%  amino  acid  identity 
between  the  catalytic  domains  of  HCoV-NL63  PLP2  and 
SARS-CoV  PLpro  in  this  region  of  nsp3  (4). 

Mutation  of  the  catalytic  residues  does  not  abolish  HCoV- 
NL63  PLP2  IFN  antagonism.  To  further  elucidate  the  inter¬ 
feron  antagonism  profile  of  PLP2,  a  dose-dependent  IFN  an¬ 
tagonism  profile  was  generated.  Concurrently,  to  determine  if 
IFN  antagonism  is  dependent  on  catalytic  activity,  plasmid 
DNAs  encoding  PLP2  or  the  C1678A  (devoid  of  enzymatic 
activity)  and  H1836A  (severely  reduced)  catalytic  mutants 
were  transfected  with  the  IFN-p  and  pRL-TK  reporters  into 
HEK293  cells,  and  IFN-p  reporter  activity  was  assessed.  N- 
RIG  was  used  to  stimulate  IFN-p  induction.  We  found  that, 
like  the  PLP2  wild  type,  both  PLP2  C1678A  and  PLP2  H1836A 
exhibit  dose-dependent  inhibition  of  IFN-p  promoter  activity; 
however,  the  catalytic  mutants  were  less  efficient  than  wild- 
type  PLP2  (Fig.  6A).  Expression  of  the  proteases  was  verified 
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FIG.  5.  PLP2  inhibits  both  RIG-I-  and  TLR3-dependent  IFN-(3  induction.  (A)  HEK293  cells  were  transfected  with  IFN-pS-Luc,  pRL-TK,  and 
either  300  ng  of  HCoV-NL63  PLP2  or  300  ng  of  SARS-CoV  PLpro.  At  24  h  posttransfection,  cells  were  either  mock  infected  or  infected  with 
Sendai  virus  for  16  h.  Following  infection,  cell  lysates  were  prepared  and  assayed  using  the  Dual-Luciferase  reporter  assay.  (B)  HEK293  cells  were 
transfected  with  IFN-0-Luc,  pRL-TK,  200  ng  N-RIG,  and  either  300  ng  of  HCoV-NL63  PLP2  or  300  ng  of  SARS-CoV  PLpro.  At  24  h 
posttransfection,  cell  lysates  were  prepared  and  assayed  using  the  Dual-Luciferase  reporter  assay.  (C)  FIEK293-TLR3  cells  were  transfected  with 
IFN-(3-Luc,  pRL-TK,  and  either  300  ng  of  HCoV-NL63  PLP2  or  300  ng  of  SARS-CoV  PLpro.  At  24  h  posttransfection,  cells  were  either  mock 
treated  or  treated  with  poly(IC)  for  6  h.  Following  poly(IC)  treatment,  cell  lysates  were  prepared  and  assayed  using  the  Dual-Luciferase  reporter 
assay.  Error  bars  indicate  standard  deviations  from  the  means  for  triplicates. 


by  Western  blotting  (Fig.  6B).  These  results  indicate  that  en¬ 
zymatic  activity  of  PLP2  is  not  strictly  required  for  inhibition  of 
antiviral  IFN  expression. 

The  transmembrane  (TM)  form  of  PLP2  is  also  an  IFN 
antagonist.  The  PLP  domains  are  part  of  a  larger  nsp3  protein 
in  SARS-CoV  and  FICoV-NL63  that  is  membrane  bound.  Pre- 
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FIG.  6.  HCoV-NL63  PLP2  and  the  catalytic  mutants  inhibit  RIG- 
I-mediated  IFN  induction  in  a  dose-dependent  manner.  PLP2  and  the 
C1678A  and  FI1836A  catalytic  mutants  were  cotransfected  with  IFN- 
(3-Luc  and  pRL-TK  reporters  into  FIEK293  cells.  A  dominant  active 
N-terminal  portion  of  RIG-I  was  used  to  stimulate  IFN-pl  induction. 
At  24  h  posttransfection,  cell  lysates  were  harvested  and  assayed  for 
luciferase  activity  via  the  Dual-Luciferase  reporter  assay.  Values  are 
expressed  as  percentages  of  N-RIG-stimulated  luciferase  controls  set 
to  100.  Error  bars  indicate  standard  deviations  from  the  means  for 
triplicates.  (B)  The  cell  lysates  described  above  were  mixed  with  2x 
sample  buffer  and  subjected  to  12.5%  SDS-PAGE.  Following  transfer 
to  nitrocellulose,  the  membrane  was  blotted  with  mouse  anti-V5  to 
detect  the  proteases  and  antiactin  as  a  loading  control. 


viously  we  showed  that  the  biologically  relevant  transmem¬ 
brane-containing  form  of  SARS-CoV  PLpro,  termed  PLpro- 
TM,  is  a  potent  IFN  antagonist  (13).  To  determine  if  the 
membrane-tethered  version  of  HCoV-NL63  PLP2  can  func¬ 
tion  as  an  interferon  antagonist,  the  NL63  TM  sequence  was 
cloned  into  the  PLP2  construct  in  frame  with  the  V5  epitope 
tag,  and  the  resulting  construct  was  designated  PLP2-TM. 
HEK293  cells  were  transfected  with  PLP2-TM  or  the  catalytic 
mutant  PLP2-TM  C1678A  or  PLP2-TM  H1836A,  along  with 
the  IFN-p  and  pRL-TK  reporters.  N-RIG  was  used  to  stimu¬ 
late  IFN-p  induction.  We  found  that  PLP2-TM  and  the  cata¬ 
lytic  mutants  were  able  to  inhibit  N-RIG-induced  IFN-p  re¬ 
porter  activity  in  a  dose-dependent  manner,  although  like  the 
soluble  version  of  PLP2,  the  catalytic  mutants  were  less  effi¬ 
cient  than  the  wild  type  (Fig.  7A).  Expression  of  the  proteases 
was  verified  by  Western  blotting  (Fig.  7B).  These  data  corrob¬ 
orate  our  previous  results  indicating  that  the  catalytic  site  is  not 
essential  for  IFN-p  antagonism  by  HCoV-NL63  PLP2. 

PLpro  IFN  antagonism  is  not  dependent  on  the  ubiquitin- 
like  domain.  In  addition  to  the  downstream  TM  domain,  a 
previous  study  suggested  that  the  upstream  ubiquitin-like  do¬ 
main  (Ubl)  plays  a  role  in  the  IFN  antagonism  of  SARS-CoV 
PLpro,  since  deletion  of  this  domain  in  the  soluble  version  of 
PLpro  results  in  a  loss  of  IFN  antagonism  (16).  Currently,  no 
analogous  domain  has  been  identified  in  HCoV-NL63  PLP2. 
To  determine  the  role  of  the  Ubl  domain  in  the  more  biolog¬ 
ically  relevant  transmembrane  form  of  PLpro,  we  generated 
identical  Ubl  deletions  (AUbl)  of  PLpro  in  both  the  soluble 
and  transmembrane  forms  and  performed  a  dose-response 
profile  of  IFN  antagonism.  In  contrast  to  results  of  the  study  by 
Frieman  and  coworkers,  we  found  that  AUbl-PLpro-Sol  was  as 
potent  as  wild-type  PLpro  in  inhibiting  N-RIG-induced  IFN-p 
reporter  activity,  as  was  AUbl-PLpro-TM  (Fig.  8B).  Expression 
of  the  proteases  was  verified  by  Western  blotting  (Fig.  8C). 
These  results  indicate  that  the  Ubl  domain  of  PLpro  has  no 
effect  on  antagonism  of  type  I  IFN  induction. 
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FIG.  7.  The  transmembrane  form  of  HCoV-NL63  PLP2  and  cata¬ 
lytic  mutants  inhibit  RIG-I-mediated  IFN-(3  induction  in  a  dose-de- 
pendent  manner.  The  native  downstream  hydrophobic  domain  was 
cloned  into  the  PLP2  plasmid,  and  the  catalytic  cysteine  or  histidine 
residue  was  mutated  to  alanine.  (A)  HEK293  cells  were  transfected 
with  the  indicated  amounts  of  PLP2-TM,  PLP2-TM  C1678A,  or 
PLP2-TM  H1836A  along  with  the  IFN-(3-Luc  and  pRL-TK  reporters. 
N-RIG  was  used  to  stimulate  IFN-(3  induction.  At  24  h  posttransfec¬ 
tion,  cell  lysates  were  harvested  and  assayed  for  luciferase  activity  via 
the  Dual-Luciferase  reporter  assay.  Values  are  expressed  as  percent¬ 
ages  of  N-RIG-stimulated  luciferase  controls  set  to  100.  Error  bars 
indicate  standard  deviations  from  the  means  for  triplicates.  (B)  The 
cell  lysates  described  above  were  mixed  with  2X  sample  buffer  and 
subjected  to  12.5%  SDS-PAGE.  Following  transfer  to  nitrocellulose, 
the  membrane  was  blotted  with  mouse  anti-V5  to  detect  the  proteases 
and  antiactin  as  a  loading  control. 


Inhibition  of  SARS-CoV  PLpro  by  protease  inhibitor  GRL- 
0617S  has  no  effect  on  IFN-p  or  IFN-stimulated  response 
element  (ISRE)  reporter  activity  but  abrogates  inhibition  of 
NF-kB  reporter  activity.  The  mutagenesis  data  for  PLP2  sug¬ 
gest  that  the  catalytic  residues  (and  thus  catalytic  activity)  are 
not  required  for  interferon  antagonism,  and  Devaraj  et  al. 
reached  a  similar  conclusion  for  SARS-CoV  PLpro  (13).  To 
further  evaluate  the  role  of  protease/DUB  activity  in  inter¬ 
feron  antagonism  of  wild-type  PLpro,  we  added  a  protease 
inhibitor  (GRL-0617S)  that  has  been  developed  and  shown  to 
specifically  and  selectively  block  protease  and  DUB  activity  of 
SARS-CoV  PLpro  (56)  and  assessed  the  ability  of  PLpro  to 
inhibit  activation  of  IFN-p-Luc,  ISRE-Luc,  or  NF-kB-Luc  re¬ 
porter  activity.  FIEK293  cells  were  transfected  with  plasmid 
DNA  encoding  the  transmembrane  form  of  PLpro-TM  (amino 
acids  1541  to  2425  of  SARS-CoV  ORF  la),  previously  shown 
to  be  a  potent  IFN  antagonist  (13),  and  pRL-TK  along  with  the 
IFN-p,  ISRE,  or  NF-kB  reporter,  and  an  nsp2/3-GFP  sub¬ 
strate.  At  the  time  of  transfection,  the  cells  were  treated  with 
100  p-M  GRL-0617S  (10  times  the  50%  effective  concentration 
[EC50])  or  DMSO  (vehicle  control).  At  24  h  after  transfection, 
cell  lysates  were  generated  and  evaluated  for  induction  of  the 
reporters  and  proteolytic  processing  of  the  nsp2/3  substrate. 
We  found  that  inhibition  of  protease  activity  by  GRL-0617S 
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FIG.  8.  The  ubiquitin-like  domain  of  SARS-CoV  PLpro  is  not  re¬ 
quired  for  IFN  antagonism.  (A)  Schematic  diagram  of  nsp3  and  the 
various  V5-tagged  deletion  constructs.  Numbers  above  the  constructs 
indicate  the  amino  acid  numbers  counting  from  ORF  la.  (B)  HEK293 
cells  were  transfected  with  the  indicated  amounts  PLpro-Sol,  AUbl- 
PLpro-Sol,  or  AUbl-PLpro-TM  along  with  the  IFN-(3-Luc  and 
pRL-TK  reporters.  A  dominant  active  N-terminal  portion  of  RIG-I 
was  used  to  stimulate  IFN-(3  induction.  At  24  h  posttransfection,  cell 
lysates  were  harvested  and  assayed  for  luciferase  activity  via  the  Dual- 
Luciferase  reporter  assay.  Values  expressed  are  relative  to  results  for 
cells  transfected  with  the  reporters  alone.  Error  bars  indicate  standard 
deviations  from  the  means  for  triplicates.  (C)  The  cell  lysates  described 
above  were  mixed  with  2X  sample  buffer  and  subjected  to  12.5% 
SDS-PAGE.  Following  transfer  to  nitrocellulose,  the  membrane  was 
blotted  with  mouse  anti-V5  to  detect  the  proteases  and  antiactin  as  a 
loading  control. 


had  no  effect  on  the  IFN-p  reporter  (Fig.  9A)  and  little  to  no 
effect  on  the  ISRE  reporter,  which  is  dependent  on  IRF-3  (Fig. 
9B).  We  did  detect  an  alleviation  of  PLpro-mediated  inhibition 
of  the  NF-kB  reporter  by  GRL-0617S  (compare  Fig.  9A  and 
9B  with  9C).  These  results  indicate  that  protease/DUB  activity 
may  be  important  for  PLpro-mediated  inhibition  of  NF-kB 
activity  but  not  essential  for  inhibition  of  IRF-3  activity.  To 
demonstrate  efficacy  of  the  protease  inhibitor,  cell  lysates  used 
in  the  reporter  assay  were  assessed  for  nsp2/3-GFP  substrate 
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detected  in  cell  lysates  that  contain  PLpro  in  the  absence  of 
GRL-0617S.  In  contrast,  the  processing  of  nsp2/3-GFP  and 
liberation  of  nsp3-GFP  was  almost  completely  abrogated  in  the 
presence  of  the  drug.  Overall,  these  data  provide  further  sup¬ 
port  that  there  is  a  catalysis-independent  component  to  type  I 
IFN  antagonism  by  the  papain-like  proteases  of  human  coro- 
naviruses. 


DISCUSSION 
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FIG.  9.  SARS-CoV  PLpro  inhibits  IFN-0  and  ISRE  but  not  NF-kB 
reporter  activity  in  a  dose-dependent  manner  in  the  presence  or  ab¬ 
sence  of  a  protease  inhibitor.  HEK293  cells  were  transfected  with  the 
indicated  amounts  of  PLpro-TM,  pRL-TK,  nsp2/3-GFP,  and  either 
IFN-(3-Luc  (A),  ISRE-Luc  (B),  or  NF-kB-Luc  (C).  N-RIG  was  used  to 
stimulate  IFN-P  and  ISRE.  TNF-a  (10  ng/ml)  was  used  to  stimulate 
the  NF-kB-Luc  reporter.  DMSO  vehicle  control  or  100  p,M  GRL- 
0617S  was  added  at  the  time  of  transfection.  At  24  h  posttransfection, 
cell  lysates  were  harvested  and  assayed  for  luciferase  activity  via  the 
Dual-Luciferase  reporter  assay.  Values  are  expressed  as  percentages  of 
N-RIG-  or  TNF-a-stimulated  luciferase  controls,  set  to  100.  Error  bars 
indicate  standard  deviations  from  the  means  from  triplicates.  (D)  The 
cell  lysates  described  above  were  mixed  with  2X  sample  buffer  and 
subjected  to  12.5%  SDS-PAGE.  Following  transfer  to  nitrocellulose, 
the  membrane  was  blotted  with  mouse  anti-V5  to  detect  PLpro-TM 
and  rabbit  anti-GFP  to  detect  the  nsp2/3-GFP  substrate  and  the  nsp3- 
GFP  cleavage  product. 


cleavage.  The  nsp2/3-GFP  substrate  contains  a  region  of  nsp2/ 
nsp3  (including  the  cleavage  site)  fused  in-frame  with  GFP  and 
was  previously  shown  to  be  a  substrate  for  PLpro  (16).  As 
shown  in  Fig.  9D,  the  nsp3-GFP  cleavage  product  was  readily 


In  this  study,  we  describe  the  multifunctional  nature  of  the 
papain-like  protease  domains  of  NL63  and  SARS  coronavi- 
ruses.  These  coronavirus  PLP  domains  act  as  viral  proteases, 
deubiquitinating/delSGylating  enzymes,  and  are  able  to  antag¬ 
onize  innate  immune  induction  of  type  I  interferon.  We  found 
that  PLP  interferon  antagonism  is  enhanced  by,  but  is  not 
strictly  dependent  on,  the  catalytic  activity  of  the  enzyme.  In¬ 
hibition  of  coronavirus  protease  and  DUB  activity  by  mutagen¬ 
esis  or  pharmacological  means  did  not  abrogate  interferon 
antagonism.  Therefore,  these  distinct  PLP  activities  provide 
multiple  targets  for  antiviral  therapies. 

The  recognition  of  the  DUB/delSGylating  activity  of  SARS- 
CoV  PLpro  and  HCoV-NL63  PLP2  provides  new  opportuni¬ 
ties  to  investigate  how  the  virus  is  modifying  the  host  cell 
environment.  Posttranslational  modification  of  proteins  by 
ubiquitin  and  ubiquitin-like  (Ubl)  molecules,  such  as  SUMO, 
ISG15,  and  Nedd8,  plays  a  critical  role  in  the  regulatory  pro¬ 
cesses  of  virtually  all  aspects  of  cell  biology  (14,  15,  20,  22,  25, 
29,  34,  41,  55).  These  modifications,  though  covalent,  are 
highly  reversible.  Deubiquitinating  enzymes  can  deconjugate 
Ub  and  Ub-like  moieties  and  thus  modulate  the  activities  of 
ubiquitinated  proteins.  There  are  about  100  DUBs  encoded  in 
the  human  genome,  and  most  of  the  known  DUBs  are  cysteine 
proteases,  characterized  by  a  Cys-His-Asp  catalytic  triad  (51). 
Several  RNA  viruses  encode  cysteine  proteases  to  generate 
mature  viral  proteins  necessary  for  replication,  and  many  have 
been  found  to  be  multifunctional  proteins.  Like  SARS-CoV 
PLpro  and  HCoV-NL63  PLP2,  the  protease  of  a  nairovirus, 
Crimean  Congo  hemorrhagic  fever  virus  (CCHFV),  and  the 
proteases  of  arteriviruses,  including  equine  arteritis  virus 
(EAV)  and  porcine  respiratory  and  reproductive  syndrome 
virus  (PRRSV),  have  DUB  and  delSGylase  activity  (4,  15,  39). 
The  contribution  of  these  enzymatic  activities  to  inhibition  of 
type  I  IFN  induction  is  currently  poorly  understood.  A  vast 
array  of  proteins  involved  in  the  type  I  IFN  signaling  cascade 
are  activated  by  ubiquitination.  Induction  of  IFN-|3,  for  exam¬ 
ple,  requires  the  activation  of  IRF-3  and  NF-kB  (45,  55). 
Ubiquitination  is  known  to  be  intimately  involved  in  the  acti¬ 
vation  of  NF-kB.  Polyubiquitination  of  receptor-interacting 
protein  (RIP),  TNF  receptor-associated  factor  6  (TRAF6), 
and  TNF  receptor-associated  factor  2  (TRAF2)  activates  these 
signaling  intermediates,  which  leads  to  the  polyubiquitination 
of  IkB.  IkB,  which  binds  to  and  holds  NF-kB  inactive  in  the 
cytoplasm,  is  degraded  via  the  proteasome,  thereby  freeing 
NF-kB  to  translocate  to  the  nucleus  and  induce  IFN-(3  tran¬ 
scription  (20).  Here  we  have  shown  that  HCoV-NL63  PLP2 
has  a  profound  and  global  deconjugation  effect  on  ubiquiti¬ 
nated  cellular  conjugates,  suggesting  that  the  DUB  activity  of 
PLP2  may  modulate  NF-kB  activation.  In  addition,  we  dem¬ 
onstrate  that  protease  inhibitors,  which  block  coronavirus 
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DUB  activity,  abrogate  the  moderate  inhibition  of  NF-kB  re¬ 
porter  activity  imposed  by  transient,  ectopic  expression  of 
PLpro.  These  results  support  a  role  for  coronavirus  DUBs  in 
modulating  the  NF-kB  response  during  coronavirus  replica¬ 
tion.  However,  further  investigation  is  needed  to  delineate  the 
physiological  effect  of  coronavirus  PLPs  on  NF-kB  signaling, 
since  SARS-CoV  PLpro  did  not  inhibit  virus-  or  double- 
stranded  RNA  (dsRNA)-induced  activation  of  two  well-char¬ 
acterized  NF-KB-dependent  genes,  encoding  interleukin  6 
(IL-6)  and  A20,  when  it  was  stably  expressed  in  HeLa  cells  in 
a  tetracycline-regulated  fashion  (13). 

The  effect  of  ISGylation  of  cellular  proteins  on  the  antiviral 
response  is  far  less  understood.  It  is  known  that  ISG15 
conjugation  is  required  for  protection  against  lethal  Sindbis 
virus  infection  of  IFN-a/p  receptor  knockout  mice  (15). 
Also,  ISGylation  has  been  shown  to  influence  the  activation  of 
the  JAK/STAT  pathway,  involved  in  type  I  IFN  signaling  (44). 
Intriguingly,  IRF-3  also  undergoes  ISGylation  during  viral  in¬ 
fections,  which  was  found  to  enhance  innate  antiviral  re¬ 
sponses  by  inhibiting  virus-induced  IRF-3  degradation  (43). 
We  found  that  HCoV-NL63  PLP2  globally  deconjugates 
ISG15  similarly  to  its  deconjugation  of  Ub,  and  this  activity 
depends  on  the  catalytic  sites  of  PLP2.  Although  the  contribu¬ 
tion  of  the  delSGylation  activity  of  PLP2  to  IFN  antagonism 
remains  to  be  further  investigated,  delSGylation  of  IRF-3  is 
unlikely  to  be  a  significant  contributing  factor,  since  the  PLP 
catalytic  mutants  devoid  of  delSGylation  activity  still  effec¬ 
tively  inhibited  IFN  induction. 

Using  human  airway  epithelial  cells,  which  represent  a  cell 
culture  model  for  respiratory  infection,  we  found  that  IFN-p 
release  induced  by  HCoV-NL63  was  weak  but  measurable. 
This  finding  is  similar  to  weak  IFN  induction  by  the  far  more 
pathogenic  human  coronavirus  SARS-CoV,  suggesting  that 
antagonism  of  type  I  IFN  is  a  common  trait  of  coronavirus 
infection  (13,  58,  63).  In  addition  to  PLpro,  SARS-CoV  en¬ 
codes  several  other  IFN  antagonists.  ORF  3b,  ORF  6,  and 
nucleocapsid  inhibit  type  I  IFN  induction  via  inhibition  of 
IRF-3  phosphorylation  and  its  subsequent  nuclear  transloca¬ 
tion  (30).  ORF  6  and  nspl  have  been  shown  to  inhibit  IFN 
signaling  by  interfering  with  the  activity  of  STAT1  (17,  27,  48). 
Mouse  hepatitis  virus  also  encodes  several  IFN  antagonists, 
including  nspl,  nucleocapsid,  and  PLP2  (58,  68,  70).  Thus, 
there  is  clear  evolutionary  pressure  to  encode  and  maintain 
multiple  IFN  antagonists.  It  has  yet  to  be  determined  if  nspl 
and  nucleocapsid  from  HCoV-NL63  are  IFN  antagonists  as 
well. 

Using  reporter  assays,  we  found  that  HCoV-NL63  PLP2  can 
antagonize  type  I  IFN  induction  independently  of  catalytic 
activity.  Catalytic  mutants  of  PLP2  can  dose-dependently  in¬ 
hibit  IFN-p  induction;  however,  this  inhibition  is  reduced  com¬ 
pared  to  equivalent  amounts  of  wild-type  PLP2.  In  addition, 
we  note  that  the  presence  of  the  transmembrane  domain  con¬ 
fers  enhanced  IFN  antagonism,  particularly  in  the  catalytic 
mutants.  We  speculate  that  the  TM  domain  may  facilitate 
either  protein  folding  or  interaction  with  cellular  protein  part¬ 
ners.  Overall,  the  antagonism  profile  of  HCoV-NL63  PLP2  is 
remarkably  similar  to  that  of  PLpro  of  SARS-CoV.  In  fact, 
using  a  known  specific  inhibitor  of  SARS-CoV  PLpro  (GRL- 
0617S),  we  found  that  PLpro  can  effectively  inhibit  type  I  IFN 
induction  despite  a  profound  reduction  in  proteolysis,  which 


corroborates  the  notion  of  a  catalysis-independent  mechanism 
for  type  I  IFN  antagonism.  Previously  we  showed  that  SARS- 
CoV  PLpro  is  able  to  inhibit  the  phosphorylation  and  nuclear 
accumulation  of  IRF-3  following  Sendai  virus  infection  (13). 
We  found  that  HCoV-NL63  PLP2  was  also  able  to  inhibit  the 
translocation  of  IRF-3  to  the  nucleus  (data  not  shown);  how¬ 
ever,  the  mechanism  of  this  inhibition  is  not  yet  clear.  We  are 
actively  searching  for  cellular  factors  that  associate  with 
HCoV-NL63  PLP2  as  well  as  SARS-CoV  PLpro. 

The  crystal  structure  of  SARS-CoV  PLpro  has  identified  a 
unique  domain  that  has  remarkable  similarity  to  ubiquitin  (57). 
Frieman  and  coworkers  reported  that  removal  of  the  Ubl  do¬ 
main  in  the  soluble  version  of  PLpro  resulted  in  a  loss  of  IFN 
antagonism  (16).  However,  we  note  that  the  authors  of  that 
study  assessed  IFN  antagonism  using  one  concentration  of 
PLpro.  In  this  study,  by  performing  a  dose-response  profile  of 
IFN  antagonism,  we  found  that  deletion  of  the  Ubl  domain 
from  both  the  soluble  and  the  transmembrane  version  of 
PLpro  had  no  effect  on  IFN  antagonism.  Thus,  it  is  critical  to 
assess  the  effect  of  IFN  antagonism  across  a  range  of  protein 
concentrations  to  fully  evaluate  the  activity  of  these  proteases. 

Many  viruses  have  been  shown  to  inhibit  the  transcriptional 
activity  of  IRF-3  in  a  wide  variety  of  ways  (9,  24,  52).  Some 
viruses  inhibit  IRF-3  phosphorylation,  dimerization,  and/or 
translocation  to  the  nucleus.  Others  can  induce  IRF-3  degra¬ 
dation  or  sequester  the  transcription  factor  (11,  23).  The  mech¬ 
anism  of  inhibition  that  can  lead  to  these  phenotypes  can  occur 
directly  on  IRF-3  or  may  affect  any  of  the  vast  array  of  proteins 
upstream  of  IRF-3  in  the  type  I  IFN  induction  cascade.  For 
example,  the  VP35  protein  of  Ebola  Zaire  virus  (EBOV)  has 
been  shown  to  impact  IRF-3  activity  by  binding  dsRNA,  thus 
preventing  detection  by  RIG-I  (5,  6).  In  addition,  VP35  was 
shown  to  interact  with  Ubc9  (SUMO  E2  enzyme),  PIAS1 
(SUMO  E3  ligase),  and  IRF-7,  leading  to  SUMOylation  of 
IRF-7  and  transcriptional  repression  of  the  IFN-p  promoter 
(10).  Respiratory  syncytial  virus  (RSV)  encodes  two  proteins, 
NS1  and  NS2,  that  act  individually  or  cooperatively  to  inhibit 
the  activity  of  IRF-3.  It  was  reported  that  these  proteins  reduce 
the  expression  of  key  kinases  involved  in  IRF-3  phosphoryla¬ 
tion  (TRAF3  and  IkB  kinase  epsilon  [IKKe]),  but  how  NS1 
and  NS2  induce  TRAF3  and  IKKe  degradation  is  still  unclear 
(64).  Coronavirus  PLPs  could  be  acting  on  an  as  yet  unidenti¬ 
fied  cellular  factor  involved  in  the  IFN  induction  cascade.  For 
example,  in  a  recent  study  by  Schroder  et  al.,  a  new  protein, 
DEAD  box  protein  3  (DDX3),  was  found  to  be  involved  in  the 
type  I  IFN  cascade.  DDX3  was  identified  by  coimmunopre- 
cipitation  with  an  IFN  antagonist  from  vaccinia  virus  (VACV) 
protein  K7  (59).  Alternatively,  it  is  intriguing  to  speculate  that 
coronavirus  PLPs  may  function  by  sequestering  polyubiquitin 
complexes  or  membrane-associated  factors,  such  as  STING/ 
MITA  (26,  67,  72).  Previously  we  showed  that  the  block  must 
be  downstream  of  TBK1  but  at  or  upstream  of  IRF-3  since  a 
constitutively  active  form  of  IRF-3  is  not  blocked  by  PLpro 
(13).  Further  studies  are  needed  to  elucidate  the  mechanism  by 
which  coronavirus  PLPs  modulate  IRF-3  activity. 

The  multiple  enzymatic  activities  of  SARS-CoV  PLpro  and 
HCoV-NL63  PLP2  may  all  influence  the  host  cell  type  I  IFN 
response.  The  DUB  activity  of  these  proteins  could  modulate 
the  activity  of  key  players  in  the  signaling  cascade  that  are 
known  to  be  activated  by  lysine-48-  or  lysine-63-linked  poly- 
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ubiquitination.  Our  data  suggest  that  catalytic  activity  may 
contribute  to  IFN  antagonism  but  ablation  of  proteolysis  does 
not  abrogate  IFN  antagonism.  Thus,  coronavirus  PLPs  also 
possess  a  catalytic  activity-independent  mechanism  that  acts  to 
inhibit  IFN  induction.  We  are  currently  working  to  delineate 
both  the  DUB  and  catalytic  activity-independent  mechanism 
by  which  coronavirus  PLPs  inhibit  type  I  IFN  induction. 

The  data  presented  in  this  study  draw  significant  parallels 
between  the  single  papain-like  protease  of  SARS-CoV  PLpro 
and  the  second  papain-like  protease  of  HCoV-NL63  PLP2. 
Despite  modest  sequence  identity  (~19%),  these  two  pro¬ 
teases  have  similar  enzymatic  activities  and  can  inhibit  type  I 
IFN  induction  independently  of  catalytic  activity.  Since  coro¬ 
navirus  PLpro/PLP2  domains  are  required  for  viral  replication, 
they  are  attractive  targets  for  antiviral  therapeutics.  Indeed, 
inhibitors  of  SARS-CoV  PLpro  have  been  shown  to  block  virus 
replication  (56).  Though  less  pathogenic  than  SARS-CoV, 
HCoV-NL63  causes  significant  morbidity  in  children,  the  el¬ 
derly,  and  immune-compromised  individuals  and  has  been 
shown  to  be  an  etiological  agent  causing  croup  (65,  66).  In 
addition,  we  now  recognize  that  bats  and  other  mammals  can 
serve  as  reservoirs  for  potentially  emerging  pandemic  corona- 
viruses  (32,  37).  Thus,  further  studies  of  these  multifunctional 
coronavirus  PLPs  are  needed  to  determine  if  both  protease 
inhibitors  and  blockers  of  interferon  antagonism  can  be  devel¬ 
oped  to  reduce  replication  and  pathogenesis  of  human  and 
zoonotic  coronaviruses. 
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